Evaluating the energetics of empty cavities and internal mutations in proteins

Clinical, genetic, and structural basis of congenital adrenal hyperplasia due to 11β-hydroxylase deficiency

Congenital adrenal hyperplasia (CAH), resulting from mutations in CYP11B1, a gene encoding 11β-hydroxylase, represents a rare autosomal recessive Mendelian disorder of aberrant sex steroid production. Unlike CAH caused by 21-hydroxylase deficiency, the disease is far more common in the Middle East and North Africa, where consanguinity is common often resulting in identical mutations. Clinically, affected female newborns are profoundly virilized (Prader score of 4/5), and both genders display significantly advanced bone ages and are oftentimes hypertensive. We find that 11-deoxycortisol, not frequently measured, is the most robust biochemical marker for diagnosing 11β-hydroxylase deficiency. Finally, computational modeling of 25 missense mutations of CYP11B1 revealed that specific modifications in the heme-binding (R374W and R448C) or substrate-binding (W116C) site of 11β-hydroxylase, or alterations in its stability (L299P and G267S), may predict severe disease. Thus, we report clinical, genetic, hormonal, and structural effects of CYP11B1 gene mutations in the largest international cohort of 108 patients with steroid 11β-hydroxylase deficiency CAH.

Peptide-Dependent Recognition of HLA-B*57:01 by KIR3DS1

Killer cell immunoglobulin-like receptors (KIRs) play an important role in the activation of natural killer (NK) cells, which in turn contribute to the effective immune control of many viral infections. In the context of HIV infection, the closely related KIR3DL1 and KIR3DS1 molecules, in particular, have been associated with disease outcome. Inhibitory signals via KIR3DL1 are disrupted by downregulation of HLA class I ligands on the infected cell surface and can also be impacted by changes in the presented peptide repertoire. In contrast, the activatory ligands for KIR3DS1 remain obscure. We used a structure-driven approach to define the characteristics of HLA class I-restricted peptides that interact with KIR3DL1 and KIR3DS1. In the case of HLA-B*57:01, we used this knowledge to identify bona fide HIV-derived peptide epitopes with similar properties. Two such peptides facilitated productive interactions between HLA-B*57:01 and KIR3DS1. These data reveal the presence of KIR3DS1 ligands within the HIV-specific peptide repertoire presented by a protective HLA class I allotype, thereby enhancing our mechanistic understanding of the processes that enable NK cells to impact disease outcome.

IMPORTANCE Natural killer (NK) cells are implicated as determinants of immune control in many viral infections, but the precise molecular mechanisms that initiate and control these responses are unclear. The activating receptor KIR3DS1 in combination with HLA-Bw4 has been associated with better outcomes in HIV infection. However, evidence of a direct interaction between these molecules is lacking. In this study, we demonstrate that KIR3DS1 recognition of HLA-Bw4 is peptide dependent. We also identify HIV-derived peptide epitopes presented by the protective HLA-B*57:01 allotype that facilitate productive interactions with KIR3DS1. Collectively, these findings suggest a mechanism whereby changes in the peptide repertoire associated with viral infection provide a trigger for KIR3DS1 engagement and NK cell activation.

Factors correlating with significant differences between X-ray structures of myoglobin

Validation of general ideas about the origins of conformational differences in proteins is critical in order to arrive at meaningful functional insights. Here, principal component analysis (PCA) and distance difference matrices are used to validate some such ideas about the conformational differences between 291 myoglobin structures from sperm whale, horse and pig. Almost all of the horse and pig structures form compact PCA clusters with only minor coordinate differences and outliers that are easily explained. The 222 whale structures form a few dense clusters with multiple outliers. A few whale outliers with a prominent distortion of the GH loop are very similar to the cluster of horse structures, which all have a similar GH-loop distortion apparently owing to intermolecular crystal lattice hydrogen bonds to the GH loop from residues near the distal histidine His64. The variations of the GH-loop coordinates in the whale structures are likely to be owing to the observed alternative intermolecular crystal lattice bond, with the change to the GH loop distorting bonds correlated with the binding of specific `unusual' ligands. Such an alternative intermolecular bond is not observed in horse myoglobins, obliterating any correlation with the ligands. Intermolecular bonds do not usually cause significant coordinate differences and cannot be validated as their universal cause. Most of the native-like whale myoglobin structure outliers can be correlated with a few specific factors. However, these factors do not always lead to coordinate differences beyond the previously determined uncertainty thresholds. The binding of unusual ligands by myoglobin, leading to crystal-induced distortions, suggests that some of the conformational differences between the apo and holo structures might not be `functionally important' but rather artifacts caused by the binding of `unusual' substrate analogs. The causes of P6 symmetry in myoglobin crystals and the relationship between crystal and solution structures are also discussed.

