Protein conformational substates from X-ray crystallography

Refinement of multiconformer ensemble models from multi-temperature X-ray diffraction data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363 K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Refinement of Multiconformer Ensemble Models from Multi-temperature X-ray Diffraction Data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Unlocking the secret of lignin-enzyme interactions: Recent advances in developing state-of-the-art analytical techniques

Bioconversion of renewable lignocellulosics to produce liquid fuels and chemicals is one of the most effective ways to solve the problem of fossil resource shortage, energy security, and environmental challenges. Among the many biorefinery pathways, hydrolysis of lignocellulosics to fermentable monosaccharides by cellulase is arguably the most critical step of lignocellulose bioconversion. In the process of enzymatic hydrolysis, the direct physical contact between enzymes and cellulose is an essential prerequisite for the hydrolysis to occur. However, lignin is considered one of the most recalcitrant factors hindering the accessibility of cellulose by binding to cellulase unproductively, which reduces the saccharification rate and yield of sugars. This results in high costs for the saccharification of carbohydrates. The various interactions between enzymes and lignin have been explored from different perspectives in literature, and a basic lignin inhibition mechanism has been proposed. However, the exact interaction between lignin and enzyme as well as the recently reported promotion of some types of lignin on enzymatic hydrolysis is still unclear at the molecular level. Multiple analytical techniques have been developed, and fully unlocking the secret of lignin-enzyme interactions would require a continuous improvement of the currently available analytical techniques. This review summarizes the current commonly used advanced research analytical techniques for investigating the interaction between lignin and enzyme, including quartz crystal microbalance with dissipation (QCM-D), surface plasmon resonance (SPR), attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy, atomic force microscopy (AFM), nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy (FLS), and molecular dynamics (MD) simulations. Interdisciplinary integration of these analytical methods is pursued to provide new insight into the interactions between lignin and enzymes. This review will serve as a resource for future research seeking to develop new methodologies for a better understanding of the basic mechanism of lignin-enzyme binding during the critical hydrolysis process.

Identifying and Visualizing Macromolecular Flexibility in Structural Biology

Structural biology comprises a variety of tools to obtain atomic resolution data for the investigation of macromolecules. Conventional structural methodologies including crystallography, NMR and electron microscopy often do not provide sufficient details concerning flexibility and dynamics, even though these aspects are critical for the physiological functions of the systems under investigation. However, the increasing complexity of the molecules studied by structural biology (including large macromolecular assemblies, integral membrane proteins, intrinsically disordered systems, and folding intermediates) continuously demands in-depth analyses of the roles of flexibility and conformational specificity involved in interactions with ligands and inhibitors. The intrinsic difficulties in capturing often subtle but critical molecular motions in biological systems have restrained the investigation of flexible molecules into a small niche of structural biology. Introduction of massive technological developments over the recent years, which include time-resolved studies, solution X-ray scattering, and new detectors for cryo-electron microscopy, have pushed the limits of structural investigation of flexible systems far beyond traditional approaches of NMR analysis. By integrating these modern methods with powerful biophysical and computational approaches such as generation of ensembles of molecular models and selective particle picking in electron microscopy, more feasible investigations of dynamic systems are now possible. Using some prominent examples from recent literature, we review how current structural biology methods can contribute useful data to accurately visualize flexibility in macromolecular structures and understand its important roles in regulation of biological processes.

X-ray refinement significantly underestimates the level of microscopic heterogeneity in biomolecular crystals

Biomolecular X-ray structures typically provide a static, time- and ensemble-averaged view of molecular ensembles in crystals. In the absence of rigid-body motions and lattice defects, B-factors are thought to accurately reflect the structural heterogeneity of such ensembles. In order to study the effects of averaging on B-factors, we employ molecular dynamics simulations to controllably manipulate microscopic heterogeneity of a crystal containing 216 copies of villin headpiece. Using average structure factors derived from simulation, we analyse how well this heterogeneity is captured by high-resolution molecular-replacement-based model refinement. We find that both isotropic and anisotropic refined B-factors often significantly deviate from their actual values known from simulation: even at high 1.0 Å resolution and Rfree of 5.9%, B-factors of some well-resolved atoms underestimate their actual values even sixfold. Our results suggest that conformational averaging and inadequate treatment of correlated motion considerably influence estimation of microscopic heterogeneity via B-factors, and invite caution in their interpretation.

