SEARCH · PubMed Health

Results for “Solvent accessibility”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33Linked to original sources

PubMed1997

The critical role of a solvent-exposed residue of an MHC class I-restricted peptide in MHC-peptide binding.

The immunodominant ovalbumin257-264 (OVA-8, SIINFEKL) and herpes simplex virus gB496-503 (HSV-8, SSIEFARL) peptides share 50% amino acid identity (residues P1, P3, P5 and P8) and bind with comparable efficacy to the murine MHC-encoded class I molecule H-2Kb. However, these two peptides bind differently to H-2Kbm8, a natural H-2Kb variant with a substitution in four amino acids on the floor of the peptide-binding site; HSV-8 binds with high and OVA-8 with a relatively low efficacy. To investigate which of the non-homologous peptide residues were responsible for this differential binding, we used substituted peptide variants and the class I thermodynamic stabilization assay. Variation at the solvent-exposed peptide residues P6 and P7 did not appreciably influence binding. By contrast, variation at the buried P2 and, surprisingly, at the solvent-exposed P4 residue was found to be important. Transplantation of the HSV-8 P2 or P4 residues onto the OVA-8 backbone created variant peptides O2S (P2I-->S) and O4E (P4N-->E) that bound considerably better to H-2Kbm8 than OVA-8. Furthermore, the double-substituted peptide, O2S4E, bound even better, revealing a cooperative effect of the two residues. The reciprocally substituted peptides H2I and H4N, generated by grafting the OVA-8 P2 and P4 residues onto the HSV-8 backbone respectively, bound to H-2Kbm8 slightly worse than HSV-8 but the double-substituted peptide H2I4N bound as poorly as OVA-8. Effects exerted by the P4 residue, which is solvent accessible and therefore available for the TCR contact, demonstrated that exposed peptide residues can, in certain situations, influence not only the TCR contact but also MHC-peptide binding.

Animals↗

PubMed2005

Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: calculation of structural and dynamic properties.

Long dynamics simulations were carried out on the B1 immunoglobulin-binding domain of streptococcal protein G (ProtG) and bovine pancreatic trypsin inhibitor (BPTI) using atomistic descriptions of the proteins and a continuum representation of solvent effects. To mimic frictional and random collision effects, Langevin dynamics (LD) were used. The main goal of the calculations was to explore the stability of tens-of-nanosecond trajectories as generated by this molecular mechanics approximation and to analyze in detail structural and dynamical properties. Conformational fluctuations, order parameters, cross correlation matrices, residue solvent accessibilities, pKa values of titratable groups, and hydrogen-bonding (HB) patterns were calculated from all of the trajectories and compared with available experimental data. The simulations comprised over 40 ns per trajectory for ProtG and over 30 ns per trajectory for BPTI. For comparison, explicit water molecular dynamics simulations (EW/MD) of 3 ns and 4 ns, respectively, were also carried out. Two continuum simulations were performed on each protein using the CHARMM program, one with the all-atom PAR22 representation of the protein force field (here referred to as PAR22/LD simulations) and the other with the modifications introduced by the recently developed CMAP potential (CMAP/LD simulations). The explicit solvent simulations were performed with PAR22 only. Solvent effects are described by a continuum model based on screened Coulomb potentials (SCP) reported earlier, i.e., the SCP-based implicit solvent model (SCP-ISM). For ProtG, both the PAR22/LD and the CMAP/LD 40-ns trajectories were stable, yielding C(alpha) root mean square deviations (RMSD) of about 1.0 and 0.8 A respectively along the entire simulation time, compared to 0.8 A for the EW/MD simulation. For BPTI, only the CMAP/LD trajectory was stable for the entire 30-ns simulation, with a C(alpha) RMSD of approximately 1.4 A, while the PAR22/LD trajectory became unstable early in the simulation, reaching a C(alpha) RMSD of about 2.7 A and remaining at this value until the end of the simulation; the C(alpha) RMSD of the EW/MD simulation was about 1.5 A. The source of the instabilities of the BPTI trajectories in the PAR22/LD simulations was explored by an analysis of the backbone torsion angles. To further validate the findings from this analysis of BPTI, a 35-ns SCP-ISM simulation of Ubiquitin (Ubq) was carried out. For this protein, the CMAP/LD simulation was stable for the entire simulation time (C(alpha) RMSD of approximately 1.0 A), while the PAR22/LD trajectory showed a trend similar to that in BPTI, reaching a C(alpha) RMSD of approximately 1.5 A at 7 ns. All the calculated properties were found to be in agreement with the corresponding experimental values, although local deviations were also observed. HB patterns were also well reproduced by all the continuum solvent simulations with the exception of solvent-exposed side chain-side chain (sc-sc) HB in ProtG, where several of the HB interactions observed in the crystal structure and in the EW/MD simulation were lost. The overall analysis reported in this work suggests that the combination of an atomistic representation of a protein with a CMAP/CHARMM force field and a continuum representation of solvent effects such as the SCP-ISM provides a good description of structural and dynamic properties obtained from long computer simulations. Although the SCP-ISM simulations (CMAP/LD) reported here were shown to be stable and the properties well reproduced, further refinement is needed to attain a level of accuracy suitable for more challenging biological applications, particularly the study of protein-protein interactions.

