On residues in the disallowed region of the Ramachandran map.

An analysis of the occurrence of nonglycyl residues in conformations disallowed in the Ramachandran plot is presented. Ser, Asn, Thr, and Cys have the highest propensities to exhibit such conformations, and the branched aliphatic residues the lowest. Residues cluster in five regions and there are some trends in the types of residues and their side-chain conformations (chi(1)) occupying these. Majority of the residues are found at the edge of helices and strands and in short loops, and are involved in different types of weak, stabilizing interactions. A structural motif has been identified where a residue in disallowed conformation occurs as the first residue of a short 3(10)-helix. On the basis of the types of neighboring residues, the location in the three-dimensional structure and accessibility, there are similarities with the occurrence of cis peptide bonds in protein structures.     

Crystal structures at 2.2 A resolution of the catalytic domains of normal ras protein and an oncogenic mutant complexed with GDP

The biological functions of ras proteins are controlled by the bound guanine nucleotide GDP or GTP. The GTP-bound conformation is biologically active, and is rapidly deactivated to the GDP-bound conformation through interaction with GAP (GTPase Activating Protein). Most transforming mutants of ras proteins have drastically reduced GTP hydrolysis rates even in the presence of GAP. The crystal structures of the GDP complexes of ras proteins at 2.2 A resolution reveal the detailed interaction between the ras proteins and the GDP molecule. All the currently known transforming mutation positions are clustered around the bound guanine nucleotide molecule. The presumed "effector" region and the GAP recognition region are both highly exposed. No significant structural differences were found between the GDP complexes of normal ras protein and the oncogenic mutant with valine at position 12, except the side-chain of the valine residue. However, comparison with GTP-analog complexes of ras proteins suggests that the valine side-chain may inhibit GTP hydrolysis in two possible ways: (1) interacting directly with the gamma-phosphate and altering its orientation or the conformation of protein residues around the phosphates; and/or (2) preventing either the departure of gamma-phosphate on GTP hydrolysis or the entrance of a nucleophilic group to attack the gamma-phosphate. The structural similarity between ras protein and the bacterial elongation factor Tu suggests that their common structural motif might be conserved for other guanine nucleotide binding proteins.     

Ramachandran plot on the web (2.0)

The Ramachandran plot displays the main chain conformation angles (Phi and Psi) of the polypeptide chain of a protein molecule. The paper reports the updated version of the Ramachandran plot web server and has several improved options for displaying the conformation angles in various regions. In addition, options are provided to display the conformation angles in various secondary structural elements and regions within the user specified Phi and Psi values in the plot. The updated version is accessible at the following URL: http://dicsoft1.physics.iisc.ernet.in/rp/.     

CADB: Conformation Angles DataBase of proteins

Conformation Angles DataBase (CADB) provides an online resource to access data on conformation angles (both main-chain and side-chain) of protein structures in two data sets corresponding to 25% and 90% sequence identity between any two proteins, available in the Protein Data Bank. In addition, the database contains the necessary crystallographic parameters. The package has several flexible options and display facilities to visualize the main-chain and side-chain conformation angles for a particular amino acid residue. The package can also be used to study the interrelationship between the main-chain and side-chain conformation angles. A web based JAVA graphics interface has been deployed to display the user interested information on the client machine. The database is being updated at regular intervals and can be accessed over the World Wide Web interface at the following URL: http://144.16.71.148/cadb/.     

Assessing side-chain perturbations of the protein backbone: a knowledge-based classification of residue Ramachandran space

