Proton nuclear magnetic resonance study of the solution distal histidine orientation in monomeric Chironomus thummi thummi cyanomet hemoglobins. Dynamic stability of the heme pocket as monitored by labile proton exchange.

The 1H nuclear magnetic resonance spectral characteristics of the cyano-Met form of Chironomus thummi thummi monomeric hemoglobins I, III and IV in 1H2O solvent are reported. A set of four exchangeable hyperfine-shifted resonances is found for each of the two heme-insertion isomers in the hyperfine-shifted region downfield of ten parts per million. An analysis of relaxation, exchange rates and nuclear Overhauser effects leads to assignments for all these resonances to histidine F8 and the side-chains of histidine E7 and arginine FG3. It is evident that in aqueous solution, the side-chain from histidine E7 does not occupy two orientations, as found for the solid state, rather the histidine E7 side-chain adopts a conformation similar to that of sperm whale myoglobin or hemoglobin A, oriented into the heme pocket and in contact with the bound ligand. Evidence is presented to show that the allosteric transition in the Chironomus thummi thummi hemoglobins arises from the "trans effect". An analysis of the exchange with bulk solvent of the assigned histidine E7 labile proton confirms that the group is completely buried within the heme pocket in a manner similar to that found for sperm whale cyano-Met myoglobin, and that the transient exposure to solvent is no more likely than in mammalian myoglobins with the "normal" distal histidine orientation. Finally, a comparison of solvent access to the heme pocket of the three monomeric C. thummi thummi hemoglobins, as measured from proton exchange rates of heme pocket protons, is made and correlated to binding studies with the diffusible small molecules such as O2.     

Analysis of the relationship between side-chain conformation and secondary structure in globular proteins

The relationship between the preferred side-chain dihedral angles and the secondary structure of a residue was examined. The structures of 61 proteins solved to a resolution of 2.0 A (1 A = 0.1 nm) or better were analysed using a relational database to store the information. The strongest feature observed was that the chi 1 distribution for most side-chains in an alpha-helix showed an absence of the g- conformation and a shift towards the t conformation when compared to the non-alpha/beta structures. The exceptions to this tendency were for short polar side-chains that form hydrogen bonds with the main-chain which prefer g+. Shifts in the chi 1 preferences for residues in the beta-sheet were observed. Other side-chain dihedral angles (chi 2, chi 3, chi 4) were found to be influenced by the main-chain. This paper presents more accurate distributions for the side-chain dihedral angles which were obtained from the increased number of proteins determined to high resolution. The means and standard deviations for chi 1 and chi 2 angles are presented for all residues according to the secondary structure of the main-chain. The means and standard deviations are given for the most popular conformations for side-chains in which chi 3 and chi 4 rotations affect the position of C atoms.     

X-ray crystal structure of the ferric sperm whale myoglobin: imidazole complex at 2.0 A resolution

The X-ray crystal structure of the ferric sperm whale (Physeter catodon) myoglobin:imidazole complex has been refined at 2.0 A resolution, to a final R-factor of 14.8%. The overall conformation of the protein is little affected by binding of the ligand. Imidazole is co-ordinated to the heme iron at the distal site, and forces distinguishable local changes in the surrounding protein residues. His64(E7) swings out of the distal pocket and becomes substantially exposed to the solvent: nevertheless, it stabilizes the exogenous ligand by hydrogen bonding. The side-chains of residues Arg45(CD3) and Asp60(E3) are also affected by imidazole association.     

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.     

On the role of the crystal environment in determining protein side-chain conformations

The role of crystal packing in determining the observed conformations of amino acid side-chains in protein crystals is investigated by (1) analysis of a database of proteins that have been crystallized in different unit cells (space group or unit cell dimensions) and (2) theoretical predictions of side-chain conformations with the crystal environment explicitly represented. Both of these approaches indicate that the crystal environment plays an important role in determining the conformations of polar side-chains on the surfaces of proteins. Inclusion of the crystal environment permits a more sensitive measurement of the achievable accuracy of side-chain prediction programs, when validating against structures obtained by X-ray crystallography. Our side-chain prediction program uses an all-atom force field and a Generalized Born model of solvation and is thus capable of modeling simple packing effects (i.e. van der Waals interactions), electrostatic effects, and desolvation, which are all important mechanisms by which the crystal environment impacts observed side-chain conformations. Our results are also relevant to the understanding of changes in side-chain conformation that may result from ligand docking and protein-protein association, insofar as the results reveal how side-chain conformations change in response to their local environment.     