Alanine↗

PubMed1997

Temperature dependence of histidine ionization constants in myoglobin.

The standard enthalpy of ionization of six titratable histidines in horse metaquomyoglobin was determined by repeating proton NMR titrations as a function of temperature and using the van't Hoff relationship. It was found that deltaH degrees varies between 16 and 37 kJ mol(-1) in the protein, compared with a value of 29 kJ mol(-1) in free histidine. The standard entropy change was evaluated by combining the enthalpy and free energy changes derived from the pKa values. Although the entropy change could not be precisely and accurately obtained by this method, it could be established that it spans a wide range, from -60 to 0 J K(-1) mol(-1), about the value of -23 J K(-1) mol(-1) for the free histidine. The entropy change was used within the framework of enthalpy-entropy compensation to partition the solvation component from the standard thermodynamic quantities for each of the titrating residues. It was shown that the partitioning of the values in the protein is not readily understood in terms of solvent accessibility or electrostatic interactions. The contribution of solvation effects to the temperature response appeared to be significant only in the case of His-119 and His-48. The standard quantities were also used to explore the energetics of proton binding in the native state at temperatures below the onset of thermal denaturation.

Animals↗

PubMedDate not supplied

Searching for patterns of amino acids in 3D protein structures.

This paper describes the program ASSAM, which has been developed to search for patterns of amino acid side-chains in the 3D structures in the Protein Data Bank. ASSAM represents an amino acid by a vector drawn from the main chain towards the functional part of the amino acid and then computes a graph representation of a protein in which the individual side-chain vectors are the nodes and the intervector distances are the edges. The presence of a query pattern in a Protein Data Bank structure can then be searched for by means of a subgraph isomorphism algorithm. Recent enhancements to ASSAM allow searches to include the following: the main-chain structure in addition to the side-chains; the secondary structure and solvent accessibility of side-chains; allowable distances from a known binding-site; disulfide bridges; and improved generic and wild-card queries. The effectiveness of these approaches is demonstrated by extensive searches of the Protein Data Bank for typical 3D query patterns.

Algorithms↗

PubMed2005

SCRATCH: a protein structure and structural feature prediction server.

SCRATCH is a server for predicting protein tertiary structure and structural features. The SCRATCH software suite includes predictors for secondary structure, relative solvent accessibility, disordered regions, domains, disulfide bridges, single mutation stability, residue contacts versus average, individual residue contacts and tertiary structure. The user simply provides an amino acid sequence and selects the desired predictions, then submits to the server. Results are emailed to the user. The server is available at http://www.igb.uci.edu/servers/psss.html.

