Conservation of inter-residue interactions and prediction of folding rates of domain repeats

Domains are the main structural and functional units of larger proteins. They tend to be contiguous in primary structure and can fold and function independently. It has been observed that 10–20% of all encoded proteins contain duplicated domains and the average pairwise sequence identity between them is usually low. In the present study, we have analyzed the structural similarity between domain repeats of proteins with known structures available in the Protein Data Bank using structure-based inter-residue interaction measures such as the number of long-range contacts, surrounding hydrophobicity, and pairwise interaction energy. We used RADAR program for detecting the repeats in a protein sequence which were further validated using Pfam domain assignments. The sequence identity between the repeats in domains ranges from 20 to 40% and their secondary structural elements are well conserved. The number of long-range contacts, surrounding hydrophobicity calculations and pairwise interaction energy of the domain repeats clearly reveal the conservation of 3-D structure environment in the repeats of domains. The proportions of mainchain–mainchain hydrogen bonds and hydrophobic interactions are also highly conserved between the repeats. The present study has suggested that the computation of these structure-based parameters will give better clues about the tertiary environment of the repeats in domains. The folding rates of individual domains in the repeats predicted using the long-range order parameter indicate that the predicted folding rates correlate well with most of the experimentally observed folding rates for the analyzed independently folded domains.     

Evaluation of direct and cooperative contributions towards the strength of buried hydrogen bonds and salt bridges

An experimental approach to evaluate the net binding free energy of buried hydrogen bonds and salt bridges is presented. The approach, which involves a modified multiple-mutant cycle protocol, was applied to selected interactions between TEM-1-beta-lactamase and its protein inhibitor, BLIP. The selected interactions (two salt bridges and two hydrogen bonds) all involving BLIP-D49, define a distinct binding unit. The penta mutant, where all side-chains constructing the binding unit were mutated to Ala, was used as a reference state to which combinations of side-chains were introduced. At first, pairs of interacting residues were added allowing the determination of interaction energies in the absence of neighbors, using double mutant cycles. Addition of neighboring residues allowed the evaluation of their cooperative effects on the interaction. The two isolated salt bridges were either neutral or repulsive whereas the two hydrogen bonds contribute 0.3 kcal mol(-1 )each. Conversely, a double mutant cycle analysis of these interactions in their native environment showed that they all stabilize the complex by 1-1.5 kcal mol(-1). Examination of the effects of neighboring residues on each of the interactions revealed that the formation of a salt bridge triad, which involves two connected salt bridges, had a strong cooperative effect on stabilizing the complex independent of the presence or absence of additional neighbors. These results demonstrate the importance of forming net-works of buried salt bridges. We present theoretical electrostatic calculations which predict the observed mode of cooperativity, and suggest that the cooperative networking effect results from the favorable contribution of the protein to the interaction. Furthermore, a good correlation between calculated and experimentally determined interaction energies for the two salt bridges, and to a lesser extent for the two hydrogen bonds, is shown. The data analysis was performed on values of DeltaDeltaG(double dagger)K(d) which reflect the strength of short range interactions, while DeltaDeltaG(o)K(D) values which include the effects of long range electrostatic forces that alter specifically DeltaDeltaG(double dagger)k(a) were treated separately.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Structure of the tandem MA-3 region of Pdcd4 protein and characterization of its interactions with eIF4A and eIF4G: molecular mechanisms of a tumor suppressor

One of the key regulatory points of translation initiation is recruitment of the 43S preinitation complex to the 5' mRNA cap by the eIF4F complex (eIF4A, eIF4E, and eIF4G). The tumor suppressor protein Pdcd4 has been shown to inhibit cap-dependent translation by interacting tightly with the RNA helicase eIF4A via its tandem MA-3 domains. The NMR studies reported here reveal a fairly extensive and well defined interface between the two MA-3 domains in solution, which appears to be stabilized by a network of interdomain salt bridges and hydrogen bonds, and reveals a unique orientation of the two domains. Characterization of the stoichiometry of the Pdcd4-eIF4A complex suggests that under physiological conditions Pdcd4 binds to a single molecule of eIF4A, which involves contacts with both Pdcd4 MA-3 domains. We also show that contacts mediated by a conserved acidic patch on the middle MA-3 domain of Pdcd4 are essential for forming a tight complex with eIF4A in vivo, whereas the equivalent region of the C-terminal MA-3 domain appears to have no role in complex formation in vivo. The formation of a 1:1 eIF4A-Pdcd4 complex in solution is consistent with the reported presence in vivo of only one molecule of eIF4A in the eIF4F complex. Pdcd4 has also been reported to interact directly with the middle region of eIF4G, however, we were unable to obtain any evidence for even a weak, transient direct interaction.     

