Orientation Of Amino Acid Side Chains Research Articles

Evaluation of Myocilin Variant Protein Structures Modeled by AlphaFold2.

Deep neural network-based programs can be applied to protein structure modeling by inputting amino acid sequences. Here, we aimed to evaluate the AlphaFold2-modeled myocilin wild-type and variant protein structures and compare to the experimentally determined protein structures. Molecular dynamic and ligand binding properties of the experimentally determined and AlphaFold2-modeled protein structures were also analyzed. AlphaFold2-modeled myocilin variant protein structures showed high similarities in overall structure to the experimentally determined mutant protein structures, but the orientations and geometries of amino acid side chains were slightly different. The olfactomedin-like domain of the modeled missense variant protein structures showed fewer folding changes than the nonsense variant when compared to the predicted wild-type protein structure. Differences were also observed in molecular dynamics and ligand binding sites between the AlphaFold2-modeled and experimentally determined structures as well as between the wild-type and variant structures. In summary, the folding of the AlphaFold2-modeled MYOC variant protein structures could be similar to that determined by the experiments but with differences in amino acid side chain orientations and geometries. Careful comparisons with experimentally determined structures are needed before the applications of the in silico modeled variant protein structures.

Structures of brain-derived 42-residue amyloid-β fibril polymorphs with unusual molecular conformations and intermolecular interactions

Fibrils formed by the 42-residue amyloid-β peptide (Aβ42), a main component of amyloid deposits in Alzheimer's disease (AD), are known to be polymorphic, i.e., to contain multiple possible molecular structures. Previous studies of Aβ42 fibrils, including fibrils prepared entirely invitro or extracted from brain tissue and using solid-state NMR (ssNMR) or cryogenic electron microscopy (cryo-EM) methods, have found polymorphs with differences in amino acid sidechain orientations, lengths of structurally ordered segments, and contacts between cross-β subunit pairs within a single filament. Despite these differences, Aβ42 molecules adopt a common S-shaped conformation in all previously described high-resolution Aβ42 fibril structures. Here we report two cryo-EM-based structures of Aβ42 fibrils that are qualitatively different, in samples derived from AD brain tissue by seeded growth. In type A fibrils, residues 12 to 42 adopt a ν-shaped conformation, with both intra-subunit and intersubunit hydrophobic contacts to form a compact core. In type B fibrils, residues 2 to 42 adopt an υ-shaped conformation, with only intersubunit contacts and internal pores. Type A and type B fibrils have opposite helical handedness. Cryo-EM density maps and molecular dynamics simulations indicate intersubunit K16-A42 salt bridges in type B fibrils and partially occupied K28-A42 salt bridges in type A fibrils. The coexistence of two predominant polymorphs, with differences in N-terminal dynamics, is supported by ssNMR data, as is faithful propagation of structures from first-generation to second-generation brain-seeded Aβ42 fibril samples. These results demonstrate that Aβ42 fibrils can exhibit a greater range of structural variations than seen in previous studies.

Analysis of Multiple Protein Alignments Using 3D-Structural Information on the Orientation of Amino Acid Side-Chains

Multiple alignment of amino acid sequences of homologous proteins is a key tool in state-of-the-art bioinformatics and evolutionary analysis. Differences in the spatial orientation of amino acid side-chains can predetermine significant functional diversity among members of one superfamily; however, this is usually not taken into account in any way when constructing alignments and during subsequent comparative analysis. First of all, this is due to the limitation of existing algorithms, which are guided by the biochemical similarity of the "alphabet" of amino acid substitutions and either do not use information about the 3D-structural organization of proteins at all, or are limited to comparing the backbone only (i.e., the atoms of the main-chain). In this work, for the first time, we introduce new software for a systematic analysis of specific orientations of amino acid side-chains in equivalent positions of homologous protein structures. The program is intended to assist the analysis of protein multiple sequence alignments. The new algorithm, based on the machine learning HDBSCAN method, can identify statistically significant differences in the side-chain orientations and classify them into subfamilies at each position of multiple alignment. The method has been tested on a wide set of real biological data. The results allow us to speak of the specific orientation of amino acid side-chains as a common phenomenon that requires further study and deserves attention in a comparative analysis of functionally diverse protein superfamilies. The software is freely available at https://github.com/LimoninaDaria/Sub-family-Specific-Sidechain-Orientations.

Efficient consideration of coordinated water molecules improves computational protein-protein and protein-ligand docking discrimination.