Disulfides↗

PubMed2004

Tyrosine residues modification studied by MALDI-TOF mass spectrometry.

Amino acid residue-specific reactivity in proteins is of great current interest in structural biology as it provides information about solvent accessibility and reactivity of the residue and, consequently, about protein structure and possible interactions. In the work presented tyrosine residues of three model proteins with known spatial structure are modified with two tyrosine-specific reagents: tetranitromethane and iodine. Modified proteins were specifically digested by proteases and the mass of resulting peptide fragments was determined using matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry. Our results show that there are only small differences in the extent of tyrosine residues modification by tetranitromethane and iodine. However, data dealing with accessibility of reactive residues obtained by chemical modifications are not completely identical with those obtained by nuclear magnetic resonance and X-ray crystallography. These interesting discrepancies can be caused by local molecular dynamics and/or by specific chemical structure of the residues surrounding.

Amino Acid Substitution↗

PubMed2003

Atom depth in protein structure and function.

Atom depth, originally defined as the distance between a protein atom and the nearest water molecule surrounding a protein, is a simple but valuable geometrical descriptor of the protein interior. It can be easily computed from the 3D structure of a protein, thus complementing the information provided by the calculation of the solvent accessible surface area and buried surface area. Depth has been found to be correlated with several molecular, residue and atomic properties, such as average protein domain size, protein stability, free energy of formation of protein complexes, amino acid type hydrophobicity, residue conservation and hydrogen/deuterium amide proton exchange rates.

Algorithms↗

PubMed2004

Recent advances in the development and application of implicit solvent models in biomolecule simulations.

Advances have recently been made in the development of implicit solvent methodologies and their application to the modeling of biomolecules, particularly with regard to generalized Born approaches, dielectric screening function formulations and models based on solvent-accessible surface areas. Interesting new developments include more refined non-polar solvation energy estimators, and implicit methods for modeling low-dielectric and heterogeneous environments such as membrane systems. These have been successfully applied to molecular dynamics simulations, the scoring of protein conformations, and the calculation of binding affinities and folding free energy landscapes.

Algorithms↗

PubMed1990

On the nature of antibody combining sites: unusual structural features that may confer on these sites an enhanced capacity for binding ligands.

A detailed analysis of the structural aspects of antibody-antigen interactions has been made possible by the availability of X-ray structures for three complexes of antilysozyme Fabs to lysozyme (reviewed by Davies et al.: J. Biol. Chem. 263:10541-10544, 1988.) Examination of the antigen-contacting residues in the three antilysozyme Fabs reveals the occurrence of a large number of aromatics, particularly tyrosines, and the absence of apolar, aliphatic residues. Calculation of the frequency of occurrence of the various amino acid types reveals that tyrosines are three times, and histidines and asparagines eight times, more likely to be found in the complementarity-determining regions than in the framework of the variable domains. Analysis of the solvent accessibility of the residues in Fvs (the modules containing variable domains of the light and heavy chains) of known three-dimensional structure indicates that tyrosines and tryptophans are more exposed when they occur in the complementarity-determining regions than when in the framework. Furthermore, many more of the asparagines in the complementarity-determining regions than in the framework are buried. These asparagines appear to have a structural role in that they hydrogen-bond through their side chains to other side chains and, even more so, to the protein backbone. The stabilizing effect of the asparagines, plus the rigidity of the framework, may serve to allow the greater exposure of the aromatic residues to solvent. In view of the greater potential contribution of aromatic side chains to the total binding energy, these results suggest that antibody combining sites have structural features that make them especially suited for interacting with ligands.

Amino Acids↗

PubMed2002

Open channel structure of MscL and the gating mechanism of mechanosensitive channels.

