Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

PLIC: protein-ligand interaction clusters

Most of the biological processes are governed through specific protein-ligand interactions. Discerning different components that contribute toward a favorable protein-ligand interaction could contribute significantly toward better understanding protein function, rationalizing drug design and obtaining design principles for protein engineering. The Protein Data Bank (PDB) currently hosts the structure of similar to 68 000 protein-ligand complexes. Although several databases exist that classify proteins according to sequence and structure, a mere handful of them annotate and classify protein-ligand interactions and provide information on different attributes of molecular recognition. In this study, an exhaustive comparison of all the biologically relevant ligand-binding sites (84 846 sites) has been conducted using PocketMatch: a rapid, parallel, in-house algorithm. PocketMatch quantifies the similarity between binding sites based on structural descriptors and residue attributes. A similarity network was constructed using binding sites whose PocketMatch scores exceeded a high similarity threshold (0.80). The binding site similarity network was clustered into discrete sets of similar sites using the Markov clustering (MCL) algorithm. Furthermore, various computational tools have been used to study different attributes of interactions within the individual clusters. The attributes can be roughly divided into (i) binding site characteristics including pocket shape, nature of residues and interaction profiles with different kinds of atomic probes, (ii) atomic contacts consisting of various types of polar, hydrophobic and aromatic contacts along with binding site water molecules that could play crucial roles in protein-ligand interactions and (iii) binding energetics involved in interactions derived from scoring functions developed for docking. For each ligand-binding site in each protein in the PDB, site similarity information, clusters they belong to and description of site attributes are provided as a relational database-protein-ligand interaction clusters (PLIC).

Tryptophan-water interaction in Monellin:Hydration patterns from molecular dynamics simulations

Femtosecond spectroscopy carried out earlier on Monellin and some other systems has given insights into the hydration dynamics of the proteins. In the present work, molecular dynamics simulations have been performed on Monellin to study the hydration dynamics. A method has been described to follow up the molecular events of the protein–water interactions in detail. The time constants of the survival correlation function match well with the reported experimental values. This validates the procedure, adapted here for Monellin, to investigate the hydration dynamics in general.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

A regulator of G protein signaling interaction surface linked to effector specificity

Proteins of the regulator of G protein signaling (RGS) family accelerate GTP hydrolysis by the α subunits (Gα) of G proteins, leading to rapid recovery of signaling cascades. Many different RGS proteins can accelerate GTP hydrolysis by an individual Gα, and GTP hydrolysis rates of different Gαs can be enhanced by the same RGS protein. Consequently, the mechanisms for specificity in RGS regulation and the residues involved remain unclear. Using the evolutionary trace (ET) method, we have identified a cluster of residues in the RGS domain that includes the RGS-Gα binding interface and extends to include additional functionally important residues on the surface. One of these is within helix α3, two are in α5, and three are in the loop connecting α5 and α6. A cluster of surface residues on Gα previously identified by ET, and composed predominantly of residues from the switch III region and helix α3, is spatially contiguous with the ET-identified residues in the RGS domain. This cluster includes residues proposed to interact with the γ subunit of Gtα's effector, cGMP phosphodiesterase (PDEγ). The proximity of these clusters suggests that they form part of an interface between the effector and the RGS-Gα complex. Sequence variations in these residues correlate with PDEγ effects on GTPase acceleration. Because ET identifies residues important for all members of a protein family, these residues likely form a general site for regulation of G protein-coupled signaling cascades, possibly by means of effector interactions.

Packing at the protein-water interface.

We have determined the packing efficiency at the protein-water interface by calculating the volumes of atoms on the protein surface and nearby water molecules in 22 crystal structures. We find that an atom on the protein surface occupies, on average, a volume approximately 7% larger than an atom of equivalent chemical type in the protein core. In these calculations, larger volumes result from voids between atoms and thus imply a looser or less efficient packing. We further find that the volumes of individual atoms are not related to their chemical type but rather to their structural location. More exposed atoms have larger volumes. Moreover, the packing around atoms in locally concave, grooved regions of protein surfaces is looser than that around atoms in locally convex, ridge regions. This as a direct manifestation of surface curvature-dependent hydration. The net volume increase for atoms on the protein surface is compensated by volume decreases in water molecules near the surface. These waters occupy volumes smaller than those in the bulk solvent by up to 20%; the precise amount of this decrease is directly related to the extent of contact with the protein.

Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

Nature of the thiamin-binding protein from chicken egg yolk

A simple, rapid and efficient procedure for the purification of thiamin-binding protein from chicken egg yolk was developed. The method involved removal, by exclusion, of lipoproteins from DEAE-cellulose and subsequent elution of water-soluble proteins held on the ion-exchanger with 1 M-NaCl, followed by treatment of the eluted protein fraction with an aqueous suspension of dextran/charcoal to generate apoprotein from the holoprotein. The resultant protein fraction was subjected to bioaffinity chromatography on thiamin pyrophosphate--AE (aminoethyl)-Sepharose. The protein eluted specifically with 10 microM-thiamin at pH 7.0, was homogeneous by the criteria of polyacrylamide-gel disc electrophoresis, had a mol.wt. of 38 000 +/- 2000 and was not a glycoprotein. The purified thiamin-binding protein specifically interacted with riboflavin-binding protein with no detectable deleterious affect on its (14C)thiamin-binding capacity. The protein bound [14C]thiamin with a molar ratio of 1.0, with dissociation constant (Kd) 0.41 microM. This protein-ligand interaction was inhibited by thiamin analogues and antagonists. The absorption spectrum of the protein in the presence of thiamin exhibited significant hypochromism at the 278 nm band, indicating the involvement of aromatic amino acid residues of the protein, during its binding to the ligand. The protein cross-reacted with the monospecific antiserum to egg-white thiamin-binding protein, showing thereby that thiamin-binding proteins present in chicken egg yolk and white are the products of the same structural gene.