Highly coordinated water molecules are frequently an integral part of protein-protein and protein-ligand interfaces. We introduce an updated energy model that efficiently captures the energetic effects of these ordered water molecules on the surfaces of proteins. A two-stage method is developed in which polar groups arranged in geometries suitable for water placement are first identified, then a modified Monte Carlo simulation allows highly coordinated waters to be placed on the surface of a protein while simultaneously sampling amino acid side chain orientations. This "semi-explicit" water model is implemented in Rosetta and is suitable for both structure prediction and protein design. We show that our new approach and energy model yield significant improvements in native structure recovery of protein-protein and protein-ligand docking discrimination tests.

Using chirality to probe the conformational dynamics and assembly of intrinsically disordered amyloid proteins

Intrinsically disordered protein (IDP) conformers occupy large regions of conformational space and display relatively flat energy surfaces. Amyloid-forming IDPs, unlike natively folded proteins, have folding trajectories that frequently involve movements up shallow energy gradients prior to the “downhill” folding leading to fibril formation. We suggest that structural perturbations caused by chiral inversions of amino acid side-chains may be especially valuable in elucidating these pathways of IDP folding. Chiral inversions are subtle in that they do not change side-chain size, flexibility, hydropathy, charge, or polarizability. They allow focus to be placed solely on the question of how changes in amino acid side-chain orientation, and the resultant alterations in peptide backbone structure, affect a peptide’s conformational landscape (Ramachandran space). If specific inversions affect folding and assembly, then the sites involved likely are important in mediating these processes. We suggest here a “focused chiral mutant library” approach for the unbiased study of amyloid-forming IDPs.

GPCRs through the keyhole: the role of protein flexibility in ligand binding to β-adrenoceptors

G protein-coupled receptors (GPCRs) are proteins of pharmaceutical importance, with over 30% of all drugs in clinical use targeting them. Increasing numbers of X-ray crystal (XRC) structures of GPCRs offer a wealth of data relating to ligand binding. For the β-adrenoceptors (β-ARs), XRC structures are available for human β2- and turkey β1-subtypes, in complexes with a range of ligands. While these structures provide insight into the origins of ligand structure-activity relationships (SARs), questions remain. The ligands in all published complexed XRC structures lack extensive substitution, with no obvious way the ligand-binding site can accommodate β1-AR-selective antagonists with extended side-chains para- to the common aryloxypropanolamine pharmacophore. Using standard computational docking tools with such ligands generally returns poses that fail to explain known SARs. Application of our Active Site Pressurisation modelling method to β-AR XRC structures and homology models, however, reveals a dynamic area in the ligand-binding pocket that, through minor changes in amino acid side chain orientations, opens a fissure between transmembrane helices H4 and H5, exposing intra-membrane space. This fissure, which we term the “keyhole”, is ideally located to accommodate extended moieties present in many high-affinity β1-AR-selective ligands, allowing the rest of the ligand structure to adopt a canonical pose in the orthosteric binding site. We propose the keyhole may be a feature of both β1- and β2-ARs, but that subtle structural differences exist between the two, contributing to subtype-selectivity. This has consequences for the rational design of future generations of subtype-selective ligands for these therapeutically important targets.

Idealized Models of Protofilaments of Human Islet Amyloid Polypeptide

Fibrils formed by assembly of human islet amyloid polypeptide (hIAPP) are found in most patients with type II diabetes. Structurally, these fibrils are composed of multiple protofilaments and are characterized by extended beta sheets, variable helical twists, and different morphologies. We have previously derived models for the hIAPP protofilament using simulations constrained by data from EPR spectroscopy. In the current work, these models were used as a basis for generating idealized hIAPP protofilaments with symmetrical geometrical properties using a new algorithm, MFIBRIL. We show good agreement of the idealized protofilaments with experimental data for amino acid side chain orientations and geometrical features including the inter-β sheet distance and the protofilament radius. These idealized protofilaments can be used in MFIBRIL to generate fibril models that may be experimentally testable at the molecular level. MFIBRIL can also be used for building structures of any repetitive molecular assembly starting with a single building block obtained from any source.

Ca2+ Binding Effects on the C2 Domain Conformation of Human Cytosolic Phospholipase A2

It has been reported that the cooperative binding of calcium ions indicated a local conformational change of the human cytosolic phospholipase A2 (cPLA2) C2 domain (Nalefski et al., (1997) Biochemistry 36, 12011-12018). However its structural evidence is less known (Malmberg et al., (2003) Biochemistry 42, 13227-13240). In this letter, life-time decay and fluorescence quenching techniques were employed to compare the calcium-induced conformational changes. The life-time decay parameters and fluorescence quenching constant changes were small between the apo- and holo-C2 domains when tryptophan residue was excited at 295 nm. In contrast, the quenching constant change was large, from 0.52 M(-1) for the apo-C2 to 8.8 M(-1) for the holo-C2 domain, when tyrosine residues were excited at 284 nm. Our results provide new information on amino acid side chain orientation change at calcium binding loop 3, which is necessary for Ca2+ binding regulated membrane targeting of human cytosolic phospholipase A2.