Mapping the conformational landscape of a dynamic enzyme by multitemperature and XFEL crystallography

Determining the interconverting conformations of dynamic proteins in atomic detail is a major challenge for structural biology. Conformational heterogeneity in the active site of the dynamic enzyme cyclophilin A (CypA) has been previously linked to its catalytic function, but the extent to which the different conformations of these residues are correlated is unclear. Here we compare the conformational ensembles of CypA by multitemperature synchrotron crystallography and fixed-target X-ray free-electron laser (XFEL) crystallography. The diffraction-before-destruction nature of XFEL experiments provides a radiation-damage-free view of the functionally important alternative conformations of CypA, confirming earlier synchrotron-based results. We monitored the temperature dependences of these alternative conformations with eight synchrotron datasets spanning 100-310 K. Multiconformer models show that many alternative conformations in CypA are populated only at 240 K and above, yet others remain populated or become populated at 180 K and below. These results point to a complex evolution of conformational heterogeneity between 180-–240 K that involves both thermal deactivation and solvent-driven arrest of protein motions in the crystal. The lack of a single shared conformational response to temperature within the dynamic active-site network provides evidence for a conformation shuffling model, in which exchange between rotamer states of a large aromatic ring in the middle of the network shifts the conformational ensemble for the other residues in the network. Together, our multitemperature analyses and XFEL data motivate a new generation of temperature- and time-resolved experiments to structurally characterize the dynamic underpinnings of protein function.

DOI:http://dx.doi.org/10.7554/eLife.07574.001

Proteins are the workhorses of the cell. The shape that a protein molecule adopts enables it to carry out its role. However, a protein’s shape, or 'conformation', is not static. Instead, a protein can shift between different conformations. This is particularly true for enzymes – the proteins that catalyze chemical reactions. The region of an enzyme where the chemical reaction happens, known as the active site, often has to change its conformation to allow catalysis to proceed. Changes in temperature can also make a protein shift between alternative conformations. Understanding how a protein shifts between conformations gives insight into how it works.

A common method for studying protein conformation is X-ray crystallography. This technique uses a beam of X-rays to figure out where the atoms of the protein are inside a crystal made of millions of copies of that protein. At room temperature or biological temperature, X-rays can rapidly damage the protein. Because of this, most crystal structures are determined at very low temperatures to minimize damage. But cooling to low temperatures changes the conformations that the protein adopts, and usually causes fewer conformations to be present.

Keedy, Kenner, Warkentin, Woldeyes et al. have used X-ray crystallography from a very low temperature (-173°C or 100 K) to above room temperature (up to 27°C or 300 K) to explore the alternative conformations of an enzyme called cyclophilin A. These alternative conformations include those that have previously been linked to this enzyme’s activity. Starting at a low temperature, parts of the enzyme were seen to shift from having a single conformation to many conformations above a threshold temperature. Unexpectedly, different parts of the enzyme have different threshold temperatures, suggesting that there isn’t a single transition across the whole protein. Instead, it appears the way a protein’s conformation changes in response to temperature is more complex than was previously realized. This result suggests that conformations in different parts of a protein are coupled to each other in complex ways.

Keedy, Kenner, Warkentin, Woldeyes et al. then performed X-ray crystallography at room temperature using an X-ray free-electron laser (XFEL). This technique can capture the protein’s structure before radiation damage occurs, and confirmed that the alternative conformations observed were not affected by radiation damage.

The combination of X-ray crystallography at multiple temperatures, new analysis methods for identifying and measuring alternative conformations, and XFEL crystallography should help future studies to characterize conformational changes in other proteins.

DOI:http://dx.doi.org/10.7554/eLife.07574.002

Protein Flexibility in Drug Discovery: From Theory to Computation

Nowadays it is widely accepted that the mechanisms of biomolecular recognition are strongly coupled to the intrinsic dynamic of proteins. In past years, this evidence has prompted the development of theoretical models of recognition able to describe ligand binding assisted by protein conformational changes. On a different perspective, the need to take into account protein flexibility in structure-based drug discovery has stimulated the development of several and extremely diversified computational methods. Herein, on the basis of a parallel between the major recognition models and the simulation strategies used to account for protein flexibility in ligand binding, we sort out and describe the most innovative and promising implementations for structure-based drug discovery.