Contributions of water transfer energy to protein-ligand association and dissociation barriers: Watermap analysis of a series of p38α MAP kinase inhibitors

In our previous work, we proposed that desolvation and resolvation of the binding sites of proteins can serve as the slowest steps during ligand association and dissociation, respectively, and tested this hypothesis on two protein-ligand systems with known binding kinetics behavior. In the present work, we test this hypothesis on another kinetically-determined protein-ligand system-that of p38α and eight Type II BIRB 796 inhibitor analogs. The kon values among the inhibitor analogs are narrowly distributed (10⁴ ≤ kon ≤ 10⁵ M⁻¹ s⁻¹), suggesting a common rate-determining step, whereas the koff values are widely distributed (10⁻¹ ≤ koff ≤ 10⁻⁶ s⁻¹), suggesting a spectrum of rate-determining steps. We calculated the solvation properties of the DFG-out protein conformation using an explicit solvent molecular dynamics simulation and thermodynamic analysis method implemented in WaterMap to predict the enthalpic and entropic costs of water transfer to and from bulk solvent incurred upon association and dissociation of each inhibitor. The results suggest that the rate-determining step for association consists of the transfer of a common set of enthalpically favorable solvating water molecules from the binding site to bulk solvent. The rate-determining step for inhibitor dissociation consists of the transfer of water from bulk solvent to specific binding site positions that are unfavorably solvated in the apo protein, and evacuated during ligand association. Different sets of unfavorable solvation are evacuated by each ligand, and the observed dissociation barriers are qualitatively consistent with the calculated solvation free energies of those sets.

General anaesthesia: from molecular targets to neuronal pathways of sleep and arousal

The mechanisms through which general anaesthetics, an extremely diverse group of drugs, cause reversible loss of consciousness have been a long-standing mystery. Gradually, a relatively small number of important molecular targets have emerged, and how these drugs act at the molecular level is becoming clearer. Finding the link between these molecular studies and anaesthetic-induced loss of consciousness presents an enormous challenge, but comparisons with the features of natural sleep are helping us to understand how these drugs work and the neuronal pathways that they affect. Recent work suggests that the thalamus and the neuronal networks that regulate its activity are the key to understanding how anaesthetics cause loss of consciousness.

Quantifying the protein core flexibility through analysis of cavity formation

We present an extensive analysis of cavity statistics in the interior of three different proteins, in liquid n-hexane, and in water performed using molecular-dynamics simulations. The heterogeneity of packing density over atomic length scales in different parts of proteins is evident in the wide range of values observed for the average cavity size, the probability of cavity formation, and the corresponding free energy of hard-sphere insertion. More interestingly, however, the distribution of cavity sizes observed at various points in the protein interior is surprisingly homogeneous in width. That width is significantly smaller than that measured for similar distributions in liquid n-hexane or water, indicating that protein interior is much less flexible than liquid hexane. The width of the cavity size distribution correlates well with the experimental isothermal compressibility data for liquids and proteins. An analysis of cavity statistics thus provides an efficient method to quantify local properties, such as packing, stiffness, or compressibility in heterogeneous condensed media.

Structural Tolerance versus Functional Intolerance to Mutation of Hydrophobic Core Residues Surrounding Cavities in a Parvovirus Capsid

The structural and functional relevance of amino acid residues surrounding cavities within the hydrophobic core of the protein subunits that form the capsid of parvoviruses has been investigated. Several of the evolutionarily conserved, hydrophobic residues that delimit these cavities in the capsid of the minute virus of mice were replaced by other hydrophobic residues that would affect the size and/or shape of the cavity. When four or more methylene-sized groups were introduced, or six or more groups removed, capsid assembly was drastically impaired. In contrast, the introduction or removal of up to three groups had no significant effect on capsid assembly or thermostability. However, many of these mutations affected a capsid conformational transition needed for viral infectivity. Replacement of some polar residues around the largest cavity showed that capsid assembly requires a carboxylate buried within this cavity, but both aspartate and glutamate are structurally accepted. Again, only the aspartate allowed the production of infectious viruses, because of a specific role in encapsidation of the viral genome. These observations provide evidence of a remarkable structural tolerance to mutation of the hydrophobic core of the protein subunits in a viral capsid, and of an involvement of core residues and internal cavities in capsid functions needed for infectivity.

Lattice models, packing density, and Boltzmann-like distribution of cavities in proteins

