Solvent Accessibility Of Residues Research Articles

Conserved salt-bridge competition triggered by phosphorylation regulates the protein interactome

Phosphorylation is a major regulator of protein interactions; however, the mechanisms by which regulation occurs are not well understood. Here we identify a salt-bridge competition or "theft" mechanism that enables a phospho-triggered swap of protein partners by Raf Kinase Inhibitory Protein (RKIP). RKIP transitions from inhibiting Raf-1 to inhibiting G-protein-coupled receptor kinase 2 upon phosphorylation, thereby bridging MAP kinase and G-Protein-Coupled Receptor signaling. NMR and crystallography indicate that a phosphoserine, but not a phosphomimetic, competes for a lysine from a preexisting salt bridge, initiating a partial unfolding event and promoting new protein interactions. Structural elements underlying the theft occurred early in evolution and are found in 10% of homo-oligomers and 30% of hetero-oligomers including Bax, Troponin C, and Early Endosome Antigen 1. In contrast to a direct recognition of phosphorylated residues by binding partners, the salt-bridge theft mechanism represents a facile strategy for promoting or disrupting protein interactions using solvent-accessible residues, and it can provide additional specificity at protein interfaces through local unfolding or conformational change.

Conformational stabilization as a strategy to prevent nucleophosmin mislocalization in leukemia

Nucleophosmin (NPM) is a nucleolar protein involved in ribosome assembly and cell homeostasis. Mutations in the C-terminal domain of NPM that impair native folding and localization are associated with acute myeloid leukemia (AML). We have performed a high-throughput screening searching for compounds that stabilize the C-terminal domain. We identified three hit compounds which show the ability to increase the thermal stability of both the C-terminal domain as well as full-length NPM. The best hit also seemed to favor folding of an AML-like mutant. Computational pocket identification and molecular docking support a stabilization mechanism based on binding of the phenyl/benzene group of the compounds to a particular hydrophobic pocket and additional polar interactions with solvent-accessible residues. Since these results indicate a chaperoning potential of our candidate hits, we tested their effect on the subcellular localization of AML-like mutants. Two compounds partially alleviated the aggregation and restored nucleolar localization of misfolded mutants. The identified hits appear promising as pharmacological chaperones aimed at therapies for AML based on conformational stabilization of NPM.

Detection of lipid-induced structural changes of the Marburg virus matrix protein VP40 using hydrogen/deuterium exchange-mass spectrometry

Marburg virus (MARV) is a lipid-enveloped virus from the Filoviridae family containing a negative sense RNA genome. One of the seven MARV genes encodes the matrix protein VP40, which forms a matrix layer beneath the plasma membrane inner leaflet to facilitate budding from the host cell. MARV VP40 (mVP40) has been shown to be a dimeric peripheral protein with a broad and flat basic surface that can associate with anionic phospholipids such as phosphatidylserine. Although a number of mVP40 cationic residues have been shown to facilitate binding to membranes containing anionic lipids, much less is known on how mVP40 assembles to form the matrix layer following membrane binding. Here we have used hydrogen/deuterium exchange (HDX) mass spectrometry to determine the solvent accessibility of mVP40 residues in the absence and presence of phosphatidylserine and phosphatidylinositol 4,5-bisphosphate. HDX analysis demonstrates that two basic loops in the mVP40 C-terminal domain make important contributions to anionic membrane binding and also reveals a potential oligomerization interface in the C-terminal domain as well as a conserved oligomerization interface in the mVP40 N-terminal domain. Lipid binding assays confirm the role of the two basic patches elucidated with HD/X measurements, whereas molecular dynamics simulations and membrane insertion measurements complement these studies to demonstrate that mVP40 does not appreciably insert into the hydrocarbon region of anionic membranes in contrast to the matrix protein from Ebola virus. Taken together, we propose a model by which association of the mVP40 dimer with the anionic plasma membrane facilitates assembly of mVP40 oligomers.

Dimethyl Labeling Coupled with Mass Spectrometry for Topographical Characterization of Primary Amines on Monoclonal Antibodies.