Distribution of charges in Bacillus intermedius 7P ribonuclease determines the number of cooperatively melting regions of the globule

A correlation between the distribution of charged side groups in the globule of Bacillus intermedius 7P ribonuclease (binase) and the process of heat denaturation was studied at different pH values in order to estimate a relation between charge distribution in globular proteins and the character of cooperative thermodynamic transitions. As was shown by comparing the results of scanning microcalorimetric analysis of heat denaturation with the three-dimensional structure of binase, at optimal pH the molecule exists as a single cooperative system stabilized by hydrogen bonds, Van der Waals' contacts, and electrostatic interactions like salt bridges. At pH lower than 4.0 (below the physiological optimum) the cooperativity type of the system was found to change due to a reversible cooperative transition in the ternary structure of the protein globule. It has been concluded that the molecular architecture and the arrangement of atoms do not change considerably in different environments; thus the thermodynamic properties of the globule vary due to the alteration of charge distribution and the consequent changes in the size and number of cooperative regions of the globule. Thus, structural and energetic domains may be non-coincident in proteins.     

Crystal structure of Leishmania major oligopeptidase B gives insight into the enzymatic properties of a trypanosomatid virulence factor

Oligopeptidase B (OPB) is a serine peptidase with dibasic substrate specificity. It is found in bacteria, plants, and trypanosomatid pathogens, where it has been identified as a virulence factor and potential drug target. In this study we expressed active recombinant Leishmania major OPB and provide the first structure of an oligopeptidase B at high resolution. The crystallographic study reveals that OPB comprises two domains, a catalytic and a propeller domain, linked together by a hinge region. The structure has been determined in complex with the oligopeptide, protease-inhibitor antipain, giving detailed information on the enzyme active site and extended substrate binding pockets. It shows that Glu-621 plays a critical role in the S1 binding pocket and, along with Phe-603, is largely responsible for the enzyme substrate specificity in P1. In the S2 binding pocket, Tyr-499 was shown to be important for substrate stability. The structure also allowed an investigation into the function of residues highlighted in other studies including Glu-623, which was predicted to be involved in the S1 binding pocket but is found forming an inter-domain hydrogen bond. Additional important salt bridges/hydrogen bonds between the two domains were observed, highlighting the significance of the domain interface in OPB. This work provides a foundation for the study of the role of OPBs as virulence factors in trypanosomatids. It could facilitate the development of specific OPB inhibitors with therapeutic potential by exploiting its unique substrate recognition properties as well as providing a model for OPBs in general.     

Probing the side chain tolerance for inhibitors of the CD95/PLCγ1 interaction

Proceeding our effort to study protein-protein interaction between the death receptor CD95 and phospholipase PLCγ1, we present in the current work chameleon-like traits of peptidomimetic inhibitors.Minute analysis of the interaction suggests that most of the binding energy relies on van der Waals contacts rather than more specific features, such as hydrogen bonds or salt bridges. The two most important positions of the peptoid for its interaction with PLCγ1 (Arg184 and Arg187) were modified to test this hypothesis. While Arg184 proves to be exchangeable for Trp, with no alteration in affinity, the nature of the amino acid replacing Arg187 is more dependent on its positive charge. However, affinity can be partially recovered by increasing van der Waals interactions. Overall, this study shows that for both positions, a subtle balance exists between hydrophobicity, surface contacts and affinity for CD95/PLCγ1, and provides information for the generation of new therapeutic compounds toward this druggable target.     