Integrative, dynamic structural biology at atomic resolution--it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing--in atomic detail--how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

Implementing electrostatic polarization cannot fill the gap between experimental and theoretical measurements for the ultrafast fluorescence decay of myoglobin

Over the past few years, time-dependent ultrafast fluorescence spectroscopy method has been applied to the study of protein dynamics. However, observations from these experiments are in a controversy with other experimental studies. Participating of theoretical methods in this debate has not reconciled the contradiction, because the predicted initial relaxation from computer simulations is one-order faster than the ultrafast fluorescence spectroscopy experiment. In those simulations, pairwise force fields are employed, which have been shown to underestimate the roughness of the free energy landscape. Therefore, the relaxation rate of protein and water molecules under pairwise force fields is falsely exaggerated. In this work, we compared the relaxations of tryptophan/environment interaction under linear response approximation employing pairwise, polarized, and polarizable force fields. Results show that although the relaxation can be slowed down to a certain extent, the large gap between experiment and theory still cannot be filled.

Vaccinia virus virulence factor N1L is a novel promising target for antiviral therapeutic intervention

The 14 kDa homodimeric N1L protein is a potent vaccinia and variola (smallpox) virulence factor. It is not essential for viral replication, but it causes a strong attenuation of viral production in culture when deleted. The N1L protein is predicted to contain the BH3-like binding domain characteristic of Bcl-2 family proteins, and it is able to bind the BH3 peptides. Its overexpression has been reported to prevent infected cells from committing apoptosis. Therefore, interfering with the N1L apoptotic blockade may be a legitimate therapeutic strategy affecting the viral growth. By using in silico ligand docking and an array of in vitro assays, we have identified submicromolar (600 nM) N1L antagonists belonging to the family of polyphenols. Their affinity is comparable to that of the BH3 peptides (70-1000 nM). We have also identified the natural polyphenol resveratrol as a moderate N1L inhibitor. Finally, we show that our ligands efficiently inhibit growth of vaccinia virus.

Modeling discrete heterogeneity in X-ray diffraction data by fitting multi-conformers

The native state of a protein is regarded to be an ensemble of conformers, which allows association with binding partners. While some of this structural heterogeneity is retained upon crystallization, reliably extracting heterogeneous features from diffraction data has remained a challenge. In this study, a new algorithm for the automatic modelling of discrete heterogeneity is presented. At high resolution, the authors' single multi-conformer model, with correlated structural features to represent heterogeneity, shows improved agreement with the diffraction data compared with a single-conformer model. The model appears to be representative of the set of structures present in the crystal. In contrast, below 2 A resolution representing ambiguous electron density by correlated multi-conformers in a single model does not yield better agreement with the experimental data. Consistent with previous studies, this suggests that variability in multi-conformer models at lower resolution levels reflects uncertainty more than coordinated motion.

Use of decoys to optimize an all-atom force field including hydration

A novel method of parameter optimization is proposed. It makes use of large sets of decoys generated for six nonhomologous proteins with different architecture. Parameter optimization is achieved by creating a free energy gap between sets of nativelike and nonnative conformations. The method is applied to optimize the parameters of a physics-based scoring function consisting of the all-atom ECEPP05 force field coupled with an implicit solvent model (a solvent-accessible surface area model). The optimized force field is able to discriminate near-native from nonnative conformations of the six training proteins when used either for local energy minimization or for short Monte Carlo simulated annealing runs after local energy minimization. The resulting force field is validated with an independent set of six nonhomologous proteins, and appears to be transferable to proteins not included in the optimization; i.e., for five out of the six test proteins, decoys with 1.7- to 4.0-A all-heavy-atom root mean-square deviations emerge as those with the lowest energy. In addition, we examined the set of misfolded structures created by Park and Levitt using a four-state reduced model. The results from these additional calculations confirm the good discriminative ability of the optimized force field obtained with our decoy sets.

Exploring structural variability in X-ray crystallographic models using protein local optimization by torsion-angle sampling