Site-specific solvent accessibility of the primary amines (mainly lysine or the N-termini) on proteins is of great interest in many research areas because amines are an important functional group for protein conjugation. In this study, we coupled dimethyl labeling via reductive amination with liquid chromatography-mass spectrometry (LC-MS) to fully characterize the solvent accessibility of lysine residues and the N-termini on human immunoglobulin G (IgG). Circular dichroism (CD) and fluorescence spectroscopy revealed that dimethyl labeling did not alter the conformation of the native IgG molecule. Based on intact protein measurements, up to 28 (light chain) and 66 (heavy chain) dimethyl tags, covering all lysine residues and the N-termini, were sequentially incorporated into IgG molecules in 1000 s. All labeled sites were identified and quantified by a bottom-up proteomics approach. Some highly exposed hot-spots (for example, the N-termini of both the heavy and the light chains) and some buried sites (for example, K415 in the heavy chain and K39 in the light chain) were unambiguously revealed. This method was also used to characterize aggregation-induced structural changes in IgGs by increasing the temperature. Substantial changes in the labeling percentage of many lysine sites were observed, indicating a non-native aggregation triggered by thermal stress. Due to high labeling yields and the van der Waals surface of the labeling reagents being comparable to that of water, dimethyl labeling is a highly promising technique for probing the amine's surface topography of proteins. It can also be used as a complementary approach to other methods for resolving the higher-order structure of proteins by LC-MS.

Improving prediction of burial state of residues by exploiting correlation among residues

BackgroundResidues in a protein might be buried inside or exposed to the solvent surrounding the protein. The buried residues usually form hydrophobic cores to maintain the structural integrity of proteins while the exposed residues are tightly related to protein functions. Thus, the accurate prediction of solvent accessibility of residues will greatly facilitate our understanding of both structure and functionalities of proteins. Most of the state-of-the-art prediction approaches consider the burial state of each residue independently, thus neglecting the correlations among residues.ResultsIn this study, we present a high-order conditional random field model that considers burial states of all residues in a protein simultaneously. Our approach exploits not only the correlation among adjacent residues but also the correlation among long-range residues. Experimental results showed that by exploiting the correlation among residues, our approach outperformed the state-of-the-art approaches in prediction accuracy. In-depth case studies also showed that by using the high-order statistical model, the errors committed by the bidirectional recurrent neural network and chain conditional random field models were successfully corrected.ConclusionsOur methods enable the accurate prediction of residue burial states, which should greatly facilitate protein structure prediction and evaluation.

Denatured State Loop Formation Thermodynamics of a Hybrid Polypeptide

Previous work with the four-helix bundle protein cytochrome c’ in Rhodopseudomonas palustris using histidine-heme loop formation thermodynamic methods revealed fold-specific deviations from random coil character in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix bundle polypeptide, the human DNA excision repair protein's second ubiquitin-associated (UBA) domain, UBA(2). We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain to the N-terminus of iso-1-cytochrome c. Using site-directed mutagenesis, we have engineered histidine into solvent accessible surface residue positions within the all-alpha fold, creating eight single histidine variants. Isothermal equilibration denaturation studies reveal that the fusion protein unfolds in a 3-state process, commencing with iso-1-cytochrome c followed by UBA(2). Thermodynamic stability experiments also demonstrate that the histidine residues in the UBA(2) domain strongly destabilize iso-1-cytochrome c. Furthermore, histidine-heme loop formation equilibria show lower apparent pKa's compared to the pseudo-wild type variant, indicating significant interactions in the denatured state. We will compare the degree of deviation of loop stability versus loop size, relative to predictions of the Jacobson-Stockmayer relationship used in our previous work on cytochrome c’. This comparison will allow evaluation of sequence-based conformational bias in the denatured state of this protein.

Surface Accessibility and Dynamics of Macromolecular Assemblies Probed by Covalent Labeling Mass Spectrometry and Integrative Modeling.