Grouping the 20 residues is a classic strategy to discover ordered patterns and insights about the fundamental nature of proteins, their structure, and how they fold. Usually, this categorization is based on the biophysical and/or structural properties of a residue's side-chain group. We extend this approach to understand the effects of side chains on backbone conformation and to perform a knowledge-based classification of amino acids by comparing their backbone phi, psi distributions in different types of secondary structure. At this finer, more specific resolution, torsion angle data are often sparse and discontinuous (especially for nonhelical classes) even though a comprehensive set of protein structures is used. To ensure the precision of Ramachandran plot comparisons, we applied a rigorous Bayesian density estimation method that produces continuous estimates of the backbone phi, psi distributions. Based on this statistical modeling, a robust hierarchical clustering was performed using a divergence score to measure the similarity between plots. There were seven general groups based on the clusters from the complete Ramachandran data: nonpolar/beta-branched (Ile and Val), AsX (Asn and Asp), long (Met, Gln, Arg, Glu, Lys, and Leu), aromatic (Phe, Tyr, His, and Cys), small (Ala and Ser), bulky (Thr and Trp), and, lastly, the singletons of Gly and Pro. At the level of secondary structure (helix, sheet, turn, and coil), these groups remain somewhat consistent, although there are a few significant variations. Besides the expected uniqueness of the Gly and Pro distributions, the nonpolar/beta-branched and AsX clusters were very consistent across all types of secondary structure. Effectively, this consistency across the secondary structure classes implies that side-chain steric effects strongly influence a residue's backbone torsion angle conformation. These results help to explain the plasticity of amino acid substitutions on protein structure and should help in protein design and structure evaluation.     

Conformational properties of trans Ac-Asn-Pro-Tyr-NHMe and trans Ac-Tyr-Pro-Asn-NHMe in dimethylsulfoxide and in water determined by multinuclear n.m.r. spectroscopy

Vicinal coupling constants between various nuclei provide backbone and side-chain conformational information for a series of asparagine- and tyrosine-containing peptides in DMSO and in H2O. By enriching Tyr of Ac-Asn-Pro-Tyr-NHMe with 15N, it has been possible to distinguish between the resonances of the two side-chain beta protons of Tyr. Analysis of the coupling constants in terms of the distributions of side-chain conformations in these peptides indicates that the addition of Asn to the Pro-Tyr sequence leads to a less random conformational distribution. When compared to the side-chain rotamer distribution of Ac-Asn-NHMe and Ac-Tyr-NHMe, particular Asn and Tyr side-chain conformations of Ac-Asn-Pro-Tyr-NHMe are stabilized in dimethylsulfoxide solution. The interaction(s) which stabilize a unique Tyr side-chain conformation of Ac-Asn-Pro-Tyr-NHMe in dimethylsulfoxide are not present in Ac-Ala-Pro-Tyr-NHMe and are unaffected by the addition of Val-Pro to the C-terminus of Asn-Pro-Tyr. In water, a preferential stabilization of one Asn side-chain conformation of Ac-Asn-Pro-Tyr-NHMe is also observed, while the Tyr side-chain rotamer distribution is similar to that of Ac-Tyr-NHMe. An interaction between the Asn side chain and the Pro-Tyr-NHMe backbone was previously shown to stabilize a beta-bend conformation at Pro-Tyr in water. Data are also presented for Ac-Tyr-Pro-Asn-NHMe, for which local interactions do not stabilize particular backbone conformations in dimethylsulfoxide or in water. The conformations of the peptides studied here are relatively insensitive to temperatures between 27 degrees and 62 degrees, both in dimethylsulfoxide and in water. The sequences Asn-Pro-Tyr and Tyr-Pro-Asn occur in ribonuclease A, and these tripeptides serve as models for the interactions involved in the folding of this protein.     

Genetic algorithm with alternating selection pressure for protein side-chain packing and pK(a) prediction

The prediction of protein side-chain conformation is central for understanding protein functions. Side-chain packing is a sub-problem of protein folding and its computational complexity has been shown to be NP-hard. We investigated the capabilities of a hybrid (genetic algorithm/simulated annealing) technique for side-chain packing and for the generation of an ensemble of low energy side-chain conformations. Our method first relies on obtaining a near-optimal low energy protein conformation by optimizing its amino-acid side-chains. Upon convergence, the genetic algorithm is allowed to undergo forward and “backward” evolution by alternating selection pressures between minimal and higher energy setpoints. We show that this technique is very efficient for obtaining distributions of solutions centered at any desired energy from the minimum. We outline the general concepts of our evolutionary sampling methodology using three different alternating selective pressure schemes. Quality of the method was assessed by using it for protein pK(a) prediction.     

A data bank merging related protein structures and sequences.