Direct observation of protonation reactions during the catalytic cycle of cytochrome c oxidase.

Cytochrome c oxidase, the terminal protein in the respiratory chain, converts oxygen into water and helps generate the electrochemical gradient used in the synthesis of ATP. The catalytic action of cytochrome c oxidase involves electron transfer, proton transfer, and O2 reduction. These events trigger specific molecular changes at the active site, which, in turn, influence changes throughout the protein, including alterations of amino acid side chain orientations, hydrogen bond patterns, and protonation states. We have used IR difference spectroscopy to investigate such modulations for the functional intermediate states E, R2,Pm, and F. These spectra reveal deprotonation of its key glutamic acid E286 in the E and in the Pm states. The consecutive deprotonation and reprotonation of E286 twice within one catalytic turnover illustrates the role of this residue as a proton shuttle. In addition, the spectra point toward deprotonation of a redox-active tyrosine, plausibly Y288, in the F intermediate. Structural insights into the molecular mechanism of catalysis based on the subtle molecular changes observed with IR difference spectroscopy are discussed.

Protalign: a 3-dimensional protein alignment assessment tool.

Protein fold recognition (sometimes called threading) is the prediction of a protein's 3-dimensional shape based on its similarity to a protein of known structure. Fold predictions are low resolution; that is, no effort is made to rotate the protein's component amino acid side chains into their correct spatial orientations. The goal is simply to recognize the protein family member that most closely resembles the target sequence of unknown structure and to create a sensible alignment of the target to the known structure (i.e., a structure-sequence alignment). To facilitate this type of structure prediction, we have designed a low resolution molecular graphics tool. ProtAlign introduces the ability to interact with and edit alignments directly in the 3-dimensional structure as well as in the usual 2-dimensional layout. It also contains several functions and features to help the user assess areas within the alignment. ProtAlign implements an open pipe architecture to allow other programs to access its molecular graphics capabilities. In addition, it is capable of "driving" other programs. Because amino acid side chain orientation is not relevant in fold recognition, we represent amino acid residues as abstract shapes or glyphs much like Lego (tm) blocks and we borrow techniques from comparative flow visualization using streamlines to provide clean depictions of the entire protein model. By creating a low resolution representation of protein structure, we are able to at least double the amount of information on the screen. At the same time, we create a view that is not as busy as the corresponding representations using traditional high resolution visualization methods which show detailed atomic structure. This eliminates distracting and possibly misleading visual clutter resulting from the mapping of protein alignment information onto a high resolution display of the known structure. This molecular graphics program is implemented in Open GL to facilitate porting to other platforms.

The synthesis of 1,6-disubstituted indanes which mimic the orientation of amino acid side-chains in a protein alpha-helix motif.

We utilize a 1,6-disubstituted indane as a template onto which two amino acid side-chains are appended in an orientation which mimics that found in a protein alpha-helix motif.

1H-NMR assignments and the dynamics of interconversion of the isomeric forms of cytochrome b5 in solution

Cytochrome b 5 reconstituted with specifically deuterated hemins has led to the assignment of the resolved 6,7β-propionate protons and heme meso protons. Freshly reconstituted cytochrome b 5 contains a mixture of two isomers in an approx. 1:1 ratio. As time proceeds the minor isomer decreases in intensity until the equilibrium ratio, approx. 8:1, of the two isomers is reached. The rate of the heme disorder kinetics was investigated for cytochrome b 5 as a function of pH, oxidation state and 2,4 heme substitutents. Comparison of the kinetic data for cytochrome b 5 with that obtained for other b-type heme proteins supports the proposal that the heme disorder arises from a 180° rotation of the heme about the α,γ- meso axis. Computer-difference methods allow the spectra of the two individual isomers to be generated. Comparison of the NMR spectral parameters for the two individual isomers indicates small structural differences for amino acid side-chain orientations.

Evidence for conformational differences in aqueous and crystalline structures of α-bungarotoxin and cobratoxin