Mass spectrometry (MS) has become an indispensable tool for investigating the architectures and dynamics of macromolecular assemblies. Here we show that covalent labeling of solvent accessible residues followed by their MS-based identification yields modeling restraints that allow mapping the location and orientation of subunits within protein assemblies. Together with complementary restraints derived from cross-linking and native MS, we built native-like models of four heterocomplexes with known subunit structures and compared them with available X-ray crystal structures. The results demonstrated that covalent labeling followed by MS markedly increased the predictive power of the integrative modeling strategy enabling more accurate protein assembly models. We applied this strategy to the F-type ATP synthase from spinach chloroplasts (cATPase) providing a structural basis for its function as a nanomotor. By subjecting the models generated by our restraint-based strategy to molecular dynamics (MD) simulations, we revealed the conformational states of the peripheral stalk and assigned flexible regions in the enzyme. Our strategy can readily incorporate complementary chemical labeling strategies and we anticipate that it will be applicable to many other systems providing new insights into the structure and function of protein complexes.

A hydrophobic spine stabilizes a surface-exposed α-helix according to analysis of the solvent-accessible surface area.

BackgroundMost of hydrophilic and hydrophobic residues are thought to be exposed and buried in proteins, respectively. In contrast to the majority of the existing studies on protein folding characteristics using protein structures, in this study, our aim was to design predictors for estimating relative solvent accessibility (RSA) of amino acid residues to discover protein folding characteristics from sequences.MethodsThe proposed 20 real-value RSA predictors were designed on the basis of the support vector regression method with a set of informative physicochemical properties (PCPs) obtained by means of an optimal feature selection algorithm. Then, molecular dynamics simulations were performed for validating the knowledge discovered by analysis of the selected PCPs.ResultsThe RSA predictors had the mean absolute error of 14.11% and a correlation coefficient of 0.69, better than the existing predictors. The hydrophilic-residue predictors preferred PCPs of buried amino acid residues to PCPs of exposed ones as prediction features. A hydrophobic spine composed of exposed hydrophobic residues of an α-helix was discovered by analyzing the PCPs of RSA predictors corresponding to hydrophobic residues. For example, the results of a molecular dynamics simulation of wild-type sequences and their mutants showed that proteins 1MOF and 2WRP_H16I (Protein Data Bank IDs), which have a perfectly hydrophobic spine, have more stable structures than 1MOF_I54D and 2WRP do (which do not have a perfectly hydrophobic spine).ConclusionsWe identified informative PCPs to design high-performance RSA predictors and to analyze these PCPs for identification of novel protein folding characteristics. A hydrophobic spine in a protein can help to stabilize exposed α-helices.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1368-z) contains supplementary material, which is available to authorized users.

A comparative analysis of the spatial structure of nonspecific porins from Yersinia ruckeri using optical spectroscopy and molecular modeling

The spatial organization of outer-membrane porins is studied by optical spectroscopy and molecular modeling. It was found that the OmpF and OmpC porins from Yеrsiniа ruckeri are β-structured membrane proteins typical of the pore-forming proteins of other Gram-negative bacteria. The spatial structures of monomers and trimers of the OmpC and OmpF porins from Y. ruckeri are simulated using methods of structural bioinformatics. It was found that the structural stability of the more thermostable OmpF trimer is sustained by a greater number of hydrogen bonds and hydrophobic interactions. The main differences of the spatial structures of the test porins are observed in the structure of their outer loops. There are three tryptophan residues in the molecules of the OmpC and OmpF porins of Y. ruckeri. It is demonstrated by moleculardynamics methods that after thermal denaturation the solvent accessibility of the Trp212 residue in OmpF porin increased by two times, while the solvent accessibility of a Trp184 residue in OmpC porin was not increased. It is hypothesized that the red-shifted tryptophan fluorescence spectrum of OmpF porin during thermal denaturation is due to the behavior of the Trp212 residue.

OR1 Mapping anti-HLA antibody epitopes by painting the HLA blue

