Side-chain Chi Research Articles

Everyone Is a Protagonist: Residue Conformational Preferences in High-Resolution Protein Structures.

In many structural bioinformatics problems, there is a broad range of unanswered questions about protein dynamics and amino acid properties. Proteins are not strictly static objects, but rather populate ensembles of conformations. One way to understand these particularities is to analyze the information available in experimental databases. The Ramachandran plot, despite being more than half a century old, remains an utterly useful tool in the study of protein conformation. Based on its assumptions, we inspected a large data set (11,130 protein structures, amounting to 5,255,768 residues) and discriminated the conformational preferences of each residue type regarding their secondary structure participation. These data were studied for phi [Formula: see text], psi [Formula: see text], and side chain chi [Formula: see text] angles, being presented in non-Ramachandranian plots. In the largest analysis of protein conformation made so far, we propose an original plot to depict conformational preferences in relation to different secondary structure elements. Despite confirming previous observations, our results strongly support a unique character for each residue type, whereas also reinforcing the observation that side chains have a major contribution to secondary structure and, by consequence, on protein conformation. This information can be further used in the development of more robust methods and computational strategies for structural bioinformatics problems.

Sidechain chi1 distribution in B3 domain of protein G from extensive sets of residual dipolar couplings

Sidechain chi1 distribution in B3 domain of protein G from extensive sets of residual dipolar couplings

Identification of natural inhibitors against angiotensin I converting enzyme for cardiac safety using induced fit docking and MM-GBSA studies.

Background:Cleistanthins A and B are isolated compounds from the leaves of Cleistanthus collinus Roxb (Euphorbiaceae). This plant is poisonous in nature which causes cardiovascular abnormalities such as hypotension, nonspecific ST-T changes and QTc prolongation. The biological activity predictions spectra of the compounds show the presence of antihypertensive, diuretic and antitumor activities.Objective:Objective of the present study was to determine the in silico molecular interaction of cleistanthins A and B with Angiotensin I- Converting Enzyme (ACE-I) using Induced Fit Docking (IFD) protocols.Materials and Methods:All the molecular modeling calculations like IFD docking, binding free energy calculation and ADME/Tox were carried out using Glide software (Schrödinger LLC 2009, USA) in CentOS EL-5 workstation.Results:The IFD complexes showed favorable docking score, glide energy, glide emodel, hydrogen bond and hydrophobic interactions between the active site residues of ACE-I and the compounds. Binding free energy was calculated for the IFD complexes using Prime MM-GBSA method. The conformational changes induced by the inhibitor at the active site of ACE-I were observed based on changes of the back bone Cα atoms and side-chain chi (x) angles. The various physicochemical properties were calculated for these compounds. Both cleistanthins A and B showed better docking score, glide energy and glide emodel when compared to captopril inhibitor.Conclusion:These compounds have successively satisfied all the in silico parameters and seem to be potent inhibitors of ACE-I and potential candidates for hypertension.

Designing of Protein Kinase C β-II Inhibitors against Diabetic complications: Structure Based Drug Design, Induced Fit docking and analysis of active site conformational changes.

Protein Kinase C β-II (PKC β-II) is an important enzyme in the development of diabetic complications like cardiomyopathy, retinopathy, neuropathy, nephropathy and angiopathy. PKC β-II is activated in vascular tissues during diabetic vascular abnormalities. Thus, PKC β-II is considered as a potent drug target and the crystal structure of the kinase domain of PKC β-II (PDB id: 2I0E) was used to design inhibitors using Structure-Based Drug Design (SBDD) approach. Sixty inhibitors structurally similar to Staurosporine were retrieved from PubChem Compound database and High Throughput Virtual screening (HTVs) was carried out with PKC β-II. Based on the HTVs results and the nature of active site residues of PKC β-II, Staurosporine inhibitors were designed using SBDD. Induced Fit Docking (IFD) studies were carried out between kinase domain of PKC β-II and the designed inhibitors. These IFD complexes showed favorable docking score, glide energy, glide emodel and hydrogen bond and hydrophobic interactions with the active site of PKC β-II. Binding free energy was calculated for IFD complexes using Prime MM-GBSA method. The conformational changes induced by the inhibitor at the active site of PKC β-II were observed for the back bone Cα atoms and side-chain chi angles. PASS prediction tool was used to analyze the biological activities for the designed inhibitors. The various physicochemical properties were calculated for the compounds. One of the designed inhibitors successively satisfied all the in silico parameters among the others and seems to be a potent inhibitor against PKC β-II.