Study of the binding environment of alpha-factor in its G protein-coupled receptor using fluorescence spectroscopy.

Mating in Saccharomyces cerevisiae is induced by the interaction of alpha-factor (W1H2W3L4Q5L6K7P8G9Q10P11M12Y13) with its cognate G protein-coupled receptor (Ste2p). Fifteen fluorescently labeled analogs of alpha-factor in which the 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD) group was placed at the alphaN-terminus and in side-chains at positions 1, 3, 4, 6, 7, 12 and 13 were synthesized and assayed for biological activity and receptor affinity. Eleven of the analogs retained 6-60% of the biological activity of the alpha-factor, as judged using a growth arrest assay. The binding affinities depended on the position of NBD attachment in the peptide and the distance of the tag from the backbone. Derivatization of the positions 3 and 7 side-chains with the NBD group resulted in analogs with affinities of 17-35% compared with that of alpha-factor. None of the other NBD-containing agonists had sufficient receptor affinity or strong enough emission for fluorescence analysis. The position 3 and 7 analogs were investigated using fluorescence spectroscopy and collisional quenching by KI in the presence of Ste2p in yeast membranes. The results showed that the lambda max of NBD in the position 7 side-chain shifted markedly to the blue (510 nm) when separated by 4 or 6 bonds from the peptide backbone and that this probe was shielded from quenching by KI. In contrast, separation by 3, 5, 10 or more bonds resulted in lambda max ( approximately 540 nm) and collisional quenching constants consistent with increasing degrees of exposure. The NBD group in the position 3 side-chain was also found to be blue shifted (lambda max=520 nm) and shielded from solvent. These results indicate that the position 7 side-chain is likely interacting with a pocket formed by extracellular domains of Ste2p, whereas the side-chain of Trp3 is in a hydrophobic pocket possibly within the transmembrane region of the receptor.     

What contributions to protein side-chain dynamics are probed by NMR experiments? A molecular dynamics simulation analysis

Molecular dynamics simulations of the structurally homologous proteins TNfn3 and FNfn10 have been used to investigate the contributions to side-chain dynamics measured by NMR relaxation experiments. The results reproduce the variation in core side-chain dynamics observed by NMR and highlight the relevance of anharmonic motion and transitions between local minima for explaining NMR side-chain order parameters. A method is described for calculating converged order parameters by use of replica exchange molecular dynamics in conjunction with an implicit solvent model. These simulations allow the influence of various factors, such as the flexibility of side-chains and their free volume, on the mobility to be tested by perturbing the system. Deletion mutations are found to have the largest effect on the more densely packed FNfn10. Some counterintuitive effects are seen, such as an increase in order parameters close to deletion mutation sites, but these can be rationalized in terms of direct interactions with the modified side-chains. A statistical analysis of published order parameters supports the conclusions drawn from the simulations.     

Determination of preferred conformations of brassinosteroids by means of NMR investigations and Boltzmann statistical analysis of simulated annealing calculations

Brassinosteroids are a class of steroidal phytohormones with high growth-promoting properties. The preferred side-chain conformations of 10 brassinosteroids were determined by means of detailed NMR investigations and molecular-modeling studies. Vacuum conformations obtained by simulated annealing calculations and Boltzmann statistical analysis were compared with solution conformations derived from NOE experiments and molecular dynamic simulations, and with X-ray structures. In general, results from simulated annealing calculations and NMR-supported molecular dynamics simulations are in good agreement. For some of the compounds investigated the conformation was less well-defined at the end of the side-chain. It could be shown that the energetically most favorable and most probable conformations also include the conformations obtained by NMR supported molecular-dynamics calculations and by X-ray analysis. For the most bioactive compound brassinolide (1) the majority of conformations show a side-chain bent towards the -face of the steroid skeleton, whereas for the less bioactive brassinosteroids, conformations with straight side-chains or side-chains bent towards the -face are more frequent.Electronic Supplementary Material available.     

Hydrophilicity of polar amino acid side-chains is markedly reduced by flanking peptide bonds