Aim Mapping anti-HLA antibody epitopes has commonly been done by comparing antibody reactivity patterns to a panel of naturally occurring HLA allelic variants such as the solid phase HLA antibody screening panels. Information including the donor’s immunizing alleles and the recipient-derived HLA further helps to refine HLA antibody epitopes. While these methods are inexpensive and widely available, they represent an indirect means of mapping those amino acids on HLA that are critical to antibody recognition and binding. Here we test a complementary direct approach to measure the ‘foot print’ of an antibody at it interacts with HLA. A recently developed method called protein painting was used to map the interaction of mouse monoclonal antibody BB7.2 with HLA-A*02:01. Methods Briefly, purified soluble HLA-A*02:01 was incubated with and without BB7.2. Protein mixtures were incubated with the trypsin inhibitor Remazol Brilliant Blue (RBB), trypsin digested, and tryptic peptide fragments of HLA-A*02:01 were analyzed using LCMS. In this manner BBB protects solvent accessible Arg and Lys residues from trypsin digestion. When HLA-A*02:01 is first incubated with BB7.2, Arg and Lys residues covered by the antibody are not painted by BBB and therefore they remain trypsin susceptible. Results A mass spectrometric analysis of the resulting trypic peptides revealed that BB7.2 protects R44, R97, R108, R111, R131, and K144 from painting. These residues are primarily located on lateral edge of the alpha-2 alpha helix as shown on the three dimensional structure of HLA-A*02:01. In a solid phase assay, BB7.2 reacts exclusively with A2 and A28 allomorphs, data that identifies the well defined 144TKH epitope. The protein painting data shown here are consistent with the 144TKH epitope as K144 was protected from RBB binding by BB7.2. Conclusion We demonstrate a new method for direct mapping of anti-HLA antibody epitopes that is complementary to solid phase assays. Furthermore, protein painting reveals not only the critical residues for antibody binding, but the surrounding amino acids covered by the antibody. R. Buchli : Employee; Company/Organization; Pure Protein LLC. W. Hildebrand: Scientific/Medical Advisor; Company/Organization; Pure Protein LLC.

A role for loop G in the β1 strand in GABAA receptor activation.

The role of the β1 strand in GABAA receptor function is unclear. It lies anti-parallel to the β2 strand, which is known to participate in receptor activation. Molecular dynamics simulation revealed solvent accessible residues within the β1 strand of the GABAA β3 homopentamer that might be amenable to analysis using the substituted Cys accessibility method. Cys substitutions from Asp43 to Thr47 in the GABAA α1 subunit showed that D43C and T47C reduced the apparent potency of GABA. F45C caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Cys43 and Cys47 were accessible to 2-aminoethyl methanethiosulphonate (MTSEA) modification, whereas Cys45 was not. Both GABA and the allosteric agonist propofol reduced MTSEA modification of Cys43 and Cys47. By contrast, modification of Cys64 in the β2 strand loop D was impeded by GABA but unaffected by propofol. These data reveal movement of β1 strand loop G residues during agonist activation of the GABAA receptor. The GABAA receptor α subunit β1 strand runs anti-parallel to the β2 strand, which contains loop D, known to participate in receptor activation and agonist binding. However, a role for the β1 strand has yet to be established. We used molecular dynamics simulation to quantify the solvent accessible surface area (SASA) of β1 strand residues in the GABAA β3 homopentamer structure. Residues in the complementary interface equivalent to those between Asp43 and Thr47 in the α1 subunit have an alternating pattern of high and low SASA consistent with a β strand structure. We investigated the functional role of these β1 strand residues in the α1 subunit by individually replacing them with Cys residues. D43C and T47C substitutions reduced the apparent potency of GABA at α1β2γ2 receptors by 50-fold and eight-fold, respectively, whereas the F45C substitution caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Receptors with D43C or T47C substitutions were sensitive to 2-aminoethyl methanethiosulphonate (MTSEA) modification. However, GABA-evoked currents mediated by α1(F45C)β2γ2 receptors were unaffected by MTSEA, suggesting that this residue is inaccessible. Both GABA and the allosteric agonist propofol reduced MTSEA modification of α1(D43C)β2γ2 and α1(T47C)β2γ2 receptors, indicating movement of the β1 strand even during allosteric activation. This is in contrast to α1(F64C)β2γ2 receptors, where only GABA, but not propofol, reduced MTSEA modification. These findings provide the first functional evidence for movement of the β1 strand during gating of the receptor and identify residues that are critical for maintaining GABAA receptor function.

