Locating missing water molecules in protein cavities by the three-dimensional reference interaction site model theory of molecular solvation

Water molecules confined in protein cavities are of great importance in understanding the protein structure and functions. However, it is a nontrivial task to locate such water molecules in protein by the ordinary molecular simulation and modeling techniques as well as experimental methods. The present study proves that the three-dimensional reference interaction site model (3D-RISM) theory, a recently developed statistical-mechanical theory of molecular solvation, has an outstanding advantage in locating such water molecules. In this paper, we demonstrate that the 3D-RISM theory is able to reproduce the structure and the number of water molecules in cavities of hen egg-white lysozyme observed commonly in the X-ray structures of different resolutions and conditions. Furthermore, we show that the theory successfully identified a water molecule in a cavity, the existence of which has been ambiguous even from the X-ray results. In contrast, we confirmed that molecular dynamics simulation is helpless at present to find such water molecules because the results substantially depend on the initial coordinates of water molecules. Possible applications of the theory to problems in the fields of biochemistry and biophysics are also discussed.     

Theoretical study of the partial molar volume change associated with the pressure-induced structural transition of ubiquitin

The partial molar volume (PMV) change associated with the pressure-induced structural transition of ubiquitin is analyzed by the three-dimensional reference interaction site model (3D-RISM) theory of molecular solvation. The theory predicts that the PMV decreases upon the structural transition, which is consistent with the experimental observation. The volume decomposition analysis demonstrates that the PMV reduction is primarily caused by the decrease in the volume of structural voids in the protein, which is partially canceled by the volume expansion due to the hydration effects. It is found from further analysis that the PMV reduction is ascribed substantially to the penetration of water molecules into a specific part of the protein. Based on the thermodynamic relation, this result implies that the water penetration causes the pressure-induced structural transition. It supports the water penetration model of pressure denaturation of proteins proposed earlier.     

Ion Counting from Explicit-Solvent Simulations and 3D-RISM