Modeling structural variability is critical for understanding protein function and for modeling reliable targets for in silico docking experiments. Because of the time-intensive nature of manual X-ray crystallographic refinement, automated refinement methods that thoroughly explore conformational space are essential for the systematic construction of structurally variable models. Using five proteins spanning resolutions of 1.0-2.8 A, it is demonstrated how torsion-angle sampling of backbone and side-chain libraries with filtering against both the chemical energy, using a modern effective potential, and the electron density, coupled with minimization of a reciprocal-space X-ray target function, can generate multiple structurally variable models which fit the X-ray data well. Torsion-angle sampling as implemented in the Protein Local Optimization Program (PLOP) has been used in this work. Models with the lowest R(free) values are obtained when electrostatic and implicit solvation terms are included in the effective potential. HIV-1 protease, calmodulin and SUMO-conjugating enzyme illustrate how variability in the ensemble of structures captures structural variability that is observed across multiple crystal structures and is linked to functional flexibility at hinge regions and binding interfaces. An ensemble-refinement procedure is proposed to differentiate between variability that is a consequence of physical conformational heterogeneity and that which reflects uncertainty in the atomic coordinates.

X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution

Crystallography supplies unparalleled detail on structural information critical for mechanistic analyses; however, it is restricted to describing low energy conformations of macromolecules within crystal lattices. Small angle X-ray scattering (SAXS) offers complementary information about macromolecular folding, unfolding, aggregation, extended conformations, flexibly linked domains, shape, conformation, and assembly state in solution, albeit at the lower resolution range of about 50 A to 10 A resolution, but without the size limitations inherent in NMR and electron microscopy studies. Together these techniques can allow multi-scale modeling to create complete and accurate images of macromolecules for modeling allosteric mechanisms, supramolecular complexes, and dynamic molecular machines acting in diverse processes ranging from eukaryotic DNA replication, recombination and repair to microbial membrane secretion and assembly systems. This review addresses both theoretical and practical concepts, concerns and considerations for using these techniques in conjunction with computational methods to productively combine solution scattering data with high-resolution structures. Detailed aspects of SAXS experimental results are considered with a focus on data interpretation tools suitable to model protein and nucleic acid macromolecular structures, including membrane protein, RNA, DNA, and protein-nucleic acid complexes. The methods discussed provide the basis to examine molecular interactions in solution and to study macromolecular flexibility and conformational changes that have become increasingly relevant for accurate understanding, simulation, and prediction of mechanisms in structural cell biology and nanotechnology.

Crystal structure analysis and solution studies of human Lck-SH3; zinc-induced homodimerization competes with the binding of proline-rich motifs

In cytosolic Src-type tyrosine kinases the Src-type homology 3 (SH3) domain binds to an internal proline-rich motif and the presence or the absence of this interaction modulates the kinase enzymatic activity. The Src-type kinase Lck plays an important role during T-cell activation and development, since it phosphorylates the T-cell antigen receptor in an early step of the activation pathway. We have determined the crystal structure of the SH3 domain from Lck kinase at a near-atomic resolution of 1.0 A. Unexpectedly, the Lck-SH3 domain forms a symmetrical homodimer in the crystal and the dimer comprises two identical zinc-binding sites in the interface. The atomic interactions formed across the dimer interface resemble strikingly those observed between SH3 domains and their canonical proline-rich ligands, since almost identical residues participate in both contacts. Ultracentrifugation experiments confirm that in the presence of zinc ions, the Lck-SH3 domain also forms dimers in solution. The Zn(2+) dissociation constant from the Lck-SH3 dimer is estimated to be lower than 100 nM. Moreover, upon addition of a proline-rich peptide with a sequence corresponding to the recognition segment of the herpesviral regulatory protein Tip, competition between zinc-induced homodimerization and binding of the peptide can be detected by both fluorescence spectroscopy and analytical ultracentrifugation. These results suggest that in vivo, too, competition between Lck-SH3 homodimerization and binding of regulatory proline-rich sequence motifs possibly represents a novel mechanism by which kinase activity is modulated. Because the residues that form the zinc-binding site are highly conserved among Lck orthologues but not in other Src-type kinases, the mechanism might be peculiar to Lck and to its role in the initial steps of T-cell activation.