A model reproducing the experimental Boltzmann-like distribution of empty cavity sizes in proteins is introduced. Proteins are represented by lattices of different dimensionalities, corresponding to different numbers of nearest neighbor contacts. Small cavities emerge and join into larger ones in a random process that can be related to random mutations. Simulations of cavity creation are performed under the constraint of a limiting total packing density. Cavities sufficiently large (20 A(3) or more), that they might accommodate at least one additional methyl group produced by a mutation, are counted and compared to the distribution of cavities according to their sizes from protein statistics. The distributions calculated with this very simple model within a realistic range of packing densities are in good agreement with the empirical cavity distribution. The results suggest that the Boltzmann-like distribution of cavities in proteins might be affected by a mechanism controlled by limiting packing density and maximum allowed protein destabilization. This supports an earlier suggestion that the agreement between the free energies of cavity formation from the mutational experiments and from the statistics of the empty cavity distribution in X-ray protein structures is nonfortuitous. A possible relation of the suggested model to the Boltzmann hypothesis is discussed.

Enzyme/non-enzyme discrimination and prediction of enzyme active site location using charge-based methods

Calculations of charge interactions complement analysis of a characterised active site, rationalising pH-dependence of activity and transition state stabilisation. Prediction of active site location through large DeltapK(a)s or electrostatic strain is relevant for structural genomics. We report a study of ionisable groups in a set of 20 enzymes, finding that false positives obscure predictive potential. In a larger set of 156 enzymes, peaks in solvent-space electrostatic properties are calculated. Both electric field and potential match well to active site location. The best correlation is found with electrostatic potential calculated from uniform charge density over enzyme volume, rather than from assignment of a standard atom-specific charge set. Studying a shell around each molecule, for 77% of enzymes the potential peak is within that 5% of the shell closest to the active site centre, and 86% within 10%. Active site identification by largest cleft, also with projection onto a shell, gives 58% of enzymes for which the centre of the largest cleft lies within 5% of the active site, and 70% within 10%. Dielectric boundary conditions emphasise clefts in the uniform charge density method, which is suited to recognition of binding pockets embedded within larger clefts. The variation of peak potential with distance from active site, and comparison between enzyme and non-enzyme sets, gives an optimal threshold distinguishing enzyme from non-enzyme. We find that 87% of the enzyme set exceeds the threshold as compared to 29% of the non-enzyme set. Enzyme/non-enzyme homologues, "structural genomics" annotated proteins and catalytic/non-catalytic RNAs are studied in this context.

Boltzmann-type distribution of side-chain conformation in proteins

We analyze packing imperfections in globular proteins as reflected in deviations of torsion angles from the equilibrium values for the isolated side chains. The distribution of conformations of methionine and lysine residues in a database of high-resolution structures is compared with energies of model compounds calculated with high-level quantum-mechanics. The distribution of the C-C and C-S torsion angles (chi(3)) correlates well with the Boltzmann factor of the torsion energy, exp(-betaE) of the model compounds C(2)H(5)-C(2)H(5) and C(2)H(5)-S-CH(3). An exponential relation was again found between the relative occurrence of g+, g- and t conformations for C(alpha)-C(beta) bonds in long side chains and the energy differences of rotamers of alpha-amino n-butyric acid, when dependence on backbone conformation was taken into account. The distribution of all 27 rotamers of methionine was correlated with the energy differences between the model's rotamers, corrected for clashes with nearby residues, the correlation being good for a set with backbone in the beta-conformation, but less clear for backbone alpha-conformation. In all correlations, the value of the coefficient beta corresponds to a temperature of circa 300 K. These results can be interpreted with a model that considers the structure of a folded protein as resulting from packing imperfectly complementary parts, with a requirement of an overall low energy. Compromises are required to optimize the fit of nonbonded contacts with surrounding groups, and side chains assume conformations away from the energy minimum. An exponential distribution is a most probable distribution, and this can be established easily under conditions other than thermal equilibrium.

Propensities, probabilities, and the Boltzmann hypothesis

The relative strengths of interactions involving polypeptide chains can be estimated with reasonable accuracy with statistical potentials, free-energy functions derived from the frequency of occurrence of structural arrangements of residues or atoms in collections of protein structures. Recent published work has shown that the energetics of side-chain/backbone interactions can be modeled by the phi/psi propensities of the 20 amino acids. In this report, the more commonly used phi/psi probabilities are demonstrated to fail in evaluating the free energies of protein conformations because of an overriding preference for all helical structures. Comparison of the hypothetical reactions implied by these two different statistics-propensities versus probabilities-leads to the conclusion that the Boltzmann hypothesis may only be applicable for the calculation of statistical potentials after the starting conformation has been specified. This conclusion supports a simple conjecture: The surprising success of the Boltzmann hypothesis in explaining the energetics of protein structures is a direct consequence of a real equilibrium, one extending over evolutionary time that has maintained the stability of each protein within a narrow range of values.

Mutagenic stabilization and/or disruption of a CD4-bound state reveals distinct conformations of the human immunodeficiency virus type 1 gp120 envelope glycoprotein