FPOP-LC-MS/MS Suggests Differences in Interaction Sites of Amphipols and Detergents with Outer Membrane Proteins

Amphipols are a class of novel surfactants that are capable of stabilizing the native state of membrane proteins. They have been shown to be highly effective, in some cases more so than detergent micelles, at maintaining the structural integrity of membrane proteins in solution, and have shown promise as vehicles for delivering native membrane proteins into the gas phase for structural interrogation. Here, we use fast photochemical oxidation of proteins (FPOP), which irreversibly labels the side chains of solvent-accessible residues with hydroxyl radicals generated by laser photolysis of hydrogen peroxide, to compare the solvent accessibility of the outer membrane protein OmpT when solubilized with the amphipol A8-35 or with n-dodecyl-β-maltoside (DDM) detergent micelles. Using quantitative mass spectrometry analyses, we show that fast photochemical oxidation reveals differences in the extent of solvent accessibility of residues between the A8-35 and DDM solubilized states, providing a rationale for the increased stability of membrane proteins solubilized with amphipol compared with detergent micelles, as a result of additional intermolecular contacts.Graphical ᅟ

Functional Sites Induce Long-Range Evolutionary Constraints in Enzymes.

Functional residues in proteins tend to be highly conserved over evolutionary time. However, to what extent functional sites impose evolutionary constraints on nearby or even more distant residues is not known. Here, we report pervasive conservation gradients toward catalytic residues in a dataset of 524 distinct enzymes: evolutionary conservation decreases approximately linearly with increasing distance to the nearest catalytic residue in the protein structure. This trend encompasses, on average, 80% of the residues in any enzyme, and it is independent of known structural constraints on protein evolution such as residue packing or solvent accessibility. Further, the trend exists in both monomeric and multimeric enzymes and irrespective of enzyme size and/or location of the active site in the enzyme structure. By contrast, sites in protein–protein interfaces, unlike catalytic residues, are only weakly conserved and induce only minor rate gradients. In aggregate, these observations show that functional sites, and in particular catalytic residues, induce long-range evolutionary constraints in enzymes.

Smyd3 open & closed lock mechanism for substrate recruitment: The hinge motion of C-terminal domain inferred from μ-second molecular dynamics simulations

BackgroundThe human lysine methyltransferase Smyd3, a member of the SET and MYND domain containing protein family, harbors methylation activity on both histone and non-histone targets in a tightly regulated manner. The mechanism of how Smyd3 dynamically regulates substrate recognition is still not fully unveiled. MethodsHere, we employed molecular dynamics simulations on full length human Smyd3, performed to a total of 1.2μ-second, in the presence (holo) and absence (apo) of the S-Adenosyl methionine (AdoMet) cofactor. The dynamical features of Smyd3 in apo and holo states have been examined and compared via examining geometrical and electrostatic properties. ResultsThe results show a distinct dynamics of the C-terminal domain (CTD) in the two states. In the apo state, the CTD undergoes a large hinge like motion and samples more opened configurations, thus acting like a loosened clamp and resulting in expanded substrate binding crevice. In the holo state, the CTD exhibits a restricted motion while the overall structure remains compact, mimicking a closed clamp. This leads to a localized increase in the negative potential at the substrate binding cleft. Further, solvent accessibility of critical residues at the target lysine access channel, important for methylation activity, is increased. ConclusionsWe postulate that AdoMet cofactor acts like a key and locks Smyd3 in a closed conformation. In effect, the cofactor binding restricts the elasticity of the CTD, presenting a compact substrate binding cleft with high negative potential, which may have implications on substrate recruitment via long range electrostatics. General significanceThe deletion of the CTD from Smyd3 has been shown to abolish the basal histone methylation activity. Our study highlights the importance of the CTD elasticity in shaping the substrate binding site for recognition and supports the previously proposed role of the CTD in stabilizing the active site for methylation activity.

Dissecting the roles of local packing density and longer-range effects in protein sequence evolution.