Computational analysis of ligand binding dynamics at the intermolecular hot spots with the aid of simulated tempering and binding free energy calculations

Equilibrium binding dynamics is studied for a panel of benzimidazole-containing compounds at the remodeled interface between human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp), engineered by mutating to glycine hot spot residues T175 from the hormone and W104 from the receptor. Binding energetics is predicted in a good agreement with the experimental data for a panel of these small molecules that complement the engineered defect and restore the binding affinity of the wild-type hGH-hGHbp complex. The results of simulated tempering ligand dynamics at the protein-protein interface reveals a diversity of ligand binding modes that is consistent with the structural orientation of the benzimidazole ring which closely mimics the position of the mutated W104 hot spot residue in the wild-type hGH-hGHbp complex. This structural positioning of the benzimidazole core motif is shown to be a critical feature of the low-energy ligand conformations binding in the engineered cavity. The binding free energy analysis provides a plausible rationale behind the experimental dissociation constants and suggests a key role of ligand-protein van der Waals interactions in restoring binding affinity.

Heterogeneity and Inaccuracy in Protein Structures Solved by X-Ray Crystallography

Proteins are dynamic molecules, exhibiting structural heterogeneity in the form of anisotropic motion and discrete conformational substates, often of functional importance. In protein structure determination by X-ray crystallography, the observed diffraction pattern results from the scattering of X-rays by an ensemble of heterogeneous molecules, ordered and oriented by packing in a crystal lattice. The majority of proteins diffract to resolutions where heterogeneity is difficult to identify and model, and are therefore approximated by a single, average conformation with isotropic variance. Here we show that disregarding structural heterogeneity introduces degeneracy into the structure determination process, as many single, isotropic models exist that explain the diffraction data equally well. The large differences among these models imply that the accuracy of crystallographic structures has been widely overestimated. Further, it suggests that analyses that depend on small differences in the relative positions of atoms may be flawed.

Protein conformer selection by sequence-dependent packing contacts in crystals of 3-phosphoglycerate kinase

In several crystal structures of 3-phosphoglycerate kinase (PGK), the two domains occupy different relative positions. It is intriguing that the two extreme (open and closed) conformations have never been observed for the enzyme from the same species. Furthermore, in certain cases, these different crystalline conformations represent the enzyme-ligand complex of the same composition, such as the ternary complex containing either the substrate 3-phosphoglycerate (3-PG) and beta,gamma-imido-adenosine-5'-triphosphate (AMP-PNP), an analogue of the substrate MgATP, or 3-PG and the product MgADP. Thus, the protein conformation in the crystal is apparently determined by the origin of the isolated enzyme: PGK from pig muscle has only been crystallized in open conformation, whereas PGK from either Thermotoga maritima or Trypanosoma brucei has only been reported in closed conformations. A systematic analysis of the underlying sequence differences at the crucial hinge regions of the molecule and in the protein-protein contact surfaces in the crystal, in two independent pairs of open and closed states, have revealed that 1) sequential differences around the molecular hinges do not explain the appearance of fundamentally different conformations and 2) the species-specific intermolecular contacts between the nonconserved residues are responsible for stabilizing one conformation over the other in the crystalline state. A direct relationship between the steric position of the contacts in the three-dimensional structure and the conformational state of the protein has been demonstrated.

Crystal structures of engrailed homeodomain mutants: implications for stability and dynamics

We report the crystal structures and biophysical characterization of two stabilized mutants of the Drosophila Engrailed homeodomain that have been engineered to minimize electrostatic repulsion. Four independent copies of each mutant occupy the crystal lattice, and comparison of these structures illustrates variation that can be partly ascribed to networks of correlated conformational adjustments. Central to one network is leucine 26 (Leu26), which occupies alternatively two side chain rotameric conformations (-gauche and trans) and different positions within the hydrophobic core. Similar sets of conformational substates are observed in other Engrailed structures and in another homeodomain. The pattern of structural adjustments can account for NMR relaxation data and sequence co-variation networks in the wider homeodomain family. It may also explain the dysfunction associated with a P26L mutation in the human ARX homeodomain protein. Finally, we observe a novel dipolar interaction between a conserved tryptophan and a water molecule positioned along the normal to the indole ring. This interaction may explain the distinctive fluorescent properties of the homeodomain family.

