Side-chain conformational changes upon protein-protein association

Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. The relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The amino acids with symmetric aromatic (Phe and Tyr) and charged (Asp and Glu) groups show the opposite trend where the near-backbone dihedral angles change the most. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr exceed the frequencies of the opposite, surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes. The results suggest that the protein association perturbs the unbound interfaces to increase the hydrophobic forces. The results facilitate better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.Comment: 21 pages, 6 figure

Experimental library screening demonstrates the successful application of computational protein design to large structural ensembles

The stability, activity, and solubility of a protein sequence are determined by a delicate balance of molecular interactions in a variety of conformational states. Even so, most computational protein design methods model sequences in the context of a single native conformation. Simulations that model the native state as an ensemble have been mostly neglected due to the lack of sufficiently powerful optimization algorithms for multistate design. Here, we have applied our multistate design algorithm to study the potential utility of various forms of input structural data for design. To facilitate a more thorough analysis, we developed new methods for the design and high-throughput stability determination of combinatorial mutation libraries based on protein design calculations. The application of these methods to the core design of a small model system produced many variants with improved thermodynamic stability and showed that multistate design methods can be readily applied to large structural ensembles. We found that exhaustive screening of our designed libraries helped to clarify several sources of simulation error that would have otherwise been difficult to ascertain. Interestingly, the lack of correlation between our simulated and experimentally measured stability values shows clearly that a design procedure need not reproduce experimental data exactly to achieve success. This surprising result suggests potentially fruitful directions for the improvement of computational protein design technology

Monte Carlo study of cooperativity in homopolypeptides

©1992 American Institute of PhysicsThe electronic version of this article is the complete one and can be found online at: http://link.aip.org/link/?JCPSA6/97/9412/1DOI:10.1063/1.463317A discretized model of globular proteins is employed in a Monte Carlo study of the helix-coil transition of polyalanine and the collapse transition of polyvaline. The present lattice realization permits real protein crystal structures to be represented at the level of 1 A resolution. Furthermore, the Monte Carlo dynamic scheme is capable of moving elements of assembled secondary and supersecondary structure. The potentials of mean force for the interactions are constructed from the statistics of a set of high resolution x-ray structures of nonhomologous proteins. The cooperativity of formation of ordered structures is found to be larger when the major contributions to the conformational energy of the low temperature states come from hydrogen bonds and short range conformational propensities. The secondary structure seen in the folded state is the result of an interplay between the short and long range interactions. Compactness itself, driven by long range, nonspecific interactions, seems to be insufficient to generate any appreciable secondary structure. A detailed examination of the dynamics of highly helical model proteins demonstrates that all elements of secondary structure are mobile in the present algorithm, and thus the folding pathways do not depend on the use of a lattice approximation. Possible applications of the present model to the prediction of protein 3D structures are briefly discussed

Subtle pH differences trigger single residue motions for moderating conformations of calmodulin

This study reveals the essence of ligand recognition mechanisms by which calmodulin (CaM) controls a variety of Ca2+ signaling processes. We study eight forms of calcium-loaded CaM each with distinct conformational states. Reducing the structure to two degrees of freedom conveniently describes main features of the conformational changes of CaM via simultaneous twist-bend motions of the two lobes. We utilize perturbation-response scanning (PRS) technique, coupled with molecular dynamics simulations. PRS is based on linear response theory, comprising sequential application of directed forces on selected residues followed by recording the resulting protein coordinates. We analyze directional preferences of the perturbations and resulting conformational changes. Manipulation of a single residue reproduces the structural change more effectively than that of single/pairs/triplets of collective modes of motion. Our findings also give information on how the flexible linker acts as a transducer of binding information to distant parts of the protein. Furthermore, by perturbing residue E31 located in one of the EF hand motifs in a specific direction, it is possible to induce conformational change relevant to five target structures. Independently, using four different pKa calculation strategies, we find this particular residue to be the charged residue (out of a total of 52), whose ionization state is most sensitive to subtle pH variations in the physiological range. It is plausible that at relatively low pH, CaM structure is less flexible. By gaining charged states at specific sites at a pH value around 7, such as E31 found in the present study, local conformational changes in the protein will lead to shifts in the energy landscape, paving the way to other conformational states. These findings are in accordance with Fluorescence Resonance Energy Transfer (FRET) measured shifts in conformational distributions towards more compact forms with decreased pH. They also corroborate mutational studies and proteolysis results which point to the significant role of E31 in CaM dynamics

Structure and substrate selectivity of the 750-kDa α6β6 holoenzyme of geranyl-CoA carboxylase.

Geranyl-CoA carboxylase (GCC) is essential for the growth of Pseudomonas organisms with geranic acid as the sole carbon source. GCC has the same domain organization and shares strong sequence conservation with the related biotin-dependent carboxylases 3-methylcrotonyl-CoA carboxylase (MCC) and propionyl-CoA carboxylase (PCC). Here we report the crystal structure of the 750-kDa α6β6 holoenzyme of GCC, which is similar to MCC but strikingly different from PCC. The structures provide evidence in support of two distinct lineages of biotin-dependent acyl-CoA carboxylases, one carboxylating the α carbon of a saturated organic acid and the other carboxylating the γ carbon of an α-β unsaturated acid. Structural differences in the active site region of GCC and MCC explain their distinct substrate preferences. Especially, a glycine residue in GCC is replaced by phenylalanine in MCC, which blocks access by the larger geranyl-CoA substrate. Mutation of this residue in the two enzymes can change their substrate preferences

Introduction to protein folding for physicists

The prediction of the three-dimensional native structure of proteins from the knowledge of their amino acid sequence, known as the protein folding problem, is one of the most important yet unsolved issues of modern science. Since the conformational behaviour of flexible molecules is nothing more than a complex physical problem, increasingly more physicists are moving into the study of protein systems, bringing with them powerful mathematical and computational tools, as well as the sharp intuition and deep images inherent to the physics discipline. This work attempts to facilitate the first steps of such a transition. In order to achieve this goal, we provide an exhaustive account of the reasons underlying the protein folding problem enormous relevance and summarize the present-day status of the methods aimed to solving it. We also provide an introduction to the particular structure of these biological heteropolymers, and we physically define the problem stating the assumptions behind this (commonly implicit) definition. Finally, we review the 'special flavor' of statistical mechanics that is typically used to study the astronomically large phase spaces of macromolecules. Throughout the whole work, much material that is found scattered in the literature has been put together here to improve comprehension and to serve as a handy reference.Comment: 53 pages, 18 figures, the figures are at a low resolution due to arXiv restrictions, for high-res figures, go to http://www.pabloechenique.co