Distribution of Charges in Bacillus intermedins 7P Ribonuclease Determines the Number of Cooperatively Melting Regions of the Globule

Abstract

A correlation between the distribution of charged side groups in the globule of Bacillus intermedius 7P ribonuclease (binase) and the process of heat denaturation was studied at different pH values in order to estimate a relation between charge distribution in globular proteins and the character of cooperative thermodynamic transitions. As was shown by comparing the results of scanning microcalorimetric analysis of heat denaturation with the three-dimentional structure of binase, at optimal pH the molecule exists as a single cooperative system stabilized by hydrogen bonds, Van der Waals' contacts, and electrostatic interactions like salt bridges. At pH lower than 4.0 (below the physiological optimum) the cooperativity type of the system was found to change due to a reversible cooperative transition in the ternary structure of the protein globule. It has been concluded that the molecular architecture and the arrangement of atoms do not change considerably in different environments; thus the thermodynamic properties of the globule vary due to the alteration of charge distribution and the consequent changes in the size and number of cooperative regions of the globule. Thus, structural and energetic domains may be non-coincident in proteins.     

Forces driving the binding of homeodomains to DNA

Homeodomains are helix-turn-helix type DNA-binding domains that exhibit sequence-specific DNA binding by insertion of their "recognition" alpha helices into the major groove and a short N-terminal arm into the adjacent minor groove without inducing substantial distortion of the DNA. The stability and DNA binding of four representatives of this family, MATalpha2, engrailed, Antennapedia, and NK-2, and truncated forms of the last two lacking their N-terminal arms have been studied by a combination of optical and microcalorimetric methods at different temperatures and salt concentrations. It was found that the stability of the free homeodomains in solution is rather low and, surprisingly, is reduced by the presence of the N-terminal arm for the Antennapedia and NK-2 domains. Their stabilities depend significantly upon the presence of salt: strongly for NaCl but less so for NaF, demonstrating specific interactions with chloride ions. The enthalpies of association of the homeodomains with their cognate DNAs are negative, at 20 degrees C varying only between -12 and -26 kJ/mol for the intact homeodomains, and the entropies of association are positive; i.e., DNA binding is both enthalpy- and entropy-driven. Analysis of the salt dependence of the association constants showed that the electrostatic component of the Gibbs energy of association resulting from the entropy of mixing of released ions dominates the binding, being about twice the magnitude of the nonelectrostatic component that results from dehydration of the protein/DNA interface, van der Waals interactions, and hydrogen bonding. A comparison of the effects of NaCl/KCl with NaF showed that homeodomain binding results in a release not only of cations from the DNA phosphates but also of chloride ions specifically associated with the proteins. The binding of the basic N-terminal arms in the minor groove is entirely enthalpic with a negative heat capacity effect, i.e., is due to sequence-specific formation of hydrogen bonds and hydrophobic interactions rather than electrostatic contacts with the DNA phosphates.     

Crystal structure of the specific DNA-binding domain of Tc3 transposase of C.elegans in complex with transposon DNA.

The crystal structure of the complex between the N-terminal DNA-binding domain of Tc3 transposase and an oligomer of transposon DNA has been determined. The specific DNA-binding domain contains three alpha-helices, of which two form a helix-turn-helix (HTH) motif. The recognition of transposon DNA by the transposase is mediated through base-specific contacts and complementarity between protein and sequence-dependent deformations of the DNA. The HTH motif makes four base-specific contacts with the major groove, and the N-terminus makes three base-specific contacts with the minor groove. The DNA oligomer adopts a non-linear B-DNA conformation, made possible by a stretch of seven G:C base pairs at one end and a TATA sequence towards the other end. Extensive contacts (seven salt bridges and 16 hydrogen bonds) of the protein with the DNA backbone allow the protein to probe and recognize the sequence-dependent DNA deformation. The DNA-binding domain forms a dimer in the crystals. Each monomer binds a separate transposon end, implying that the dimer plays a role in synapsis, necessary for the simultaneous cleavage of both transposon termini.     

Interdomain salt bridges modulate ligand-induced domain motion of the sulfate receptor protein for active transport.