The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl_(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20–25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na⁺ binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability.     

Theoretical study of volume changes associated with the helix-coil transition of an alanine-rich peptide in aqueous solution

The changes in the partial molar volume (PMV) associated with the conformational transition of an alanine-rich peptide AK16 from the alpha-helix structure to various random coil structures are calculated by the three-dimensional interaction site model (3D-RISM) theory coupled with the Kirkwood-Buff theory. The volume change is analyzed by decomposing it into contributions from geometry and hydration: the changes in the van der Waals, void, thermal, and interaction volume. The total change in the PMV is positive. This is primarily due to the growth of void space within the peptide, which is canceled in part by the volume reduction resulting from the increase in the electrostatic interaction between the peptide and water molecules. The changes in the void and thermal volume of the coil structures are widely distributed and tend to compensate each other. Additionally, the relations between the hydration volume components and the surface properties are investigated. We categorize coil structures into extended coils with the PMV smaller than helix and general coils with the PMV larger than helix. The pressure therefore can both stabilize and destabilize the coil structures. The latter seems to be a more proper model of random coil structures of the peptide.     

Ion Counting from Explicit-Solvent Simulations and 3D-RISM

The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl_(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20–25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na⁺ binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability.     

THE COLLOIDAL BEHAVIOR OF SERUM GLOBULIN

1. The globulin prepared from ox serum by dilution and precipitation with carbon dioxide has been found, by electrometric titration experiments, to behave like an amphoteric electrolyte, reacting stoichiometrically with acids and bases. 2. The potential difference developed between a solution of globulin chloride, phosphate, or acetate and a solution of the corresponding acid, free from protein, separated from the globulin by a collodion membrane, was found to be influenced by hydrogen ion concentration and salt concentration in the way predicted by Donnan''s theory of membrane equilibrium. In experiments with sodium globulinate and sodium hydroxide it was found that the potential difference could be similarly explained. 3. The osmotic pressure of such solutions could be qualitatively accounted for by the Donnan theory, but exhibited a discrepancy which is explicable by analogy with certain experiments of Loeb on gelatin. 4. The application of Loeb''s theory of colloidal behavior, which had previously been found to hold in the case of gelatin, casein, egg albumin, and edestin, has thus been extended to another protein, serum globulin.     

A simple, continuous-flow ultrafiltration apparatus

A continuous-flow, inverted ultrafiltration apparatus has been made, which permits protein solutions to be concentrated rapidly and has allowed the fractionation of some binary model mixtures. A major problem, which prevents such fractionations in conventional apparatus, has been identified. The design avoids shearing and overheating of labile solutions and appears to be suitable for both micro- and macro-scale applications.     

Theoretical study for volume changes associated with the helix-coil transition of peptides.

We calculate the partial molar volumes and their changes associated with the coil(extended)-to-helix transition of two types of peptide, glycine-oligomer and glutamic acid-oligomer, in aqueous solutions by using the Kirkwood-Buff solution theory coupled with the three-dimensional reference interaction site model (3D-RISM) theory. The volume changes associated with the transition are small and positive. The volume is analyzed by decomposing it into five contributions following the procedure proposed by Chalikian and Breslauer: the ideal volume, the van der Waals volume, the void volume, the thermal volume, and the interaction volume. The ideal volumes and the van der Waals volumes do not change appreciably upon the transition. In the both cases of glycine-peptide and glutamic acid-peptide, the changes in the void volumes are positive, while those in the thermal volumes are negative, and tend to balance those in the void volumes. The change in the interaction volume of glycine-peptide does not significantly contribute, while that of glutamic acid-peptide makes a negative contribution.     

Computational approaches to RNA structure prediction, analysis, and design

RNA molecules are important cellular components involved in many fundamental biological processes. Understanding the mechanisms behind their functions requires RNA tertiary structure knowledge. Although modeling approaches for the study of RNA structures and dynamics lag behind efforts in protein folding, much progress has been achieved in the past two years. Here, we review recent advances in RNA folding algorithms, RNA tertiary motif discovery, applications of graph theory approaches to RNA structure and function, and in silico generation of RNA sequence pools for aptamer design. Advances within each area can be combined to impact many problems in RNA structure and function.     

High hydrostatic pressure dissociates early aggregates of TTR105-115, but not the mature amyloid fibrils

A range of disorders such as Alzheimer's disease and type II diabetes have been linked to protein misfolding and aggregation. Transthyretin is an amyloidogenic protein which is involved in familial amyloid polyneuropathy, the most common form of systemic amyloid disease. A peptide fragment of this protein, TTR105-115, has been shown to form well-defined amyloid fibrils in vitro. In this study, the stability of amyloid fibrils towards high hydrostatic pressure has been investigated by Fourier transform infrared spectroscopy. Information on the morphology of the species exposed to high hydrostatic pressure was obtained by atomic force microscopy. The species formed early in the aggregation process were found to be dissociated by relatively low hydrostatic pressure (220 MPa), whereas mature fibrils are pressure insensitive up to 1.3 GPa. The pressure stability of the mature fibrils is consistent with a fibril structure in which there is an extensive hydrogen bond network in a tightly packed environment from which water is excluded. The fact that early aggregates can be dissociated by low pressure suggests, however, that hydrophobic and electrostatic interactions are the dominant factors stabilizing the species formed in the early stages of fibril formation.     

Protein interactions in 3D: From interface evolution to drug discovery

Over the past 10years, much research has been dedicated to the understanding of protein interactions. Large-scale experiments to elucidate the global structure of protein interaction networks have been complemented by detailed studies of protein interaction interfaces. Understanding the evolution of interfaces allows one to identify convergently evolved interfaces which are evolutionary unrelated but share a few key residues and hence have common binding partners. Understanding interaction interfaces and their evolution is an important basis for pharmaceutical applications in drug discovery. Here, we review the algorithms and databases on 3D protein interactions and discuss in detail applications in interface evolution, drug discovery, and interface prediction.     

Ku, Artemis, and ataxia-telangiectasia-mutated: signalling networks in DNA damage

Cell death linked to DNA damage has been implicated in various diseases caused by environmental stress and infection. Severe DNA damage, which is beyond the capacity of the DNA repair proteins, triggers apoptosis. Accumulation of DNA damage has been proposed to be a principal mechanism of infection, inflammation, cancer, and aging. The most deleterious form of DNA damage is double-strand breaks (DSBs), where ataxia-telangiectasia-mutated (ATM) is the main transducer of the double-strand DNA break signal. Once the DNA is damaged, the DNA repair protein Ku70/80 translocates into the nucleus, a process which may be mediated by ataxia-telangiectasia-mutated, a member of the phosphoinositide-3-kinase-like family. The function and stability of Artemis may also be regulated by ataxia-telangiectasia-mutated through its phosphorylation upon the occurrence of DNA damage. Interestingly, both Artemis and Ku70/80 are substrates of DNA-dependent protein kinase (DNA-PK), another member of the phosphoinositide-3-kinase-like family. In this review, we show how Ku and Artemis function in the DNA damage response and the ataxia-telangiectasia-mutated signaling pathway and discuss potential applications of agents targeting these DNA damage response molecules in the treatment of inflammation and cancer.     

Chiral counter-current chromatography: A survey of its instrument,mechanism, procedure,and applications

The chiral separation by counter-current chromatography has made great progress in the past three decades. It has become increasingly popular in the field of chiral separation, and many applications have been introduced during the last years. This review mainly focuses on the current topics, applications, and trends in chiral separation by counter-current chromatography. It contains the development of modern counter-current chromatography apparatus, theory of counter-current chromatography, overview of applications of chiral counter-current chromatography enantioseparation, its current situation, and challenges. At last, some conclusions and perspectives also have been discussed in this review.     

Genetic engineering of streptavidin,a versatile affinity tag

Streptavidin, a tetrameric protein produced by Streptomyces avidinii, has been used as a useful, versatile affinity tag in a variety of biological applications. The efficacy of streptavidin is derived from its extremely high binding affinity for the vitamin biotin. For the last several years, we have used genetic engineering as a primary means to enhance the properties of streptavidin and to expand the application of streptavidin as an affinity tag. In this review, we describe several genetically engineered streptavidin variants, which include a streptavidin with a reduced biotin-binding affinity, a dimeric streptavidin, and a fusion protein between streptavidin and protein A, along with their potential applications in biological science.     

Structure,aromaticity, stability,and energetic performance of the analogues and derivatives of s-heptazine

The reaction mechanism of diastereoselective oxidative dearomatization by iodoxybenzoic acid of key compounds involved in the total synthesis of epicocconone analogs, which are efficient fluorophores with a wide range of applications in protein staining and separation, was studied using density functional theory. In particular, the conformational space was investigated, as was the role of the so-called hypervalent twist move, which is thought to be the rate-determining step. Both kinetic and thermodynamical aspects of the mechanism were considered from static and dynamic viewpoints, including solvent effects. The results were then rationalized using conceptual density functional theory and Bader’s atoms-in-molecules framework, which demonstrated how complementary these two approaches are when studying organic chemistry reactions theoretically.