Mechanosensitive channels act as membrane-embedded mechano-electrical switches, opening a large water-filled pore in response to lipid bilayer deformations. This process is critical to the response of living organisms to direct physical stimulation, such as in touch, hearing and osmoregulation. Here, we have determined the structural rearrangements that underlie these events in the large prokaryotic mechanosensitive channel (MscL) using electron paramagnetic resonance spectroscopy and site-directed spin labelling. MscL was trapped in both the open and in an intermediate closed state by modulating bilayer morphology. Transition to the intermediate state is characterized by small movements in the first transmembrane helix (TM1). Subsequent transitions to the open state are accompanied by massive rearrangements in both TM1 and TM2, as shown by large increases in probe dynamics, solvent accessibility and the elimination of all intersubunit spin-spin interactions. The open state is highly dynamic, supporting a water-filled pore of at least 25 A, lined mostly by TM1. These structures suggest a plausible molecular mechanism of gating in mechanosensitive channels.

Computer Simulation↗

PubMed1992

The effect of viscosity on the accessibility of the single tryptophan in human serum albumin.

When reactions take place with one of the reactants tied to protein matrix, movements along the reaction coordinate towards the transition state can become coupled to structural fluctuations of the protein matrix. This investigation aims to test the assumptions underlying the arguments supporting such a coupling. A coupling is allowed only if the activation barrier is high and broad enough as shown to be the case for the proton catalyzed isotope exchange at Trp-63 of lysozyme. In the present investigation the activation barrier for the same reaction has been lowered radically in an effort to show that the coupling, as measured by the dependence of rate on solution viscosity, will diminish and ideally vanish, despite the unchanged effects of cosolvents on the chemical activities of all the reactants. The isotope exchange rate at the indole nitrogen of the single tryptophan residue of human serum albumin was measured with UV. This residue is rigidly held to the protein surface and the solvent access, although restricted, corresponds to a partially exposed residue. As a consequence, the isotope exchange rates and the bimolecular quenching rate of fluorescence by acrylamide, also measured, are high. The experiments were carried out at pH 5.2 where the molecule is in the N-form and the exchange is catalyzed by OH- ions. The activation energy of the hydroxyl catalyzed reaction is 22 kJ lower than for the proton catalyzed process. Under these conditions the exchange rate is viscosity independent both in the case of glycerol and in ethylene glycol.(ABSTRACT TRUNCATED AT 250 WORDS)

Deuterium↗

PubMed1995

Conformation of the circular dumbbell d : structure determination and molecular dynamics.

The circular DNA decamer 5'-d -3' was studied in solution by means of NMR spectroscopy and molecular dynamics in H2O. At a temperature of 269 K, a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4) is present. The L2L2 form contains three Watson-Crick C-G base pairs and two two-residue loops is opposite parts of the molecule. On raising the temperature from 269 K to 314 K, the L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted G(anti)-C(syn) closing base pair in the 5'-GTTC-3' loop with only one remaining (solvent-accessible) hydrogen bond between NH alpha of the cytosine dC(1) and O6 of the guanine dG(8). The opposite 5'-CTTG-3' loop remains stable. The two conformers occur in slow equilibrium (rate constant 2-20 s-1). Structure determination of the L2L2 and L2L4 forms was performed with the aid of a full relaxation matrix approach (IRMA) in combination with restrained MD. Torsional information was obtained from coupling constants. Coupling constant analysis (3JHH, 3JHP, 3JCP) gave detailed information about the local geometry around backbone torsion angles beta, gamma, delta, and epsilon, revealing a relatively high flexibility of the 5'-GTTC-3' loop. The values of the coupling constants are virtually temperature-independent. 'Weakly constrained' molecular dynamics in solvent was used to sample the conformational space of the dumbbell. The relaxation matrices from the MD simulation were averaged over to predict dynamic NOE volumes. In order to account for the 1:1 conformational mixture of L2L2 and L2L4 present at 271 K, we also included S2 factors and averaging of the -averaged relaxation matrices. On matrix averaging, the agreement of NOE volumes with experiment improved significantly for protons located in the thermodynamically less stable 5'-GTTC-3' loop. The difference in stability of the 5'-CTTG-3' and 5'-GTTC-3' loops is mainly caused by differences in the number of potential hydrogen bonds in the minor groove and differences in stacking overlap of the base pairs closing the minihairpin loops. The syn conformation for dC(1), favored at high temperature, is stabilized by solvation in the major groove. However, the conformational properties of the dC(1) base, as deduced from R-factor analysis and MD simulations, include a large flexibility about torsion angle chi.