The refined crystal structure of the liganded form of the Salmonella typhimurium sulfate-binding protein, a periplasmic receptor of active transport, is made up of two globular domains bisected by a deep cleft wherein the dehydrated sulfate is completely engulfed and bound by hydrogen bonds and van der Waals' forces. Two salt bridges (between Glu15 and Arg174 and between Asp68 and Arg134) span the cleft opening. To elucidate the role of the inter-domain salt bridges in the ligand-induced domain motion, the acidic residues were changed (singly and together) to their corresponding amide side-chains by site-directed mutagenesis of the recombinant Escherichia coli sulfate-binding protein. Rapid kinetics and equilibrium measurements of sulfate binding to the purified mutant proteins demonstrate that these salt bridges stabilize the closed liganded form of the receptor and modulate the rate of cleft opening. Our results have new implications in understanding the dynamics of many other multidomain proteins that undergo similar large-scale domain motions.     

Interdomain interactions in hinge-bending transitions.

BACKGROUND: The mechanisms that allow or constrain protein movement have not been understood. Here we study interdomain interactions in proteins to investigate hinge-bending motions. RESULTS: We find a limited number of salt bridges and hydrogen bonds at the interdomain interface, in both the "closed" and the "open" conformations. Consistently, analysis of 222 salt bridges in an independently selected database indicates that most salt bridges form within rather than between independently folding hydrophobic units. Calculations show that these interdomain salt bridges either destabilize or only marginally stabilize the closed conformation in most proteins. In contrast, the nonpolar buried surface area between the moving parts can be extensive in the closed conformations. However, when the nonpolar buried surface area is large, we find that at the interdomain interface in the open conformation it may be as large or larger than in the closed conformation. Hence, the energetic penalty of opening the closed conformation is overcome. Consistently, a large nonpolar surface area buried in the closed interdomain interface accompanies limited opening of the domains, yielding a larger interface. CONCLUSIONS: Short-range electrostatic interactions are largely absent between moving domains. Interdomain nonpolar buried surface area may be large in the closed conformation, but it is largely offset by the area buried in the open conformation. In such cases the opening of the domains appears to be relatively small. This may allow prediction of the extent of domain opening. Such predictions may have implications for the shape and size of the binding pockets in drug/protein design.     

Stability and rigidity/flexibility—Two sides of the same coin?

Protein molecules require both flexibility and rigidity for functioning. The fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. We have determined flexible regions for four homologous pairs from thermophilic and mesophilic organisms by two methods: the fast FoldUnfold which uses amino acid sequence and the time consuming MDFirst which uses three-dimensional structures. We demonstrate that both methods allow determining flexible regions in protein structure. For three of the four thermophile–mesophile pairs of proteins, FoldUnfold predicts practically the same flexible regions which have been found by the MD/First method. As expected, molecular dynamics simulations show that thermophilic proteins are more rigid in comparison to their mesophilic homologues. Analysis of rigid clusters and their decomposition provides new insights into protein stability. It has been found that the local networks of salt bridges and hydrogen bonds in thermophiles render their structure more stable with respect to fluctuations of individual contacts. Such network includes salt bridge triads Agr-Glu-Lys and Arg-Glu-Arg, or salt bridges (such as Arg-Glu) connected with hydrogen bonds. This ionic network connects alpha helices and rigidifies the structure. Mesophiles can be characterized by stand alone salt bridges and hydrogen bonds or small ionic clusters. Such difference in the network of salt bridges results in different flexibility of homologous proteins. Combining both approaches allows characterizing structural features in atomic detail that determine the rigidity/flexibility of a protein structure. This article is a part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

The Ras G Domain Lacks the Intrinsic Propensity to Form Dimers