FDS: flexible ligand and receptor docking with a continuum solvent model and soft-core energy function

The docking of flexible small molecule ligands to large flexible protein targets is addressed in this article using a two-stage simulation-based method. The methodology presented is a hybrid approach where the first component is a dock of the ligand to the protein binding site, based on deriving sets of simultaneously satisfied intermolecular hydrogen bonds using graph theory and a recursive distance geometry algorithm. The output structures are reduced in number by cluster analysis based on distance similarities. These structures are submitted to a modified Monte Carlo algorithm using the AMBER-AA molecular mechanics force field with the Generalized Born/Surface Area (GB/SA) continuum model. This solvent model is not only less expensive than an explicit representation, but also yields increased sampling. Sampling is also increased using a rotamer library to direct some of the protein side-chain movements along with large dihedral moves. Finally, a softening function for the nonbonded force field terms is used, enabling the potential energy function to be slowly turned on throughout the course of the simulation. The docking procedure is optimized, and the results are presented for a single complex of the arabinose binding protein. It was found that for a rigid receptor model, the X-ray binding geometry was reproduced and uniquely identified based on the associated potential energy. However, when side-chain flexibility was included, although the X-ray structure was identified, it was one of three possible binding geometries that were energetically indistinguishable. These results suggest that on relaxing the constraint on receptor flexibility, the docking energy hypersurface changes from being funnel-like to rugged. A further 14 complexes were then examined using the optimized protocol. For each complex the docking methodology was tested for a fully flexible ligand, both with and without protein side-chain flexibility. For the rigid protein docking, 13 out of the 15 test cases were able to find the experimental binding mode; this number was reduced to 11 for the flexible protein docking. However, of these 11, in the majority of cases the experimental binding mode was not uniquely identified, but was present in a cluster of low energy structures that were energetically indistinguishable. These results not only support the presence of a rugged docking energy hypersurface, but also suggest that it may be necessary to consider the possibility of more than one binding conformation during ligand optimization.

Protein side-chain rearrangement in regions of point mutations

A major problem in predicting amino acid side-chain rearrangements following point mutations is the potentially large search space. We analyzed a nonredundant data set of 393 Protein Data Bank protein pairs, each consisting of structures differing in one amino acid, to determine the number of residues changing conformation in the region of mutation. In 91-95% of cases, two or fewer residues underwent side-chain conformational change. If mutation sites with backbone displacements were excluded, the number increased to 97%. The majority of rearrangements (over 60%) were due to the inherent flexibility of side-chains, as derived from analysis of a control set of protein subunits whose crystal structures were determined more than once. Different amino acids demonstrated different degrees of flexibility near mutation sites. Large polar or charged residues, and serine, are more flexible, while the aromatic amino acids, and cysteine, are less so. This pattern is common to the inherent side-chain flexibility, as well as the increased flexibility at ligand binding sites and mutation sites. The probability for conformational change was correlated with B-factor, frequency of the side-chain conformation in proteins and solvent accessibility. The last trend was stronger for aromatic and hydrophilic residues than for hydrophobic ones. We conclude that the search space for predicting side-chain conformations in the region of mutation can be effectively restricted. However, the overall ability to predict a particular side-chain conformation, or to check predictions according to individual existing structures, is limited. These findings may be useful in deriving empirical rules for modeling side-chain conformations.

An electrostatic model for the frequency shifts in the carbonmonoxy stretching band of myoglobin: correlation of hydrogen bonding and the stark tuning rate