A data collection which merges protein structural and sequence information is described. Structural superpositions amongst proteins with similar main-chain fold were performed or collected from the literature. Sequences taken from the protein primary structure databases were associated with the multiple structural alignments providing they were at least 50% homologous in residue identity to one of the structural sequences and at least 50% of the structural sequence residues were alignable. Such restrictions allow reasonable confidence that the primary sequences share the conformation of the tertiary structural templates, except in the less conserved loop regions. Multiple structural superpositions were collected for 38 familial groups containing a total of 209 tertiary structures; 45 structures had no superposable mates and were used individually. Other information is also provided as main-chain and side-chain conformational angles, secondary structural assignments and the like. Wedding the primary and tertiary structural data resulted in an 8-fold increase of data bank sequence entries over those associated with the known three-dimensional architectures alone.     

Tensile Mechanics of Alanine-Based Helical Polypeptide: Force Spectroscopy versus Computer Simulations

In nature, an α-helix is commonly used to build thermodynamically stable and mechanically rigid protein conformations. In view of growing interest in the mechanical rigidity of proteins, we measured the tensile profile of an alanine-based α-helical polypeptide on an atomic-force microscope to investigate the basic mechanics of helix extension with minimal interference from side-chain interactions. The peptide was extended to its maximum contour length with much less force than in reported cases of poly-L-Glu or poly-L-Lys, indicating that chain stiffness strongly depended on the physicochemical properties of side chains, such as their bulkiness. The low tensile-force extension originated presumably in locally unfolded parts because of spontaneous structural fluctuations. In 50% trifluoroethanol, the well-known helix-promoting agent, the rigidity of the sample polypeptide was markedly increased. Computer simulations of the peptide-stretching process showed that a majority of constituent residues underwent a transition from an α-helical to an extended conformation by overcoming an energy barrier around ψ ∼0° on the Ramachandran plot. The observed lability of an isolated helix signified the biological importance of the lateral bundling of helices to maintain a rigid protein structure.     

Structure of myoglobin-ethyl isocyanide. Histidine as a swinging door for ligand entry

The structure of myoglobin(Fe II)-ethyl isocyanide has been solved at 1.68 A resolution by X-ray crystallography. The isocyano group of the ligand is distorted from the linear conformation observed in solution and in model compounds. Local changes in the protein conformation are also seen. The side-chain of Arg-CD3 moves out into the solvent, and the side-chain of His-E7 swings up and away from the ligand. Both of these side-chains show disorder indicative of dynamic behavior. These outward movements of His-E7 and Arg-CD3 side-chains clear a path from the solvent to the heme iron, suggesting a mechanism for ligand entry.     

One nanosecond molecular dynamics simulation of the N-terminal domain of the lambda repressor protein

We have carried out molecular dynamics simulation of the N-terminal domain of the lambda repressor protein in a surrounding environment including explicit waters and ions. We observe two apparent dynamics substates in the nanosecond protein simulation, the transition occurring around 500 ps. The existence of these two apparent substates results from a high flexibility of the arm in each monomer, a relative flexibility of both arms with respect to each other, and a relative displacement of the recognition helices from 30 to 40 A of interhelical distance. Many amino acid residues, including those involved in DNA recognition, undergo a simultaneous transition in their side-chain conformations, consistent with the relationship between side-chain conformation and secondary structural elements, as observed in protein crystal structures. This result suggests plausible conformational changes experienced by the protein upon DNA binding. On the whole, the non-consensus monomer appears to be more flexible than its consensus counterpart.     

Optimization of van der Waals energy for protein side-chain placement and design