KING (Kinemage, Next Generation): A versatile interactive molecular and scientific visualization program

Proper visualization of scientific data is important for understanding spatial relationships. Particularly in the field of structural biology, where researchers seek to gain an understanding of the structure and function of biological macromolecules, it is important to have access to visualization programs which are fast, flexible, and customizable. We present KiNG, a Java program for visualizing scientific data, with a focus on macromolecular visualization. KiNG uses the kinemage graphics format, which is tuned for macromolecular structures, but is also ideal for many other kinds of spatially embedded information. KiNG is written in cross-platform, open-source Java code, and can be extended by end users through simple or elaborate "plug-in" modules. Here, we present three such applications of KiNG to problems in structural biology (protein backbone rebuilding), bioinformatics of high-dimensional data (e.g., protein sidechain chi angles), and classroom education (molecular illustration). KiNG is a mature platform for rapidly creating and capitalizing on scientific visualizations. As a research tool, it is invaluable as a test bed for new methods of visualizing scientific data and information. It is also a powerful presentation tool, whether for structure browsing, teaching, direct 3D display on the web, or as a method for creating pictures and videos for publications. KiNG is freely available for download at http://kinemage.biochem.duke.edu.

Solution structure of a hydrophobic analogue of the winter flounder antifreeze protein.

The solution structure of a synthetic mutant type I antifreeze protein (AFP I) was determined in aqueous solution at pH 7.0 using nuclear magnetic resonance (NMR) spectroscopy. The mutations comprised the replacement of the four Thr residues by Val and the introduction of two additional Lys-Glu salt bridges. The antifreeze activity of this mutant peptide, VVVV2KE, has been previously shown to be similar to that of the wild type protein, HPLC6 (defined here as TTTT). The solution structure reveals an alphahelix bent in the same direction as the more bent conformer of the published crystal structure of TTTT, while the side chain chi1 rotamers of VVVV2KE are similar to those of the straighter conformer in the crystal of TTTT. The Val side chains of VVVV2KE assume the same orientations as the Thr side chains of TTTT, confirming the conservative nature of this mutation. The combined data suggest that AFP I undergoes an equilibrium between straight and bent helices in solution, combined with independent equilibria between different side chain rotamers for some of the amino acid residues. The present study presents the first complete sequence-specific resonance assignments and the first complete solution structure determination by NMR of any AFP I protein.