DNA, Circular↗

PubMed2006

Role of solvation barriers in protein kinetic stability.

The stability of several protein systems of interest has been shown to have a kinetic basis. Besides the obvious biotechnological implications, the general interest of understanding protein kinetic stability is emphasized by the fact that some emerging molecular approaches to the inhibition of amyloidogenesis focus on the increase of the kinetic stability of protein native states. Lipases are among the most important industrial enzymes. Here, we have studied the thermal denaturation of the wild-type form, four single-mutant variants and two highly stable, multiple-mutant variants of lipase from Thermomyces lanuginosa. In all cases, thermal denaturation was irreversible, kinetically controlled and conformed to the two-state irreversible model. This result supports that the novel molecular-dynamics-focused, directed-evolution approach involved in the preparation of the highly stable variants is successful likely because it addresses kinetic stability and, in particular, because heated molecular dynamics simulations possibly identify regions of disrupted native interactions in the transition state for irreversible denaturation. Furthermore, we find very large mutation effects on activation enthalpy and entropy, which were not accompanied by similarly large changes in kinetic urea m-value. From this we are led to conclude that these mutation effects are associated to some structural feature of the transition state for the irreversible denaturation process that is not linked to large changes in solvent accessibility. Recent computational studies have suggested the existence of solvation/desolvation barriers in at least some protein folding/unfolding processes. We thus propose that a solvation barrier (arising from the asynchrony between breaking of internal contacts and water penetration) may contribute to the kinetic stability of lipase from T. lanuginosa (and, possibly, to the kinetic stability of other proteins as well).

Ascomycota↗

PubMed2005

Three-dimensional computation of atom depth in complex molecular structures.

MOTIVATION: For a complex molecular system the delineation of atom-atom contacts, exposed surface and binding sites represents a fundamental step to predict its interaction with solvent, ligands and other molecules. Recently, atom depth has been also considered as an additional structural descriptor to correlate protein structure with folding and functional properties. The distance between an atom and the nearest water molecule or the closest surface dot has been proposed as a measure of the atom depth, but, in both cases, the 3D character of depth is largely lost. In the present study, a new approach is proposed to calculate atom depths in a way that the molecular shape can be taken into account. RESULTS: An algorithm has been developed to calculate intersections between the molecular volume and spheres centered on the atoms whose depth has to be quantified. Many proteins with different size and shape have been chosen to compare the results obtained from distance-based and volume-based depth calculations. From the wealth of experimental data available for hen egg white lysozyme, H/D exchange rates and TEMPOL induced paramagnetic perturbations have been analyzed both in terms of depth indexes and of atom distances to the solvent accessible surface. The algorithm here proposed yields better correlations between experimental data and atom depth, particularly for those atoms which are located near to the protein surface. AVAILABILITY: Instructions to obtain source code and the executable program are available either from http://sienabiografix.com or http://sadic.sourceforge.net CONTACT: niccolai@unisi.it SUPPLEMENTARY INFORMATION: http://www.Sienabiogzefix.com/publication.

Algorithms↗

PubMed2006

Topological and dynamical properties of Azurin anchored to a gold substrate as investigated by molecular dynamics simulation.