Site-specific N-Linked Glycosylation of Receptor Guanylyl Cyclase C Regulates Ligand Binding, Ligand-mediated Activation and Interaction with Vesicular Integral Membrane Protein 36, VIP36

Guanylyl cyclase C (GC-C) is a multidomain, membrane-associated receptor guanylyl cyclase. GC-C is primarily expressed in the gastrointestinal tract, where it mediates fluid-ion homeostasis, intestinal inflammation, and cell proliferation in a cGMP-dependent manner, following activation by its ligands guanylin, uroguanylin, or the heat-stable enterotoxin peptide (ST). GC-C is also expressed in neurons, where it plays a role in satiation and attention deficiency/hyperactive behavior. GC-C is glycosylated in the extracellular domain, and differentially glycosylated forms that are resident in the endoplasmic reticulum (130 kDa) and the plasma membrane (145 kDa) bind the ST peptide with equal affinity. When glycosylation of human GC-C was prevented, either by pharmacological intervention or by mutation of all of the 10 predicted glycosylation sites, ST binding and surface localization was abolished. Systematic mutagenesis of each of the 10 sites of glycosylation in GC-C, either singly or in combination, identified two sites that were critical for ligand binding and two that regulated ST-mediated activation. We also show that GC-C is the first identified receptor client of the lectin chaperone vesicular integral membrane protein, VIP36. Interaction with VIP36 is dependent on glycosylation at the same sites that allow GC-C to fold and bind ligand. Because glycosylation of proteins is altered in many diseases and in a tissue-dependent manner, the activity and/or glycan-mediated interactions of GC-C may have a crucial role to play in its functions in different cell types.

Function of the hydration layer around an antifreeze protein revealed by atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations are used to investigate the mechanism by which the antifreeze protein from the spruce budworm, Choristoneura fumiferana, binds to ice. Comparison of structural and dynamic properties of the water around the three faces of the triangular prism-shaped protein in aqueous solution reveals that at low temperature the water structure is ordered and the dynamics slowed down around the ice-binding face of the protein, with a disordering effect observed around the other two faces. These results suggest a dual role for the solvation water around the protein. The preconfigured solvation shell around the ice-binding face is involved in the initial recognition and binding of the antifreeze protein to ice by lowering the barrier for binding and consolidation of the protein:ice interaction surface. Thus, the antifreeze protein can bind to the molecularly rough ice surface by becoming actively involved in the formation of its own binding site. Also, the disruption of water structure around the rest of the protein helps prevent the adsorbed protein becoming covered by further ice growth.

Amino-terminal protein-protein interaction motif (POZ-domain) is responsible for activities of the promyelocytic leukemia zinc finger-retinoic acid receptor-alpha fusion protein.

Promyelocytic leukemia zinc finger-retinoic acid receptor a (PLZF-RARalpha), a fusion receptor generated as a result of a variant t(11;17) chromosomal translocation that occurs in a small subset of acute promyelocytic leukemia (APL) patients, has been shown to display a dominant-negative effect against the wild-type RARalpha/retinoid X receptor alpha (RXRalpha). We now show that its N-terminal region (called the POZ-domain), which mediates protein-protein interaction as well as specific nuclear localization of the wild-type PLZF and chimeric PLZF-RARalpha proteins, is primarily responsible for this activity. To further investigate the mechanisms of PLZF-RARalpha action, we have also studied its ligand-receptor, protein-protein, and protein-DNA interaction properties and compared them with those of the promyelocytic leukemia gene (PML)-RARalpha, which is expressed in the majority of APLs as a result of t(15;17) translocation. PLZF-RARalpha and PML-RARalpha have essentially the same ligand-binding affinities and can bind in vitro to retinoic acid response elements (RAREs) as homodimers or heterodimers with RXRalpha. PLZF-RARalpha homodimerization and heterodimerization with RXRalpha were primarily mediated by the POZ-domain and RARalpha sequence, respectively. Despite having identical RARalpha sequences, PLZF-RARalpha and PML-RARalpha homodimers recognized with different affinities distinct RAREs. Furthermore, PLZF-RARalpha could heterodimerize in vitro with the wild-type PLZF, suggesting that it may play a role in leukemogenesis by antagonizing actions of not only the retinoid receptors but also the wild-type PLZF and possibly other POZ-domain-containing regulators. These different protein-protein interactions and the target gene specificities of PLZF-RARalpha and PML-RARalpha may underlie, at least in part, the apparent resistance of APL with t(11;17) to differentiation effects of all-trans-retinoic acid.

Quantitative description of the interaction between folate and the folate-binding protein from cow's milk

A detailed study has been carried out on the dependence of folate binding on the concentration of FBP (folate-binding protein) at pH 5.0, conditions selected to prevent complications arising from the pre-existing self-association of the acceptor. In contrast with the mandatory requirement that reversible interaction of ligand with a single acceptor site should exhibit a unique, rectangular hyperbolic binding curve, results obtained by ultrafiltration for the FBP-folate system required description in terms of (i) a sigmoidal relationship between concentrations of bound and free folate and (ii) an inverse dependence of affinity on FBP concentration. These findings have been attributed to the difficulties in determining the free ligand concentration in the FBP-folate mixtures for which reaction is essentially stoichiometric. This explanation also accounts for the similar published behaviour of the FBP-folate system at neutral pH, which had been attributed erroneously to acceptor self-association, a phenomenon incompatible with the experimental findings because of its prediction of a greater affinity for folate with increasing FBP concentration.