Efficient measurement of (3)J(N,Cgamma) and (3)J(C',Cgamma) coupling constants of aromatic residues in (13)C, (15)N-labeled proteins.

An NMR pulse sequence is proposed for the simultaneous determination of side chain chi1 torsion-angle related (3)J(N,Cgamma) and (3)J(C', Cgamma) couplings in aromatic amino acid spin systems. The method is of the quantitative J correlation type and takes advantage of attenuated (15)N and (1)H transverse relaxation by means of the TROSY principle. Unlike previously developed schemes for the measurement of either of the two coupling types, spectra contain internal reference peaks that are usually recorded in separate experiments. Therefore, the desired information is extracted from a single rather than four data sets. The new method is demonstrated with uniformly (13)C/(15)N labeled Desulfovibrio vulgaris flavodoxin, which contains 14 aromatic out of 147 total amino acid residues.

Heteronuclear relayed E.COSY revisited: determination of 3J(H(alpha),C(gamma)) couplings in Asx and aromatic residues in proteins.

Constant-time 3D heteronuclear relayed E.COSY [Schmidt et al. (1996) J. Biomol. NMR, 7, 142-152], as based on generic 2D small-flip-angle HMQC-COSY [Schmidt et al. (1995) J. Biomol. NMR, 6, 95-105], has been modified to allow for quantitative determination of heteronuclear three-bond 3J(H(alpha),C(gamma)) couplings. The method is applicable to amino acid spin topologies with carbons in the gamma position which lack attached protons, i.e. to asparagine, aspartate, and aromatic residues in uniformly 13C-enriched proteins. The pulse sequence critically exploits heteronuclear triple-quantum coherence (HTQC) of CH2 moieties involving geminal H(beta) proton pairs, taking advantage of improved multiple-quantum relaxation properties, at the same time avoiding scalar couplings between those spins involved in multiple-quantum coherence, thus yielding E.COSY-type multiplets with a splitting structure that is simpler than with the original scheme. Numerical least-squares 2D line-shape simulation is used to extract 3J(H(alpha),C(gamma)) coupling constants which are of relevance to side-chain chi1 dihedral-angle conformations in polypeptides. Methods are demonstrated with recombinant 15N,13C-enriched ribonuclease T1 and Desulfovibrio vulgaris flavodoxin with bound oxidized FMN.

Conformational mimicry: synthesis and solution conformation of a cyclic somatostatin hexapeptide containing a tetrazole cis amide bond surrogate.

Potent, cyclic hexapeptide analogues of somatostatin are generally believed to adopt some common secondary structural features: a II' beta turn at one end of the cycle, and a type VI turn with a cis amide bond at the other. A proposed cis amide surrogate, the 1,5-disubstituted tetrazole, has been placed into a cyclic hexapeptide analog of somatostatin in order to constrain the putative cis amide bond. The final cyclization was done by either chemical or enzymatic means. The product, cyclo(Ala6-Tyr7-D-Trp8-Lys9-Val10-Phe11-psi[CN4] ), was found to have 83% of the activity of somatostatin. Solution nmr analysis in DMSO/water revealed that the backbone as well as side chain chi1 and chi2 were well ordered. Relaxation matrix methods were used to extract distance restraints from the nuclear Overhauser effect spectroscopy data set, and these were used in a systematic search of torsional space to identify structures consistent with the nmr data. Restrained minimizations of these structures using a number of different force fields produced structures having the expected beta II' turn at D-Trp8-Lys9 and a beta VIa turn in the Phe11-psi[CN4]-Ala6 portion of the molecule. The similarity of the minimized structures to those previously reported for cyclic hexapeptide analogues of somatostatin confirms the similarity of the tetrazole geometry to that of the cis amide in solution.

Three-dimensional structure of bovine heart fatty-acid-binding protein with bound palmitic acid, determined by multidimensional NMR spectroscopy.

The three-dimensional structure of the holo form of recombinant cellular bovine heart fatty-acid-binding protein (H-FABPc), a polypeptide of 133 amino acid residues with a molecular mass of 15 kDa, has been determined by multidimensional homonuclear and heteronuclear NMR spectroscopy applied to uniformly 15N-labeled and unlabeled protein. A nearly complete set of 1H and 15N chemical shift assignments was obtained. A total of 2329 intramolecular distance constraints and 42 side-chain chi 1 dihedral-angle constraints were derived from cross-relaxation and J coupling information. 3D nuclear Overhauser enhancement and exchange spectroscopy combined with heteronuclear multiple-quantum coherence (NOESY-HMQC) experiments, performed on a sample of uniformly 13C-labeled palmitic acid bound to unlabeled cellular heart fatty-acid-binding protein revealed 10 intermolecular contacts that determine the orientation of the bound fatty acid. An ensemble of protein conformations was calculated with the distance-geometry algorithm for NMR applications (DIANA) using the redundant dihedral-angle constraint (REDAC) strategy. After docking the fatty acid into the protein, the protein-ligand arrangement was subject to distance-restrained energy minimization. The overall conformation of the protein is a beta-barrel consisting of 10 antiparallel beta-strands which form two nearly orthogonal beta-sheets of five strands each. Two short helices form a helix-turn-helix motif in the N-terminal region of the polypeptide chain. The palmitic acid is bound within the protein in a U-shaped conformation close to the two helices. The obtained solution structure of the protein is consistent with a number of fatty-acid-binding-protein crystal structures.

Three-Dimensional Structure of Bovine Heart Fatty-acid-binding Protein with Bound Palmitic Acid, Determined by Multidimensional NMR Spectroscopy

The three-dimensional structure of the holo form of recombinant cellular bovine heart fatty-acid-binding protein (H-FABPc), a polypeptide of 133 amino acid residues with a molecular mass of 15 kDa, has been determined by multidimensional homonuclear and heteronuclear NMR spectroscopy applied to uniformly 15N-labeled and unlabeled protein. A nearly complete set of 1H and 15N chemical shift assignments was obtained. A total of 2329 intramolecular distance constraints and 42 side-chain chi 1 dihedral-angle constraints were derived from cross-relaxation and J coupling information. 3D nuclear Overhauser enhancement and exchange spectroscopy combined with heteronuclear multiple-quantum coherence (NOESY-HMQC) experiments, performed on a sample of uniformly 13C-labeled palmitic acid bound to unlabeled cellular heart fatty-acid-binding protein revealed 10 intermolecular contacts that determine the orientation of the bound fatty acid. An ensemble of protein conformations was calculated with the distance-geometry algorithm for NMR applications (DIANA) using the redundant dihedral-angle constraint (REDAC) strategy. After docking the fatty acid into the protein, the protein-ligand arrangement was subject to distance-restrained energy minimization. The overall conformation of the protein is a beta-barrel consisting of 10 antiparallel beta-strands which form two nearly orthogonal beta-sheets of five strands each. Two short helices form a helix-turn-helix motif in the N-terminal region of the polypeptide chain. The palmitic acid is bound within the protein in a U-shaped conformation close to the two helices. The obtained solution structure of the protein is consistent with a number of fatty-acid-binding-protein crystal structures.

Three-dimensional solution structure of Cucurbita maxima trypsin inhibitor-V determined by NMR spectroscopy.

The solution structure of Cucurbita maxima trypsin inhibitor-V (CMTI-V), which is also a specific inhibitor of the blood coagulation protein, factor XIIa, was determined by 1H NMR spectroscopy in combination with a distance-geometry and simulated annealing algorithm. Sequence-specific resonance assignments were made for all the main-chain and most of the side-chain hydrogens. Stereospecific assignments were also made for some of the beta-, gamma-, delta-, and epsilon-hydrogens and valine methyl hydrogens. The ring conformations of all six prolines in the inhibitor were determined on the basis of 1H-1H vicinal coupling constant patterns; most of the proline ring hydrogens were stereospecifically assigned on the basis of vicinal coupling constant and intraresidue nuclear Overhauser effect (NOE) patterns. Distance constraints were determined on the basis of NOEs between pairs of hydrogens. Dihedral angle constraints were determined from estimates of scalar coupling constants and intraresidue NOEs. On the basis of 727 interproton distance and 111 torsion angle constraints, which included backbone phi angles and side-chain chi 1, chi 2, chi 3, and chi 4 angles, 22 structures were calculated by a distance geometry algorithm and refined by energy minimization and simulated annealing methods. Both main-chain and side-chain atoms are well-defined, except for a loop region, two terminal residues, and some side-chain atoms located on the molecular surface. The average root mean squared deviation in the position for equivalent atoms between the 22 individual structures and the mean structure obtained by averaging their coordinates is 0.58 +/- 0.06 A for the main-chain atoms and 1.01 +/- 0.07 A for all the non-hydrogen atoms of residues 3-40 and 49-67. These structures were compared to the X-ray crystallographic structure of another protein of the same inhibitor family-chymotrypsin inhibitor-2 from barley seeds [CI-2; McPhalen, C. A., & James, M. N. G. (1987) Biochemistry 26, 261-269]. The main-chain folding patterns are highly similar for the two proteins, which possess 62% sequence differences. However, major differences are noted in the N- and C-terminal segments, which may be due to the presence of a disulfide bridge in CMTI-V, but not in CI-2.

Conformation and sweet tastes of L-aspartyl dipeptide methyl esters.

In order to investigate the conformational preferences to elicit tastes, conformational free energy calculations using an empirical potential (ECEPP/2) and the hydration shell model were carried out on the L-aspartyl dipeptide methyl esters, L-(+)HAsp(-)-L-Xaa-OMe, in the hydrated state, where Xaa includes sweet (Phe, Tyr, Met, and Gly), bitter (Ala, Trp, Val, Leu, and Ile), and tasteless (Ser, Thr, and Abu) residues. The refined preferred conformation of the Phe dipeptide (aspartame) with side chain chi 1/2 conformation g- is g-Fg- in the hydrated state, which is consistent with the structure deduced from 1H-nmr experiments. Irrespective of the Xaa and taste, all the dipeptides have the same conformation for the Asp residue, which is attributable to the hydrogen bond between protonated amino hydrogen and carboxylate oxygen and the favored hydration of the carboxylate group. This implies that the L-aspartyl residue is a necessary factor for the dipeptides to be sweet not a sufficient factor. The computed conformational preferences for sweet, bitter, and tasteless dipeptides in the hydrated state indicate to us that the conformation about the N--C alpha bond of the Xaa residue, i.e., the orientation of the hydrophobic moiety with respect to the AH/B functionalities in the aspartyl moiety, seems to be crucial to elicit the tastes. In addition, the hydrophobicity and the size of the Xaa residue are found to play a major role in determining the tastes. These well accord with the related works reported previously.

Fast and simple Monte Carlo algorithm for side chain optimization in proteins: application to model building by homology.

An unknown protein structure can be predicted with fair accuracy once an evolutionary connection at the sequence level has been made to a protein of known 3-D structure. In model building by homology, one typically starts with a backbone framework, rebuilds new loop regions, and replaces nonconserved side chains. Here, we use an extremely efficient Monte Carlo algorithm in rotamer space with simulated annealing and simple potential energy functions to optimize the packing of side chains on given backbone models. Optimized models are generated within minutes on a workstation, with reasonable accuracy (average of 81% side chain chi 1 dihedral angles correct in the cores of proteins determined at better than 2.5 A resolution). As expected, the quality of the models decreases with decreasing accuracy of backbone coordinates. If the back-bone was taken from a homologous rather than the same protein, about 70% side chain chi 1 angles were modeled correctly in the core in a case of strong homology and about 60% in a case of medium homology. The algorithm can be used in automated, fast, and reproducible model building by homology.

Conformational studies of stereoisomeric 14-membered cyclic enkephalin analogues containing 1-naphthylalanine at the fourth position: chirality effect of leucine at the fifth position on biological activity and receptor selectivity.

In order to study structure-activity relationships of enkephalin-related analogues, we report the biological activity and conformational analysis of four 14-membered cyclic enkephalin analogues with beta-(1-naphthyl) alanine in place of phenylalanine at the fourth position, Tyr-c[D-A2bu-Gly-(L and D)-beta Nal(1)-(L and D)-Leu]. The L-beta Nal(1)-containing analogues display higher activity at both the mu and delta receptors than the corresponding analogues with the L-Phe residue. In contrast to the linear enkephalins, the cyclic analogues with the D-beta Nal(1) residue are also active at the mu receptor since the relative spatial arrangement of functional groups required for biological activity is achieved by the constrained nature of the cyclic molecules. A comparison of the findings from the conformational analysis and biological assays establishes that relatively extended structures, in which the two aromatic side chains are oriented in opposite directions with a approximately 14 A separation, is required for activity at the mu receptor. On the other hand, folded conformations with nearly parallel orientations and a close proximity (less than 10A) of the aromatic rings of the Tyr and beta Nal(1) residues are required for activity at the delta receptor. It should be noted that the overall structures and thus the biological profiles of the 14-membered cyclic enkephalin analogues are strongly dependent on the conformation of the second residue. The folded conformations with parallel orientation of the two aromatic side chains of Tyr-c[D-A2bu-Gly-L-beta Nal(1)-D-Leu] is stabilized by an interaction between the Tyr phenolic OH proton and beta Nal(1) C*O groups. This analogue, which shows the highest activity at both the mu and delta receptors among the four stereoisomers studied, displays an increase of the fraction of the side-chain chi 1 = t conformer for the beta Nal(1) residue. It is concluded that the incorporation of the D-Leu residue at the fifth position increases the relative fraction of the folded conformations with parallel orientation of the aromatic side chains, and hence enhances activity at the delta receptor as compared to the corresponding L-Leu containing analogue.