Computational determination of optimal side-chain conformations in protein structures has been a long-standing and challenging problem. Solving this problem is important for many applications including homology modeling, protein docking, and for placing small molecule ligands on protein-binding sites. Programs available as of this writing are very fast and reasonably accurate, as measured by deviations of side-chain dihedral angles; however, often due to multiple atomic clashes, they produce structures with high positive energies. This is problematic in applications where the energy values are important, for example when placing small molecules in docking applications; the relatively small binding energy of the small molecule is drowned by the large energy due to atomic clashes that hampers finding the lowest energy state of the docked ligand. To address this we have developed an algorithm for generating a set of side-chain conformations that is dense enough that at least one of its members would have a root mean-square deviation of no more than R Å from any possible side-chain conformation of the amino acid. We call such a set a side-chain cover set of order R for the amino acid. The size of the set is constrained by the energy of the interaction of the side chain to the backbone atoms. Then, side-chain cover sets are used to optimize the conformation of the side chains given the coordinates of the backbone of a protein. The method we use is based on a variety of dead-end elimination methods and the recently discovered dynamic programming algorithm for this problem. This was implemented in a computer program called Octopus where we use side-chain cover sets with very small values for R, such as 0.1 Å, which ensures that for each amino-acid side chain the set contains a conformation with a root mean-square deviation of, at most, R from the optimal conformation. The side-chain dihedral-angle accuracy of the program is comparable to other implementations; however, it has the important advantage that the structures produced by the program have negative energies that are very close to the energies of the crystal structure for all tested proteins.     

13C-NMR chemical shifts and the conformations of rigid polypeptides

¹³C-nmr chemical shifts of backbone carbonyl and side-chain β-carbons in polypeptides provide structural information. Recent utilization of substituent effects on ¹³C-nmr chemical shifts (principally γ-effects) has permitted the rationalization of their sequence and conformation dependence when observed in linear, flexible polypeptides. In this report, we apply the γ-effect method to study the ¹³C-nmr chemical shifts observed in solution and in the solid state for the backbone carbonyl and side-chain β-carbons in conformationally rigid polypeptides, which are usually cyclic. As found previously for flexible, linear polypetides, the relative ¹³C-nmr chemical shifts observed for the backbone carbonyl and side-chain β-carbons in conformationally rigid polypeptides are predictable from knowledge of their peptide residue sequence (primary structure) and conformation (secondary structure) via the γ-effect method.     

DNA binding mediates conformational changes and metal ion coordination in the active site of PcrA helicase.

Based upon the crystal structures of PcrA helicase, we have made and characterised mutations in a number of conserved helicase signature motifs around the ATPase active site. We have also determined structures of complexes of wild-type PcrA with ADPNP and of a mutant PcrA complexed with ADPNP and Mn2+. The kinetic and structural data define roles for a number of different residues in and around the ATP binding site. More importantly, our results also show that there are two functionally distinct conformations of ATP in the active site. In one conformation, ATP is hydrolysed poorly whereas in the other (activated) conformation, ATP is hydrolysed much more rapidly. We propose a mechanism to explain how the stimulation of ATPase activity afforded by binding of single-stranded DNA stabilises the activated conformation favouring Mg2+binding and a consequent repositioning of the gamma-phosphate group which promotes ATP hydrolysis. A part of the associated conformational change in the protein forces the side-chain of K37 to vacate the Mg2+binding site, allowing the cation to bind and interact with ATP.     

Trans-substitution of the proximal hydrogen bond in myoglobin: I. Structural consequences of hydrogen bond deletion

The structural role of a side-chain to side-chain protein hydrogen bond is examined using trans-substitution of the proximal histidine of myoglobin with methylimidazoles (Barrick, Biochemistry 1994;33:6546-6554). Modification of the chemical structure of exogenous ligands allows this hydrogen bond to be disrupted. Comparison of the crystal structures of H93G myoglobin complexed 4-methylimidazole (4meimd; methylation at carbon 4) and 1-methylimidazole (1meimd; methylation at the adjacent nitrogen, preventing hydrogen bonding between the imidazole ligand and the protein) shows that the polypeptide, heme, and methylimidazole orientations are the same within error. For 4meimd there appear to be major and minor conformations corresponding to different tautomeric states of the ligand. Conformational heterogeneity is also seen in the hyperfine-shifted region of the NMR spectrum of 4meimd complexed with high-spin H93G deoxyMb. The major conformation of the 4meimd ligand and the 1meimd ligand, as seen in the respective crystal structures, are quite similar except that the proximal ligand NH-to-Ser92-OH hydrogen bond is eliminated in the 1meimd complex, and instead the proximal ligand CH is adjacent to the Ser92-OH. Thus, this system provides a means to eliminate the Mb proximal hydrogen bond in a chemically and structurally conservative way.     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