Ras GTPase is a molecular switch controlling a number of cellular pathways including growth, proliferation, differentiation, and apoptosis. Recent reports indicated that Ras undergoes dimerization at the membrane surface through protein-protein interactions. If firmly established this property of Ras would require profound reassessment of a large amount of published data and modification of the Ras signaling paradigm. One proposed mechanism of dimerization involves formation of salt bridges between the two GTPase domains (G domains) leading to formation of a compact dimer as observed in Ras crystal structures. In this work, we interrogated the intrinsic ability of Ras to self-associate in solution by creating conditions of high local concentration through irreversibly tethering the two G domains together at their unstructured C-terminal tails. We evaluated possible self-association in this inverted tandem conjugate via analysis of the time-domain fluorescence anisotropy and NMR chemical shift perturbations. We did not observe the increased rotational correlation time expected for the G domain dimer. Variation of the ionic strength (to modulate stability of the salt bridges) did not affect the rotational correlation time in the tandem further supporting independent rotational diffusion of two G domains. In a parallel line of experiments to detect and map weak self-association of the G domains, we analyzed NMR chemical shifts perturbations at a number of sites near the crystallographic dimer interface. The nearly complete lack of chemical shift perturbations in the tandem construct supported a simple model with the independent G domains repelled from each other by their overall negative charge. These results lead us to the conclusion that self-association of the G domains cannot be responsible for homodimerization of Ras reported in the literature.     

Structural basis of the LOV1 dimerization of Arabidopsis phototropins 1 and 2

Phototropin (phot) is a blue-light receptor protein that triggers phototropic responses, chloroplast relocation, and stomata opening to maximize the efficiency of photosynthesis in higher plants. Phot is composed of three functional domains. The N-terminal half folds into two light-oxygen-voltage-sensing domains called LOV1 and LOV2, each binding a flavin mononucleotide to absorb blue light. The C-terminal half is a serine/threonine kinase domain that causes light-dependent autophosphorylation leading to cellular signaling cascades. LOV2 domain is primarily responsible for activation of the kinase, and LOV1 domain is thought to act as a dimerization site and to regulate sensitivity to activation by blue light. Here we show the crystal structures of LOV1 domains of Arabidopsis phot1 and phot2 in the dark at resolutions of 2.1 Å and 2.0 Å, respectively. Either LOV1 domain forms a dimer through face-to-face association of β-scaffolds in the crystallographic asymmetric unit. Three types of interactions stabilizing the dimer structures found are as follows: contacts of side chains in their β-scaffolds, hydrophobic interactions of a short helix found in the N-terminus of a subunit with the β-scaffolds of both subunits, and hydrogen bonds mediated by hydration water molecules filling the dimer interface. The critical residues for dimerization are Cys261, forming a disulfide bridge between subunits in phot1-LOV1 domain, and Thr217 and Met232 in phot2-LOV1. The topology in homodimeric associations of the LOV1 domains is discussed when referring to those of homodimers or heterodimers of light-oxygen-voltage-sensing or Per-ARNT-Sim domains. The present results also provide clues to understanding structural basis in dimeric interactions of Per-ARNT-Sim protein modules in cellular signaling.     

Structural basis of the Na+/H+ exchanger regulatory factor PDZ1 interaction with the carboxyl-terminal region of the cystic fibrosis transmembrane conductance regulator

The PDZ1 domain of the Na(+)/H(+) exchanger regulatory factor (NHERF) binds with nanomolar affinity to the carboxyl-terminal sequence QDTRL of the cystic fibrosis transmembrane conductance regulator (CFTR) and plays a central role in the cellular localization and physiological regulation of this chloride channel. The crystal structure of human NHERF PDZ1 bound to the carboxyl-terminal peptide QDTRL has been determined at 1.7-A resolution. The structure reveals the specificity and affinity determinants of the PDZ1-CFTR interaction and provides insights into carboxyl-terminal leucine recognition by class I PDZ domains. The peptide ligand inserts into the PDZ1 binding pocket forming an additional antiparallel beta-strand to the PDZ1 beta-sheet, and an extensive network of hydrogen bonds and hydrophobic interactions stabilize the complex. Remarkably, the guanido group of arginine at position -1 of the CFTR peptide forms two salt bridges and two hydrogen bonds with PDZ1 residues Glu(43) and Asn(22), respectively, providing the structural basis for the contribution of the penultimate amino acid of the peptide ligand to the affinity of the interaction.     

Electrostatic contributions drive the interaction between Staphylococcus aureus protein Efb-C and its complement target C3d