What are the structural determinants of protein sequence evolution? A number of site-specific structural characteristics have been proposed, most of which are broadly related to either the density of contacts or the solvent accessibility of individual residues. Most importantly, there has been disagreement in the literature over the relative importance of solvent accessibility and local packing density for explaining site-specific sequence variability in proteins. We show that this discussion has been confounded by the definition of local packing density. The most commonly used measures of local packing, such as contact number and the weighted contact number, represent the combined effects of local packing density and longer-range effects. As an alternative, we propose a truly local measure of packing density around a single residue, based on the Voronoi cell volume. We show that the Voronoi cell volume, when calculated relative to the geometric center of amino-acid side chains, behaves nearly identically to the relative solvent accessibility, and each individually can explain, on average, approximately 34% of the site-specific variation in evolutionary rate in a data set of 209 enzymes. An additional 10% of variation can be explained by nonlocal effects that are captured in the weighted contact number. Consequently, evolutionary variation at a site is determined by the combined effects of the immediate amino-acid neighbors of that site and effects mediated by more distant amino acids. We conclude that instead of contrasting solvent accessibility and local packing density, future research should emphasize on the relative importance of immediate contacts and longer-range effects on evolutionary variation. Proteins 2016; 84:841-854. © 2016 Wiley Periodicals, Inc.

An in silico approach to investigate the source of the controversial interpretations about the phenotypic results of the human AhR-gene G1661A polymorphism

Aryl hydrocarbon receptor (AhR) acts as an enhancer binding ligand-activated intracellular receptor. Chromatin remodeling components and general transcription factors such as TATA-binding protein (TBP) are evoked on AhR-target genes by interaction with its flexible transactivation domain (TAD). AhR-G1661A single nucleotide polymorphism (SNP: rs2066853) causes an arginine to lysine substitution in the acidic sub-domain of TAD at position 554 (R554K). Although, numerous studies associate the SNP with some abnormalities such as cancer, other reliable investigations refuse the associations. Consequently, the interpretation of the phenotypic results of G1661A-transition has been controversial. In this study, an in silico analysis were performed to investigate the possible effects of the transition on AhR-mRNA, protein structure, interaction properties and modifications. The analysis revealed that the R554K substitution affects secondary structure and solvent accessibility of adjacent residues. Also, it causes to decreasing of the AhR stability; altering the hydropathy features of the local sequence and changing the pattern of the residues at the binding site of the TAD-acidic sub-domain. Generating of new sites for ubiquitination and acetylation for AhR-K554 variant respectively at positions 544 and 560 was predicted. Our findings intensify the idea that the AhR-G1661A transition may affects AhR-TAD interactions, especially with the TBP, which influence AhR-target genes expression. However, the previously reported flexibility of the modular TAD could act as an intervening factor, moderate the SNP effects and causes distinct outcomes in different individuals and tissues.

Highly accurate sequence-based prediction of half-sphere exposures of amino acid residues in proteins.

Solvent exposure of amino acid residues of proteins plays an important role in understanding and predicting protein structure, function and interactions. Solvent exposure can be characterized by several measures including solvent accessible surface area (ASA), residue depth (RD) and contact numbers (CN). More recently, an orientation-dependent contact number called half-sphere exposure (HSE) was introduced by separating the contacts within upper and down half spheres defined according to the Cα-Cβ (HSEβ) vector or neighboring Cα-Cα vectors (HSEα). HSEα calculated from protein structures was found to better describe the solvent exposure over ASA, CN and RD in many applications. Thus, a sequence-based prediction is desirable, as most proteins do not have experimentally determined structures. To our best knowledge, there is no method to predict HSEα and only one method to predict HSEβ. This study developed a novel method for predicting both HSEα and HSEβ (SPIDER-HSE) that achieved a consistent performance for 10-fold cross validation and two independent tests. The correlation coefficients between predicted and measured HSEβ (0.73 for upper sphere, 0.69 for down sphere and 0.76 for contact numbers) for the independent test set of 1199 proteins are significantly higher than existing methods. Moreover, predicted HSEα has a higher correlation coefficient (0.46) to the stability change by residue mutants than predicted HSEβ (0.37) and ASA (0.43). The results, together with its easy Cα-atom-based calculation, highlight the potential usefulness of predicted HSEα for protein structure prediction and refinement as well as function prediction. The method is available at http://sparks-lab.org yuedong.yang@griffith.edu.au or yaoqi.zhou@griffith.edu.au Supplementary data are available at Bioinformatics online.