Several amino acid side-chain hydropathy scales have been devised on the basis of solubility and water/organic solvent partitioning data obtained with free amino acids or side-chain analogs. In nearly all cases, these scales are based upon the structure-additivity assumption; it has been assumed that the transfer free energies of the amino acid side-chains are the same in these model compounds as they are in a polypeptide. This assumption is probably wrong. In the present study, deviations from additivity for amino acid side-chains are demonstrated by comparing a theoretically derived scale, which N-acetylamino acid amides. The results show that the flanking peptide bonds dramatically reduce the hydrophilicity of the polar side-chains, with deviations up to several kilocalories (1 kcal = 4.184 kJ) for the charged side-chains at pH 7.0. Further calculation shows that these deviations are due to reductions of 40 to 85% in the unfavorable transfer free energy of the polar functional groups. In addition, proximity of the neighboring amide bonds in the parent molecule (N-acetylglycine amide) decreases the hydrophilicity of the -CONH-backbone unit by 36%. This decrease is expected to be twice as large for -CONH- units in the interior of a polypeptide backbone. The significance of these observations is: (1) valid hydropathy scales can be obtained only with model peptides; (2) deviations from additivity are expected in all solvent systems, including non-polar solvents that are thought to mimic the interior of a membrane; (3) the spontaneous insertion of polypeptides into membranes is likely to occur much more readily than has been previously thought. In order to estimate the free energy of transferring the side-chains and the polypeptide backbone from water to the interior of a lipid bilayer, the results of this study are used to construct a hydropathy scale based upon the partitioning of solutes between water and non-polar solvents. The validity of hydropathy scales that are based on criteria other than solubility and water/organic solvent partitioning data is also discussed.     

Refinement of the structure of recombinant rat intestinal fatty acid-binding apoprotein at 1.2-A resolution.

The three-dimensional structure of the 131-residue rat intestinal fatty acid-binding protein, without bound ligand (apoI-FABP), has been refined with x-ray diffraction data to a nominal resolution of 1.19 A. The final model has a conventional crystallographic R-factor of 16.9% for 34,290 unique reflections [a root mean square (r.m.s.) deviation for bond length of 0.012 A and a r.m.s. deviation of 2.368 degrees for bond angles]. Ninety-two residues are present as components of the protein's 10 anti-parallel beta-strands while 14 residues are part of its two short alpha-helices. The beta-strands and alpha-helices are organized into two nearly orthogonal beta-sheets. Particular attention has been placed in defining solvent structure and the structures of discretely disordered groups in this protein. Two hundred thirty-seven solvent molecules have been identified; 24 are located within apoI-FABP. The refined model includes alternate conformers for 228 protein atoms (109 main-chain, 119 side-chain) and 63 solvent molecules. We have found several aromatic side-chains with multiple conformations located near, or in, the protein's ligand binding site. This observation, along with the fact that these side-chains have a temperature factor that is relatively higher than that of other aromatic residues, suggests that they may be involved in the process of noncovalent binding of fatty acid. The absence of a true hydrophobic core in I-FABP suggests that its structural integrity may be maintained primarily by a hydrogen bonding network involving protein and solvent atoms.     

Prediction of protein side-chain conformation by packing optimization

We have developed a rapid and completely automatic method for prediction of protein side-chain conformation, applying the simulated annealing algorithm to optimization of side-chain packing (van der Waals) interactions. The method directly attacks the combinatorial problem of simultaneously predicting many residues' conformation, solving in 8 to 12 hours problems for which the systematic search would require over 10(300) central processing unit years. Over a test set of nine proteins ranging in size from 46 to 323 residues, the program's predictions for side-chain atoms had a root-mean-square (r.m.s.) deviation of 1.77 A overall versus the native structures. More importantly, the predictions for core residues were especially accurate, with an r.m.s. value of 1.25 A overall: 80 to 90% of the large hydrophobic side-chains dominating the internal core were correctly predicted, versus 30 to 40% for most current methods. The predictions' main errors were in surface residues poorly constrained by packing and small residues with greater steric freedom and hydrogen bonding interactions, which were not included in the program's potential function. van der Waals interactions appear to be the supreme determinant of the arrangement of side-chains in the core, enforcing a unique allowed packing that in every case so far examined matches the native structure.     

SPINFAST: using structure alignment profiles to enhance the accuracy and assess the reliability of protein side-chain modeling