The C3-inhibitory domain of Staphylococcus aureus extracellular fibrinogen-binding protein (Efb-C) defines a novel three-helix bundle motif that regulates complement activation. Previous crystallographic studies of Efb-C bound to its cognate subdomain of human C3 (C3d) identified Arg-131 and Asn-138 of Efb-C as key residues for its activity. In order to characterize more completely the physical and chemical driving forces behind this important interaction, we employed in this study a combination of structural, biophysical, and computational methods to analyze the interaction of C3d with Efb-C and the single-point mutants R131A and N138A. Our results show that while these mutations do not drastically affect the structure of the Efb-C/C3d recognition complex, they have significant adverse effects on both the thermodynamic and kinetic profiles of the resulting complexes. We also characterized other key interactions along the Efb-C/C3d binding interface and found an intricate network of salt bridges and hydrogen bonds that anchor Efb-C to C3d, resulting in its potent complement inhibitory properties.     

Understanding Thermostability Factors of <Emphasis Type="Italic">Aspergillus niger</Emphasis> PhyA Phytase: A Molecular Dynamics Study

Molecular dynamics simulation was used to study the dynamic differences between native Aspergillus niger PhyA phytase and a mutant with 20 % greater thermostability. Atomic root mean square deviation, radius of gyration, and number of hydrogen bonds and salt bridges are examined to determine thermostability factors. The results suggest that, among secondary structure elements, loops have the most impact on the thermal stability of A. niger phytase. In addition, the location rather than the number of hydrogen bonds is found to have an important contribution to thermostability. The results also show that salt bridges may have stabilizing or destabilizing effect on the enzyme and influence its thermostability accordingly.     

Significant role of electrostatic interactions for stabilization of protein assemblies

Contribution of electrostatic interactions to stability of BPTI orthorhombic, pig-insulin cubic crystals, and horse L ferritin crystals was evaluated with numerical calculation of Poisson-Boltzmann equation based on a dielectric model. The stability of a ferritin molecule (24-mer) composed of 24 subunits was also evaluated. It was found that the surface charge-charge interactions at separation distances (< 5 Å) were insensitive to variations in the ionic strength, and thus stabilized assembled states of the proteins (i.e., crystalline state and oligomeric state). It was also revealed that the charge density and the packing of the protein crystals were largely responsible for the ionic strength dependence of the crystal stability. The stability of the 5PTI crystalline state with a high charge density drastically increased as the concentration of the solvent ions increased. In contrast, that of the insulin crystal with a low charge density and large solvent region was insensitive to changes in the ionic concentration. The electrostatic interaction between ferritin 24-mers was attributed to two salt bridges mediated by Cd ion. For the stability of the ferritin 24-mer, which is evolutionally designed, the electrostatic stabilization between the subunits was attributed to polar bonds such as buried salt bridges or hydrogen bonds, which occasionally yielded more than 5 kcal/mol and were numerous and very strong compared with the bonds between molecules in the 5PTI and 9INS crystals.By analyzing the atomic charge-charge interactions in detail, it was found that charge pairs separated by less than 3 Å, such as hydrogen bonds, dominantly stabilize the assembled states, and that pairs 3 to 5 Å apart were also important. The stability of the assembled states evaluated by the total EET was determined by the fine balance between the two competing contributions arising from the stabilizing atoms and the destabilizing atoms.Changes of the ASA and hydration free energy were also evaluated in accordance with the process of the subunit assembly. The change of hydration free energy, which was very large (i.e., ~+ 100 kcal/mol/subunit) and unfavorable for the assembly, was proportional to the electrostatic hydration energy (i.e., Born energy change in the hydration process). Hydrophobic groups were likely to appear more frequently than hydrophilic groups at the interfaces.This study offers a method which can improve the stability of protein crystals by introducing polar or charged residues that are properly designed to form specific hydrogen bonds or salt bridges between neighboring protein molecules. This method is also applicable to crystallography, because it improves refinement of protein structures in crystals by taking the inter-protein interactions into account.     

Hydration of Nucleic Acid Bases Studied Using Novel Atom-Atom Potential Functions

Abstract

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contain terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases—uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecule interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interaction play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by “bridges” formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.