Logistic regression models to predict solvent accessible residues using sequence- and homology-based qualitative and quantitative descriptors applied to a domain-complete X-ray structure learning set.

A working example of relative solvent accessibility (RSA) prediction for proteins is presented. Novel logistic regression models with various qualitative descriptors that include amino acid type and quantitative descriptors that include 20- and six-term sequence entropy have been built and validated. A domain-complete learning set of over 1300 proteins is used to fit initial models with various sequence homology descriptors as well as query residue qualitative descriptors. Homology descriptors are derived from BLASTp sequence alignments, whereas the RSA values are determined directly from the crystal structure. The logistic regression models are fitted using dichotomous responses indicating buried or accessible solvent, with binary classifications obtained from the RSA values. The fitted models determine binary predictions of residue solvent accessibility with accuracies comparable to other less computationally intensive methods using the standard RSA threshold criteria 20 and 25% as solvent accessible. When an additional non-homology descriptor describing Lobanov-Galzitskaya residue disorder propensity is included, incremental improvements in accuracy are achieved with 25% threshold accuracies of 76.12 and 74.79% for the Manesh-215 and CASP(8+9) test sets, respectively. Moreover, the described software and the accompanying learning and validation sets allow students and researchers to explore the utility of RSA prediction with simple, physically intuitive models in any number of related applications.

Structure of the ordered hydration of amino acids in proteins: analysis of crystal structures.

Crystallography provides unique information about the arrangement of water molecules near protein surfaces. Using a nonredundant set of 2818 protein crystal structures with a resolution of better than 1.8 Å, the extent and structure of the hydration shell of all 20 standard amino-acid residues were analyzed as function of the residue conformation, secondary structure and solvent accessibility. The results show how hydration depends on the amino-acid conformation and the environment in which it occurs. After conformational clustering of individual residues, the density distribution of water molecules was compiled and the preferred hydration sites were determined as maxima in the pseudo-electron-density representation of water distributions. Many hydration sites interact with both main-chain and side-chain amino-acid atoms, and several occurrences of hydration sites with less canonical contacts, such as carbon-donor hydrogen bonds, OH-π interactions and off-plane interactions with aromatic heteroatoms, are also reported. Information about the location and relative importance of the empirically determined preferred hydration sites in proteins has applications in improving the current methods of hydration-site prediction in molecular replacement, ab initio protein structure prediction and the set-up of molecular-dynamics simulations.

The utility of protein structure as a predictor of site-wise dN/dS varies widely among HIV-1 proteins.

Protein structure acts as a general constraint on the evolution of viral proteins. One widely recognized structural constraint explaining evolutionary variation among sites is the relative solvent accessibility (RSA) of residues in the folded protein. In influenza virus, the distance from functional sites has been found to explain an additional portion of the evolutionary variation in the external antigenic proteins. However, to what extent RSA and distance from a reference site in the protein can be used more generally to explain protein adaptation in other viruses and in the different proteins of any given virus remains an open question. To address this question, we have carried out an analysis of the distribution and structural predictors of site-wise dN/dS in HIV-1. Our results indicate that the distribution of dN/dS in HIV follows a smooth gamma distribution, with no special enrichment or depletion of sites with dN/dS at or above one. The variation in dN/dS can be partially explained by RSA and distance from a reference site in the protein, but these structural constraints do not act uniformly among the different HIV-1 proteins. Structural constraints are highly predictive in just one of the three enzymes and one of three structural proteins in HIV-1. For these two proteins, the protease enzyme and the gp120 structural protein, structure explains between 30 and 40% of the variation in dN/dS. Finally, for the gp120 protein of the receptor-binding complex, we also find that glycosylation sites explain just 2% of the variation in dN/dS and do not explain gp120 evolution independently of either RSA or distance from the apical surface.