Secondary and tertiary structures of the transmembrane domains of the translocator protein TSPO determined by NMR. Stabilization of the TSPO tertiary fold upon ligand binding

Numerous biological functions are attributed to the peripheral-type benzodiazepine receptor (PBR) recently renamed translocator protein (TSPO). The best characterized function is the translocation of cholesterol from the outer to inner mitochondrial membrane, which is a rate-determining step in steroid biosynthesis. TSPO drug ligands have been shown to stimulate pregnenolone formation by inducing TSPO-mediated translocation of cholesterol. Until recently, no direct structural data on this membrane protein was available. In a previous paper, we showed that a part of the mouse TSPO (mTSPO) C-terminal region adopts a helical conformation, the side-chain distribution of which provides a groove able to fit a cholesterol molecule. We report here on the overall structural properties of mTSPO. This study was first undertaken by dissecting the protein sequence and studying the conformation of five peptides encompassing the five putative transmembrane domains from ¹H-NMR data. The secondary structure of the recombinant protein in micelles was then studied using CD spectroscopy. In parallel, the stability of its tertiary fold was probed using ¹H-¹⁵N NMR. This study provides the first experimental evidence for a five-helix fold of mTSPO and shows that the helical conformation of each transmembrane domain is mainly formed through local short-range interactions. Our data show that, in micelles, mTSPO exhibits helix content close to what is expected but an unstable tertiary fold. They reveal that the binding of a drug ligand that stimulates cholesterol translocation is able to stabilize the mTSPO tertiary structure.     

Accuracy of side-chain prediction upon near-native protein backbones generated by ab initio folding methods

The ab initio folding problem can be divided into two sequential tasks of approximately equal computational complexity: the generation of native-like backbone folds and the positioning of side chains upon these backbones. The prediction of side-chain conformation in this context is challenging, because at best only the near-native global fold of the protein is known. To test the effect of displacements in the protein backbones on side-chain prediction for folds generated ab initio, sets of near-native backbones (≤ 4 Å Cα RMS error) for four small proteins were generated by two methods. The steric environment surrounding each residue was probed by placing the side chains in the native conformation on each of these decoys, followed by torsion-space optimization to remove steric clashes on a rigid backbone. We observe that on average 40% of the χ1 angles were displaced by 40° or more, effectively setting the limits in accuracy for side-chain modeling under these conditions. Three different algorithms were subsequently used for prediction of side-chain conformation. The average prediction accuracy for the three methods was remarkably similar: 49% to 51% of the χ1 angles were predicted correctly overall (33% to 36% of the χ1+2 angles). Interestingly, when the inter-side-chain interactions were disregarded, the mean accuracy increased. A consensus approach is described, in which side-chain conformations are defined based on the most frequently predicted χ angles for a given method upon each set of near-native backbones. We find that consensus modeling, which de facto includes backbone flexibility, improves side-chain prediction: χ1 accuracy improved to 51–54% (36–42% of χ1+2). Implications of a consensus method for ab initio protein structure prediction are discussed. Proteins 33:204–217, 1998. © 1998 Wiley-Liss, Inc.     

Predicting structural effects in HIV-1 protease mutant complexes with flexible ligand docking and protein side-chain optimization

We present a computational approach for predicting structures of ligand-protein complexes and analyzing binding energy landscapes that combines Monte Carlo simulated annealing technique to determine the ligand bound conformation with the dead-end elimination algorithm for side-chain optimization of the protein active site residues. Flexible ligand docking and optimization of mobile protein side-chains have been performed to predict structural effects in the V32I/I47V/V82I HIV-1 protease mutant bound with the SB203386 ligand and in the V82A HIV-1 protease mutant bound with the A77003 ligand. The computational structure predictions are consistent with the crystal structures of these ligand-protein complexes. The emerging relationships between ligand docking and side-chain optimization of the active site residues are rationalized based on the analysis of the ligand-protein binding energy landscape. Proteins 33:295–310, 1998. © 1998 Wiley-Liss, Inc.