We present a novel, knowledge-based method for the side-chain addition step in protein structure modeling. The foundation of the method is a conditional probability equation, which specifies the probability that a side-chain will occupy a specific rotamer state, given a set of evidence about the rotamer states adopted by the side-chains at aligned positions in structurally homologous crystal structures. We demonstrate that our method increases the accuracy of homology model side-chain addition when compared with the widely employed practice of preserving the side-chain conformation from the homology template to the target at conserved residue positions. Furthermore, we demonstrate that our method accurately estimates the probability that the correct rotamer state has been selected. This interesting result implies that our method can be used to understand the reliability of each and every side-chain in a protein homology model.     

Fast photochemical oxidation of proteins coupled with mass spectrometry

Fast photochemical oxidation of proteins (FPOP) is a hydroxyl radical footprinting approach whereby radicals, produced by UV laser photolysis of hydrogen peroxide, induce oxidation of amino acid side-chains. Mass Spectrometry (MS) is employed to locate and quantify the resulting irreversible, covalent oxidations to use as a surrogate for side-chain solvent accessibility. Modulation of oxidation levels under different conditions allows for the characterisation of protein conformation, dynamics and binding epitopes. FPOP has been applied to structurally diverse and biopharmaceutically relevant systems from small, monomeric aggregation-prone proteins to proteome-wide analysis of whole organisms. This review evaluates the current state of FPOP, the progress needed to address data analysis bottlenecks, particularly for residue-level analysis, and highlights significant developments of the FPOP platform that have enabled its versatility and complementarity to other structural biology techniques.     

Sampling the conformation of protein surface residues for flexible protein docking

Background  

The problem of determining the physical conformation of a protein dimer, given the structures of the two interacting proteins in their unbound state, is a difficult one. The location of the docking interface is determined largely by geometric complementarity, but finding complementary geometry is complicated by the flexibility of the backbone and side-chains of both proteins. We seek to generate candidates for docking that approximate the bound state well, even in cases where there is backbone and/or side-chain difference from unbound to bound states.     

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.     

New insights into intracellular lipid binding proteins: The role of buried water

The crystal structures of most intracellular lipid binding proteins (LBPs) show between 5 and 20 internally bound water molecules, depending on the presence or the absence of ligand inside the protein cavity. The structural and functional significance of these waters has been discussed for several LBPs based on studies that used various biophysical techniques. The present work focuses on two very different LBPs, heart-type fatty acid binding protein (H-FABP) and ileal lipid binding protein (ILBP). Using high-resolution nuclear magnetic resonance spectroscopy, certain resonances belonging to side-chain protons that are located inside the water-filled lipid binding cavity were observed. In the case of H-FABP, the pH- and temperature-dependent behavior of selected side-chain resonances (Ser82 OgH and the imidazole ring protons of His93) indicated an unusually slow exchange with the solvent, implying that the intricate hydrogen-bonding network of amino-acid side-chains and water molecules in the protein interior is very rigid. In addition, holo H-FABP appeared to display a reversible self-aggregation at physiological pH. For ILBP, on the other hand, a more solvent-accessible protein cavity was deduced based on the pH titration behavior of its histidine residues. Comparison with data from other LBPs implies that the evolutionary specialization of LBPs for certain ligand types was not only because of mutations of residues directly involved in ligand binding but also to a refinement of the internal water scaffold.     

Modeling side-chains using molecular dynamics improve recognition of binding region in CAPRI targets

The CAPRI-II experiment added an extra level of complexity to the problem of predicting protein-protein interactions by including 5 targets for which participants had to build or complete the 3-dimensional (3D) structure of either the receptor or ligand based on the structure of a close homolog. In this article, we describe how modeling key side-chains using molecular dynamics (MD) in explicit solvent improved the recognition of the binding region of a free energy- based computational docking method. In particular, we show that MD is able to predict with relatively high accuracy the rotamer conformation of the anchor side-chains important for molecular recognition as suggested by Rajamani et al. (Proc Natl Acad Sci USA 2004;101:11287-11292). As expected, the conformations are some of the most common rotamers for the given residue, while latch side-chains that undergo induced fit upon binding are forced into less common conformations. Using these models as starting conformations in conjunction with the rigid-body docking server ClusPro and the flexible docking algorithm SmoothDock, we produced valuable predictions for 6 of the 9 targets in CAPRI-II, missing only the 3 targets that underwent significant structural rearrangements upon binding. We also show that our free energy- based scoring function, consisting of the sum of van der Waals, Coulombic electrostatic with a distance-dependent dielectric, and desolvation free energy successfully discriminates the nativelike conformation of our submitted predictions. The latter emphasizes the critical role that thermodynamics plays on our methodology, and validates the generality of the algorithm to predict protein interactions.     