Laser Raman spectroscopy has been employed to investigate the structures of α-bungarotoxin ( Bungarus multicinctus) and cobratoxin ( Naja naja siamensis) in H 2O and D 2O solutions. Structures of the aqueous neurotoxins are compared with one another and with the X-ray crystal structures. The results indicate that the solution and crystal molecular structures of cobratoxin are in substantial agreement with one another, but those of α-bungarotoxin are not. Raman data provide no evidence for strained disulfides in aqueous α-bungarotoxin, although strained CSSC dihedral angles are indicated for the X-ray crystal structure. The data are interpreted as evidence for a strained molecular conformation of α-bungarotoxin in the crystal, which converts to a relaxed, more energetically favorable conformation in aqueous solution. Raman spectra also suggest more β-strand secondary structure in aqueous α-bungarotoxin (47 ± 5%) than in the crystalline form ( < 10%). The high β-strand content measured by Raman spectroscopy could be due to either a secondary structure in solution that is appreciably different than that of the crystal, or to the imprecision of the Raman method in distinguishing peptide configurations that are vibrationally equivalent but conformationally inequivalent. Aqueous α-bungarotoxin and cobratoxin also differ from one another in amino acid side chain orientations and interactions, though not in main chain conformations. Different geometries are indicated for cystine CCSS dihedral angles, and different hydrogen bonding states are indicated for internal tyrosines. Tyrosine-24 of α-bungarotoxin is shown to donate a strong hydrogen bond to a negative acceptor, deduced to be glutamate-41, whereas the equivalently positioned residue of cobratoxin is apparently hydrogen bonded to solvent molecules.

Stereochemical control of ribosomal peptidyltransferase reaction. Role of amino acid side-chain orientation of acceptor substrate.

The substrate specificity of the acceptor site of peptidyltransferase of Escherichia coli 70S ribosomes was investigated in the fMet-tRNA.A-U-G.70S ribosome and AcPhe-tRNA.poly(U).70S ribosome systems by using a series of 2'- and 3'-aminoacyldinucleoside phosphates as acceptors. These chemically synthesized compounds are analogues of the 3' termini of either 2'(3')-, 2'-, or 3'-aminoacyl transfer ribonucleic acids (AA-tRNAs) of the types C-A-aa, C-2'-dA-aa, C-3'-dA-aa, C-3'-dA-3'-NH-aa, and C-2'-dA-2'-NH-aa (aa = Phe, D-Phe, Lys, Leu, Ala, Glu, Pro, Gly, Asp, Met, and alpha-aminoisobutyryl). It was found that the 3'-aminoacyl derivatives of optically active amino acids are much better acceptors of N-formyl-L-methionine (fMet) or N-acetyl-L-phenylalanine (AcPhe) residues than the isomeric 2'-aminoacyl derivatives with affinity constant ratios (KM 2'/3') greater than 100. Likewise, C-A(D-Phe) is a weaker acceptor than the corresponding L derivative C-A-Phe. In contrast, all glycyl derivatives (C-2'-dA-Gly, C-3'-dA-Gly, C-3'-dA-3'-NH-Gly and C-2'-dA-2'-NH-Gly) are good acceptors of the fMet residue, with ratios (KM 2'/3') of approximately 2. On the basis of these results, a model for the stereochemical control of the peptidyl-transferase reaction is proposed. It assigns a major role to the orientation of the amino acid side chain in 2'- or 3'-AA-tRNA. A detailed model of the interaction of the acceptor terminus of 3'-AA-tRNA with the acceptor site of peptidyltransferase is also proposed. The model is strikingly similar to those for the active sites of proteolytic enzymes.

Stereochemical control of the ribosomal peptidyltransferase reaction: The role of acceptor substrate amino acid side chain orientation

The aminoacyl residue of aa-tRNA is bound by an ester linkage to the 2’or 3’-hydroxyl groups of the ribose residue of the 3’-terminal adenosine [l] and is able to migrate rapidly between these two positions [2]. Thus, the determination of an ‘active’ aa-tRNA isomer in various subreactions of protein synthesis has been a topic of considerable interest. Modified aa-tRNAs, or their 3’-terminal fragments, in which the 3%’ transacylation cannot occur, have been recently used for studies of isomer-specificity of the acceptor site of peptidyltransferase [3]. Most of the studies have concluded that 3’-PhetRNA, is the only acceptor in the peptide bondforming step [3], although a significant acceptor activity was also observed [4] with 2’-Phe-tRNA. The comparison of reaction rates of both isomers was not performed, but it was thought that the 3’Phe-tRNA is probably the preferred substrate [5]. From these experiments important questions have evolved: 1. Is the peptidyltransferase really specific for 3’-aatRNA as acceptor substrate and is this specificity