The human immunodeficiency virus type 1 (HIV-1) gp120 exterior envelope glycoprotein is conformationally flexible. Upon binding to the host cell receptor CD4, gp120 assumes a conformation that is recognized by the second receptor, CCR5 and/or CXCR4, and by the CD4-induced (CD4i) antibodies. Guided by the X-ray crystal structure of a gp120-CD4-CD4i antibody complex, we introduced changes into gp120 that were designed to stabilize or disrupt this conformation. One mutant, 375 S/W, in which the tryptophan indole group is predicted to occupy the Phe 43 cavity in the gp120 interior, apparently favors a gp120 conformation closer to that of the CD4-bound state. The 375 S/W mutant was recognized as well as or better than wild-type gp120 by CD4 and CD4i antibodies, and the large decrease in entropy observed when wild-type gp120 bound CD4 was reduced for the 375 S/W mutant. The recognition of the 375 S/W mutant by CD4BS antibodies, which are directed against the CD4-binding region of gp120, was markedly reduced compared with that of the wild-type gp120. Compared with the wild-type virus, viruses with the 375 S/W envelope glycoproteins were resistant to neutralization by IgG1b12, a CD4BS antibody, were slightly more sensitive to soluble CD4 neutralization and were neutralized more efficiently by the 2G12 antibody. Another mutant, 423 I/P, in which the gp120 bridging sheet was disrupted, did not bind CD4, CCR5, or CD4i antibodies, even though recognition by CD4BS antibodies was efficient. These results indicate that CD4BS antibodies recognize conformations of gp120 different from that recognized by CD4 and CD4i antibodies.

Composites of local structure propensities: evidence for local encoding of long-range structure

To estimate how extensively the ensemble of denatured-state conformations is constrained by local side-chain-backbone interactions, propensities of each of the 20 amino acids to occur in mono- and dipeptides mapped to discrete regions of the Ramachandran map are computed from proteins of known structure. In addition, propensities are computed for the trans, gauche-, and gauche+ rotamers, with or without consideration of the values of phi and psi. These propensities are used in scoring functions for fragment threading, which estimates the energetic favorability of fragments of protein sequence to adopt the native conformation as opposed to hundreds of thousands of incorrect conformations. As finer subdivisions of the Ramachandran plot, neighboring residue phi/psi angles, and rotamers are incorporated, scoring functions become better at ranking the native conformation as the most favorable. With the best composite propensity function, the native structure can be distinguished from 300,000 incorrect structures for 71% of the 2130 arbitrary protein segments of length 40, 48% of 2247 segments of length 30, and 20% of 2368 segments of length 20. A majority of fragments of length 30-40 are estimated to be folded into the native conformation a substantial fraction of the time. These data suggest that the variations observed in amino acid frequencies in different phi/psi/chi1 environments in folded proteins reflect energetically important local side-chain-backbone interactions, interactions that may severely restrict the ensemble of conformations populated in the denatured state to a relatively small subset with nativelike structure.

PoPMuSiC, an algorithm for predicting protein mutant stability changes: application to prion proteins

A novel tool for computer-aided design of single-site mutations in proteins and peptides is presented. It proceeds by performing in silico all possible point mutations in a given protein or protein region and estimating the stability changes with linear combinations of database-derived potentials, whose coefficients depend on the solvent accessibility of the mutated residues. Upon completion, it yields a list of the most stabilizing, destabilizing or neutral mutations. This tool is applied to mouse, hamster and human prion proteins to identify the point mutations that are the most likely to stabilize their cellular form. The selected mutations are essentially located in the second helix, which presents an intrinsic preference to form beta-structures, with the best mutations being T183-->F, T192-->A and Q186-->A. The T183 mutation is predicted to be by far the most stabilizing one, but should be considered with care as it blocks the glycosylation of N181 and this blockade is known to favor the cellular to scrapie conversion. Furthermore, following the hypothesis that the first helix might induce the formation of hydrophilic beta-aggregates, several mutations that are neutral with respect to the structure's stability but improve the helix hydrophobicity are selected, among which is E146-->L. These mutations are intended as good candidates to undergo experimental tests.

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.

Molecular mechanisms of resistance: free energy calculations of mutation effects on inhibitor binding to HIV-1 protease

The changes in the inhibitor binding constants due to the mutation of isoleucine to valine at position 84 of HIV-1 protease are calculated using molecular dynamics simulations. The calculations are done for three potent inhibitors--KNI-272, L-735,524 (indinavir or MK-639), and Ro 31-8959 (saquinavir). The calculations agree with the experimental data both in terms of an overall trend and in the magnitude of the resulting free energy change. HIV-1 protease is a homodimer, so each mutation causes two changes in the enzyme. The decrease in the binding free energy from each mutated side chain differs among the three inhibitors and correlates well with the size of the cavities induced in the protein interior near the mutated residue. The cavities are created as a result of a mutation to a smaller side chain, but the cavities are less than would be predicted from the wild-type structures, indicating that there is significant relaxation to partially fill the cavities.