Stabilization of bound polycyclic aromatic hydrocarbons by a pi-cation interaction

Proteins can use aromatic side-chains to stabilize bound cationic ligands through cation-pi interactions. Here, we report the first example of the reciprocal process, termed pi-cation, in which a cationic protein side-chain stabilizes a neutral aromatic ligand. Site-directed mutagenesis revealed that an arginine side-chain located in the deep binding pocket of a monoclonal antibody (4D5) is essential for binding the neutral polynuclear aromatic hydrocarbon benzo[a]pyrene. This Arg was very likely selected for in the primary response, further underscoring the importance of the pi-cation interaction for ligand binding, which should be considered in protein analysis and design when ligands include aromatic groups.     

Shape and energetics of a cavity in c-Myb probed by natural and non-natural amino acid mutations.

The shape and the energetics of a functional cavity in the R2 subdomain (90-141) of the c-Myb DNA-binding domain were investigated by spectroscopy and thermodynamic analysis. We focused on the valine 103 residue located in front of the cavity. Nine mutants, in which valine 103 was substituted with alanine, 2-aminobutyric acid, norvaline, norleucine, leucine, isoleucine, allo -isoleucine, cyclohexylglycine, and cyclohexylalanine, were chemically synthesized and analyzed. These mutants provided a wide distribution of sizes which ranged from forming additional cavity space to filling and overflowing the cavity space. Temperature-scanning circular dichroism measurements and differential scanning calorimetry revealed a linear relationship between the van't Hoff enthalpy and the thermal transition temperature for the cavity-filling mutations. On the other hand, the mutants with side-chains larger than the side-chain of leucine resulted in a relatively low transition enthalpy and temperature, most likely due to the exposure of the side-chain to solvent and the increase in the entropy of the folded states. Branching at the beta-carbon atom reduced the unfolding free energy due to the steric constraint in the cavity. In particular, the mutational elongation of the side-chain from beta-carbon to the trans -to-CO direction proved to be more hindered than that from beta-carbon to the trans -to-NH. The unfolding free energy versus side-chain volume formed a bell-shaped plot with a maximum free energy for the leucine mutant. The difference in the transition free energy for cavity-filling mutants with beta-unbranched side-chains were two to four times larger than the difference in the transfer energy from organic solvent to water. Therefore, the increase in unfolding free energy would most likely be attributed to van der Waals interactions in the cavity wall, which would be a origin of stabilization by the sliding of tryptophan 95 into the cavity upon DNA binding.     

Neutron diffraction study of carbonmonoxymyoglobin.

Neutron diffraction data from a crystal of carbonmonoxymyoglobin were refined by PROLSQ, a modern restrained least-squares procedure in reciprocal space, in conjunction with a solvent analysis technique, to a final R-factor of 11.3%. The ligand CO occupies two sites and its binding conformations are distorted from the linear conformation. The N epsilon atom of the distal histidine residue is deprotonated (not deuterated), and a water molecule is bound to the N delta atom of the distal histidine. The side-chain of Lys56 (D6) exists in two alternative charge-binding sites. His24 (B5) and His119 (GH1) share a hydrogen atom. His12 (A10) and His36 (C1) are deprotonated. The deprotonated imidazole ring of His12 (A10) may act as a hydrogen-bond acceptor. The heme group is planar within 0.09 A root-mean-square (r.m.s.) deviation from planarity. The solvent environments for the two propionic acid groups are different. The side-chain of Arg45 (CD3) forms hydrogen bonds with the side-chain of Asp60 (E3) and one of the two propionic acid groups. An average N-2H . . . O angle in helical regions is 147 (+/- 11) degrees. Eleven main-chain amide hydrogen atoms from hydrophobic residues do not exchange with deuterium. The overall atomic occupancy factors for the main-chain and side-chain atoms are quite uniform, at 0.97 (+/- 0.07) and 0.93 (+/- 0.10), respectively, as shown by an occupancy analysis made at the end of the refinement procedure.