The effect of internal and applied external electric fields on the vibrational stretching frequency for bound CO (nu(CO)) in myoglobin mutants was studied using density functional theory. Geometry optimization and frequency calculations were carried out for an imidazole-iron-porphine-carbonmonoxy adduct with various small molecule hydrogen-bonding groups. Over 70 vibrational frequency calculations of different model geometries and hydrogen-bonding groups were compared to derive overall trends in the C-O stretching frequency (nu(CO)) in terms of the C-O bond length and Mulliken charge. Simple linear functions were derived to predict the Stark tuning rate using an approach analogous to the vibronic theory of activation.(1) Potential energy calculations show that the strongest interaction occurs for C-H or N-H hydrogen bonding nearly perpendicular to the Fe-C-O bond axis. The calculated frequencies are compared to the structural data available from 18 myoglobin crystal structures, supporting the hypothesis that the vast majority of hydrogen-bonding interactions with CO occur from the side, rather than the end, of the bound CO ligand. The nu(CO) frequency shifts agree well with experimental frequency shifts for multiple bands, known as A states, and site-directed mutations in the distal pocket of myoglobin. The model calculations quantitatively explain electrostatic effects in terms of specific hydrogen-bonding interactions with bound CO in heme proteins.

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.

Flexibility and packing in proteins

Structural flexibility is an essential attribute, without which few proteins could carry out their biological functions. Much information about protein flexibility has come from x-ray crystallography, in the form of atomic mean-square displacements (AMSDs) or B factors. Profiles showing the AMSD variation along the polypeptide chain are usually interpreted in dynamical terms but are ultimately governed by the local features of a highly complex energy landscape. Here, we bypass this complexity by showing that the AMSD profile is essentially determined by spatial variations in local packing density. On the basis of elementary statistical mechanics and generic features of atomic distributions in proteins, we predict a direct inverse proportionality between the AMSD and the contact density, i.e., the number of noncovalent neighbor atoms within a local region of approximately 1.5 nm(3) volume. Testing this local density model against a set of high-quality crystal structures of 38 nonhomologous proteins, we find that it accurately and consistently reproduces the prominent peaks in the AMSD profile and even captures minor features, such as the periodic AMSD variation within alpha helices. The predicted rigidifying effect of crystal contacts also agrees with experimental data. With regard to accuracy and computational efficiency, the model is clearly superior to its predecessors. The quantitative link between flexibility and packing density found here implies that AMSDs provide little independent information beyond that contained in the mean atomic coordinates.

Probing substates in sperm whale myoglobin using high-pressure crystallography

Pressures in the 100 MPa range are known to have an enormous number of effects on the action of proteins, but straightforward means for determining the structural basis of these effects have been lacking. Here, crystallography has been used to probe effects of pressure on sperm whale myoglobin structure. A comparison of pressure effects with those seen at low pH suggests that structural changes under pressure are interpretable as a shift in the populations of conformational substates. Furthermore, a novel high-pressure protein crystal-cooling method has been used to show low-temperature metastability, providing an alternative to room temperature, beryllium pressure cell-based techniques. The change in protein structure due to pressure is not purely compressive and involves conformational changes important to protein activity. Correlation with low-pH structures suggests observed structural changes are associated with global conformational substates. Methods developed here open up a direct avenue for exploration of the effects of pressure on proteins.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.

The 1.0 A crystal structure of Ca(2+)-bound calmodulin: an analysis of disorder and implications for functionally relevant plasticity

Calmodulin (CaM) is a highly conserved 17 kDa eukaryotic protein that can bind specifically to over 100 protein targets in response to a Ca(2+) signal. Ca(2+)-CaM requires a considerable degree of structural plasticity to accomplish this physiological role; however, the nature and extent of this plasticity remain poorly characterized. Here, we present the 1.0 A crystal structure of Paramecium tetraurelia Ca(2+)-CaM, including 36 discretely disordered residues and a fifth Ca(2+) that mediates a crystal contact. The 36 discretely disordered residues are located primarily in the central helix and the two hydrophobic binding pockets, and reveal correlated side-chain disorder that may assist target-specific deformation of the binding pockets. Evidence of domain displacements and discrete backbone disorder is provided by translation-libration-screw (TLS) analysis and multiconformer models of protein disorder, respectively. In total, the evidence for disorder at every accessible length-scale in Ca(2+)-CaM suggests that the protein occupies a large number of hierarchically arranged conformational substates in the crystalline environment and may sample a quasi-continuous spectrum of conformations in solution. Therefore, we propose that the functionally distinct forms of CaM are less structurally distinct than previously believed, and that the different activities of CaM in response to Ca(2+) may result primarily from Ca(2+)-mediated alterations in the dynamics of the protein.