A classical molecular dynamics study of the electron transfer protein azurin, covalently bound to a gold substrate through its native disulphide group, is carried out at full hydration. With the aim of investigating the effects on the protein structure and dynamics as induced by the presence of an electric field, simulations are performed on neutral, positively and negatively charged substrates. A number of parameters, such as the average structure, the root mean square deviations and fluctuations, the intraprotein hydrogen bonds and solvent accessible surface of the protein, are monitored during 10 ns of run. The orientation, the height and the lateral size of the protein, with respect to the substrate are evaluated and compared with the experimental data obtained by scanning probe nanoscopies. The electron transfer properties between the copper redox center and the disulphide bridge bound to the substrate are investigated and briefly discussed.

Azurin↗

PubMed2003

Motif refinement of the peroxisomal targeting signal 1 and evaluation of taxon-specific differences.

Eukaryote peroxisomes, plant glyoxysomes and trypanosomal glycosomes belong to the microbody family of organelles that compartmentalise a variety of biochemical processes. The interaction between the PTS1 signal and its cognate receptor Pex5 initiates the major import mechanism for proteins into the matrix of these organelles. Relying on the analysis of amino acid sequence variability of known PTS1-targeted proteins and PTS1-containing peptides that interact with Pex5 in the yeast two-hybrid assay, on binding site studies of the Pex5-ligand complex crystal structure, 3D models and sequences of Pex5 proteins from various taxa, we derived the requirements for a C-terminal amino acid sequence to interact productively with Pex5. We found evidence that, at least the 12 C-terminal residues of a given substrate protein are implicated in PTS1 signal recognition. This motif can be structurally and functionally divided into three regions: (i) the C-terminal tripeptide, (ii) a region interacting with the surface of Pex5 (about four residues further upstream), and (iii) a polar, solvent-accessible and unstructured region with linker function (the remaining five residues). Specificity differences are confined to taxonomic subgroups (metazoa and fungi) and are connected with amino acid type preferences in region 1 and deviating hydrophobicity patterns in region 2.

Amino Acid Motifs↗

PubMed1998

Determinants of protein hydrogen exchange studied in equine cytochrome c.

The exchange of a large number of amide hydrogens in oxidized equine cytochrome c was measured by NMR and compared with structural parameters. Hydrogens known to exchange through local structural fluctuations and through larger unfolding reactions were separately considered. All hydrogens protected from exchange by factors greater than 10(3) are in defined H-bonds, and almost all H-bonded hydrogens including those at the protein surface were measured to exchange slowly. H-exchange rates do not correlate with H-bond strength (length) or crystallographic B factors. It appears that the transient structural fluctuation necessary to bring an exchangeable hydrogen into H-bonding contact with the H-exchange catalyst (OH(-)-ion) involves a fairly large separation of the H-bond donor and acceptor, several angstroms at least, and therefore depends on the relative resistance to distortion of immediately neighboring structure. Accordingly, H-exchange by way of local fluctuational pathways tends to be very slow for hydrogens that are neighbored by tightly anchored structure and for hydrogens that are well buried. The slowing of buried hydrogens may also reflect the need for additional motions that allow solvent access once the protecting H-bond is separated, although it is noteworthy that burial in a protein like cytochrome c does not exceed 4 angstroms. When local fluctuational pathways are very slow, exchange can become dominated by a different category of larger, cooperative, segmental unfolding reactions reaching up to global unfolding.

Animals↗

PubMed1996

Dimer asymmetry in superoxide dismutase studied by molecular dynamics simulation.

Molecular dynamics (MD) simulations of 100 ps have been carried out to study the active-site behaviour of the Cu,Zn superoxide dismutase dimer (SOD) in water. The active site of each subunit was monitored during the whole simulation by calculating the distances between functional residues and the catalytic copper. The results indicate that charge orientation is maintained at each active site but the solvent accessibility varies. Analysis of the MD simulation, carried out by using the atomic displacement covariance matrix, has shown a different intra-subunit correlation pattern for the two monomers and the presence of inter-subunit correlations. The MD simulation presented here indicates an asymmetry in the two active sites and different dynamic behaviour of the two SOD subunits.

Animals↗