Theoretical investigation on the glycan‐binding specificity of Agrocybe cylindracea galectin using molecular modeling and molecular dynamics simulation studies

Galectins are β‐galactoside binding proteins which have the ability to serve as potent antitumor, cancer biomarker, and induce tumor cell apoptosis. Agrocybe cylindracea galectin (ACG) is a fungal galectin which specifically recognizes α(2,3)‐linked sialyllactose at the cell surface that plays extensive roles in the biological recognition processes. To investigate the change in glycan‐binding specificity upon mutations, single point and double point site‐directed in silico mutations are performed at the binding pocket of ACG. Molecular dynamics (MD) simulation studies are carried out for the wild‐type (ACG) and single point (ACG1) and double point (ACG2) mutated ACGs to investigate the dynamics of substituted mutants and their interactions with the receptor sialyllactose. Plausible binding modes are proposed for galectin–sialylglycan complexes based on the analysis of hydrogen bonding interactions, total pair‐wise interaction energy between the interacting binding site residues and sialyllactose and binding free energy of the complexes using molecular mechanics–Poisson–Boltzmann surface area. Our result shows that high contribution to the binding in different modes is due to the direct and water‐mediated hydrogen bonds. The binding specificity of double point mutant Y59R/N140Q of ACG2 is found to be high, and it has 26 direct and water‐mediated hydrogen bonds with a relatively low‐binding free energy of −47.52 ± 5.2 kcal/mol. We also observe that the substituted mutant Arg59 is crucial for glycan‐binding and for the preference of α(2,3)‐linked sialyllactose at the binding pocket of ACG2 galectin. When compared with the wild‐type and single point mutant, the double point mutant exhibits enhanced affinity towards α(2,3)‐linked sialyllactose, which can be effectively used as a model for biological cell marker in cancer therapeutics. Copyright © 2015 John Wiley & Sons, Ltd.     

Origins of the specificity of inhibitor P218 toward wild-type and mutant PfDHFR: a molecular dynamics analysis

Molecular dynamics simulations were performed to evaluate the origin of the antimalarial effect of the lead compound P218. The simulations of the ligand in the cavities of wild-type, mutant Plasmodium falciparum Dihydrofolate Reductase (PfDHFR) and the human DHFR revealed the differences in the atomic-level interactions and also provided explanation for the specificity of this ligand toward PfDHFR. The binding free energy estimation using Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that P218 has higher binding affinity (~ ?30 to ?35 kcal/mol) toward PfDHFR (both in wild-type and mutant forms) than human DHFR (~ ?22 kcal/mol), corroborating the experimental observations. Intermolecular hydrogen bonding analysis of the trajectories showed that P218 formed two stable hydrogen bonds with human DHFR (Ile7 and Glu30), wild-type and double-mutant PfDHFR’s (Asp54 and Arg122), while it formed three stable hydrogen bonds with quadruple-mutant PfDHFR (Asp54, Arg59, and Arg122). Additionally, P218 binding in PfDHFR is stabilized by hydrogen bonds with residues Ile14 and Ile164. It was found that mutant residues do not reduce the binding affinity of P218 to PfDHFR, in contrast, Cys59Arg mutation strongly favors inhibitor binding to quadruple-mutant PfDHFR. The atomistic-level details explored in this work will be highly useful for the design of non-resistant novel PfDHFR inhibitors as antimalarial agents.     

Interaction of chitosans and their N-acylated derivatives with lipopolysaccharide of gram-negative bacteria

The interactions of lipopolysaccharide (LPS) with the natural polycation chitosan and its derivatives--high molecular weight chitosans (80 kD) with different degree of acetylation, low molecular weight chitosan (15 kD), acylated oligochitosan (5.5 kD) and chitooligosaccharides (biose, triose, and tetraose)--were studied using ligand-enzyme solid-phase assay. The LPS-binding activity of chitosans (80 kD) decreased with increase in acetylation degree. Affinity of LPS interaction with chitosans increased after introduction of a fatty acid residue at the reducing end of chitosan. Activity of N-monoacylated chitooligosaccharides decreased in the order: oligochitosan --> tetra- > tri- --> disaccharides. The three-dimensional structures of complexes of R-LPS and chitosans with different degree of acetylation, chitooligosaccharides, and their N-monoacylated derivatives were generated by molecular modeling. The number of bonds stabilizing the complexes and the energy of LPS binding with chitosans decreased with increase in acetate group content in chitosans and resulted in changing of binding sites. It was shown that binding sites of chitooligosaccharides on R-LPS overlapped and chitooligosaccharide binding energies increased with increase in number of monosaccharide residues in chitosan molecules. The input of the hydrophobic fragment in complex formation energy is most prominent for complexes in water phase and is due to the hydrophobic interaction of chitooligosaccharide acyl fragment with fatty acid residues of LPS.     

pH modulates the TGF‐β ligands binding to the receptors: a computational analysis

In spite of showing high sequence similarity and forming structurally similar ternary complex in vitro, the in vivo role of TGF‐β1 and TGF‐β3 ligands suggests against their functional redundancy and necessitates the importance for the study of the specificity of these ligands. A comparative computational analysis of binary and ternary complexes of these two ligands shows that anchor residues of ligand and receptor at TGF‐β:TβR2 interface are similar in both complexes. However, the potential anchor residues of TGF‐β at TGF‐β:TβR1 interface are different, Tyr50 and Lys51 in TGF‐β3 complex and Lys60 and Tyr6 in TGF‐β1 complex. Pro55 and Asp57 of TβRI may act as anchor residues in complexes of both ligands along with Ile54 for TGF‐β3 complex and Val61 for TGF‐β1 complex. Arg58 of TβR1 acts as a potential hot residue for TGF‐β3 ternary complex but not for TGF‐β1 ternary complex formation whereas Pro55 and Phe60 may act as hot residues for both complexes. The Delphi analysis of the pH dependence of the binding energy indicates that pH has a remarkable effect on the binding energy of TβR2 to the open form of TGF‐β3. Lowering of pH from 7 to 4 favors binding of the open form of TGF‐β3 to TβR2. Now, apart from the residues at pH 7, residues Arg25, Lys31 and Arg94 of TGF‐β3 and Asp118 and Glu119 of TβR2 also contribute significantly to the binding energy. Contrary to the binding energy of TβR2 to TGF‐β3/TGF‐β1, TβR1 shows appreciable pH dependence for its binding in ternary complex of TGF‐β3/TGF‐β1. In TGF‐β3 ternary complex, the TβR1 electrostatic interaction energy disfavors complex formation at pH 7 while it is favored at pH 4. Copyright © 2014 John Wiley & Sons, Ltd.     

Structure-Based Analysis of the Ligand-Binding Mechanism for DhelOBP21, a C-minus Odorant Binding Protein,from Dastarcus helophoroides (Fairmaire; Coleoptera: Bothrideridae)

Odorant binding proteins (OBPs) transport hydrophobic odor molecules across the sensillar lymph to trigger a neuronal response. Herein, the Minus-C OBP (DhelOBP21) was characterized from Dastarcus helophoroides, the most important natural parasitic enemy insect that targets Monochamus alternatus. Homology modeling and molecular docking were conducted on the interaction between DhelOBP21 and 17 volatile molecules (including volatiles from pine bark, the larva of M. alternatus, and the faeces of the larva). The predicted three-dimensional structure showed only two disulfide bridges and a hydrophobic binding cavity with a short C-terminus. Ligand-binding experiments using N-phenylnaphthylamine (1-NPN) as a fluorescent probe showed that DhelOBP21 exhibited better binding affinities against those ligands with a molecular volume between 100 and 125 ų compared with ligands with a molecular volume between 160 and 185 ų. Molecules that are too big or too small are not conducive for binding. We mutated the amino acid residues of the binding cavity to increase either hydrophobicity or hydrophilia. Ligand-binding experiments and cyber molecular docking assays indicated that hydrophobic interactions are more significant than hydrogen-bonding interactions. Although hydrogen-bond interactions could be predicted for some binding complexes, the hydrophobic interactions had more influence on binding following hydrophobic changes that affected the cavity. The orientation of ligands affects binding by influencing hydrophobic interactions. The binding process is controlled by multiple factors. This study provides a basis to explore the ligand-binding mechanisms of Minus-C OBP.     

Functional role of putative critical residues in Mycobacterium tuberculosis RNase P protein

RNase P is involved in processing the 5⿲ end of pre-tRNA molecules. Bacterial RNase P contains a catalytic RNA subunit and a protein subunit. In this study, we have analyzed the residues in RNase P protein of M. tuberculosis that differ from the residues generally conserved in other bacterial RNase Ps. The residues investigated in the current study include the unique residues, Val27, Ala70, Arg72, Ala77, and Asp124, and also Phe23 and Arg93 which have been found to be important in the function of RNase P protein components of other bacteria. The selected residues were individually mutated either to those present in other bacterial RNase P protein components at respective positions or in some cases to alanine. The wild type and mutant M. tuberculosis RNase P proteins were expressed in E. coli, purified, used to reconstitute holoenzymes with wild type RNA component in vitro, and functionally characterized. The Phe23Ala and Arg93Ala mutants showed very poor catalytic activity when reconstituted with the RNA component. The catalytic activity of holoenzyme with Val27Phe, Ala70Lys, Arg72Leu and Arg72Ala was also significantly reduced, whereas with Ala77Phe and Asp124Ser the activity of holoenzyme was similar to that with the wild type protein. Although the mutants did not suffer from any binding defects, Val27Phe, Ala70Lys, Arg72Ala and Asp124Ser were less tolerant towards higher temperatures as compared to the wild type protein. The K_m of Val27Phe, Ala70Lys, Arg72Ala and Ala77Phe were >2-fold higher than that of the wild type, indicating the substituted residues to be involved in substrate interaction. The study demonstrates that residues Phe23, Val27 and Ala70 are involved in substrate interaction, while Arg72 and Arg93 interact with other residues within the protein to provide it a functional conformation.     

Synthesis,characterization, and binding assessment with human serum albumin of three bipyridine lanthanide(III) complexes

In this work, the terbium(III), dysprosium(III), and ytterbium(III) complexes containing 2, 2′-bipyridine (bpy) ligand have been synthesized and characterized using CHN elemental analysis, FT-IR, UV–Vis and ¹H-NMR techniques and their binding behavior with human serum albumin (HSA) was studied by UV–Vis, fluorescence and molecular docking examinations. The experimental data indicated that all three lanthanide complexes have high binding affinity to HSA with effective quenching of HSA fluorescence via static mechanism. The binding parameters, the type of interaction, the value of resonance energy transfer, and the binding distance between complexes and HSA were estimated from the analysis of fluorescence measurements and Förster theory. The thermodynamic parameters suggested that van der Waals interactions and hydrogen bonds play an important role in the binding mechanism. While, the energy transfer from HSA molecules to all these complexes occurs with high probability, the order of binding constants (BpyTb > BpyDy > BpyYb) represents the importance of radius of Ln³⁺ ion in the complex-HSA interaction. The results of molecular docking calculation and competitive experiments assessed site 3 of HSA, located in subdomain IB, as the most probable binding site for these ligands and also indicated the microenvironment residues around the bound mentioned complexes. The computational results kept in good agreement with experimental data.     

Structure-function analysis of SH3 domains: SH3 binding specificity altered by single amino acid substitutions.

SH3 domains mediate intracellular protein-protein interactions through the recognition of proline-rich sequence motifs on cellular proteins. Structural analysis of the Src SH3 domain (Src SH3) complexed with proline-rich peptide ligands revealed three binding sites involved in this interaction: two hydrophobic interactions (between aliphatic proline dipeptides in the SH3 ligand and highly conserved aromatic residues on the surface of the SH3 domain), and one salt bridge (between Asp-99 of Src and an Arg three residues upstream of the conserved Pro-X-X-Pro motif in the ligand). We examined the importance of the arginine binding site of SH3 domains by comparing the binding properties of wild-type Src SH3 and Abl SH3 with those of a Src SH3 mutant containing a mutated arginine binding site (D99N) and Abl SH3 mutant constructs engineered to contain an arginine binding site (T98D and T98D/F91Y). We found that the D99N mutation diminished binding to most Src SH3-binding proteins in whole cell extracts; however, there was only a moderate reduction in binding to a small subset of Src SH3-binding proteins (including the Src substrate p68). p68 was shown to contain two Arg-containing Asp-99-dependent binding sites and one Asp-99-independent binding site which lacks an Arg. Moreover, substitution of Asp for Thr-98 in Abl SH3 changed the binding specificity of this domain and conferred the ability to recognize Arg-containing ligands. These results indicate that Asp-99 is important for Src SH3 binding specificity and that Asp-99-dependent binding interactions play a dominant role in Src SH3 recognition of cellular binding proteins, and they suggest the existence of two Src SH3 binding mechanisms, one requiring Asp-99 and the other independent of this residue.     

In silico analyses of substrate interactions with human serum paraoxonase 1

Human paraoxonase (HuPON1) is a serum enzyme that exhibits a broad spectrum of hydrolytic activities, including the hydrolysis of various organophosphates, esters, and recently identified lactone substrates. Despite intensive site-directed mutagenesis and other biological studies, the structural basis for the specificity of substrate interactions of HuPON1 remains elusive. In this study, we apply homology modeling, docking, and molecular dynamic (MD) simulations to probe the binding interactions of HuPON1 with representative substrates. The results suggest that the active site of HuPON1 is characterized by two distinct binding regions: the hydrophobic binding site for arylesters/lactones, and the paraoxon binding site for phosphotriesters. The unique binding modes proposed for each type of substrate reveal a number of key residues governing substrate specificity. The polymorphic residue R/Q192 interacts with the leaving group of paraoxon, suggesting it plays an important role in the proper positioning of this substrate in the active site. MD simulations of the optimal binding complexes show that residue Y71 undergoes an "open-closed" conformational change upon ligand binding, and forms strong interactions with substrates. Further binding free energy calculations and residual decomposition give a more refined molecular view of the energetics and origin of HuPON1/substrate interactions. These studies provide a theoretical model of substrate binding and specificity associated with wild type and mutant forms of HuPON1, which can be applied in the rational design of HuPON1 variants as bioscavengers with enhanced catalytic activity.     

Roles of the P1, P2, and P3 residues in determining inhibitory specificity of kallistatin toward human tissue kallikrein

Kallistatin is a serpin with a unique P1 Phe, which confers an excellent inhibitory specificity toward tissue kallikrein. In this study, we investigated the P3-P2-P1 residues (residues 386-388) of human kallistatin in determining inhibitory specificity toward human tissue kallikrein by site-directed mutagenesis and molecular modeling. Human kallistatin mutants with 19 different amino acid substitutions at each P1, P2, or P3 residue were created and purified to compare their kallikrein binding activity. Complex formation assay showed that P1 Arg, P1 Phe (wild type), P1 Lys, P1 Tyr, P1 Met, and P1 Leu display significant binding activity with tissue kallikrein among the P1 variants. Kinetic analysis showed the inhibitory activities of the P1 mutants toward tissue kallikrein in the order of P1 Arg > P1 Phe > P1 Lys >/= P1 Tyr > P1 Leu >/= P1 Met. P1 Phe displays a better selectivity for human tissue kallikrein than P1 Arg, since P1 Arg also inhibits several other serine proteinases. Heparin distinguishes the inhibitory specificity of kallistatin toward kallikrein versus chymotrypsin. For the P2 and P3 variants, the mutants with hydrophobic and bulky amino acids at P2 and basic amino acids at P3 display better binding activity with tissue kallikrein. The inhibitory activities of these mutants toward tissue kallikrein are in the order of P2 Phe (wild type) > P2 Leu > P2 Trp > P2 Met and P3 Arg > P3 Lys (wild type). Molecular modeling of the reactive center loop of kallistatin bound to the reactive crevice of tissue kallikrein indicated that the P2 residue required a long and bulky hydrophobic side chain to reach and fill the hydrophobic S2 cleft generated by Tyr(99) and Trp(219) of tissue kallikrein. Basic amino acids at P3 could stabilize complex formation by forming electrostatic interaction with Asp(98J) and hydrogen bond with Gln(174) of tissue kallikrein. Our results indicate that tissue kallikrein is a specific target proteinase for kallistatin.     

Steered molecular dynamics simulations of ligand–receptor interaction in lipocalins

Retinol binding protein (RBP) and an engineered lipocalin, DigA16, have been studied using molecular dynamics simulations. Special emphasis has been placed on explaining the ligand–receptor interaction in RBP–retinol and DigA16–digoxigenin complexes, and steered molecular dynamics simulations of 10–20 ns have been carried out for the ligand expulsion process. Digoxigenin is bound deep inside the cavity of DigA16 and forms several stable hydrogen bonds in addition to the hydrophobic van der Waals interaction with the aromatic side-chains. Four crystalline water molecules inside the ligand-binding cavity remain trapped during the simulations. The strongly hydrophobic receptor site of RBP differs considerably from DigA16, and the main source of ligand attraction comes from the phenyl side-chains. The hydrogen bonds between digoxigenin and DigA16 cause the rupture forces on ligand removal in DigA16 and RBP to differ. The mutated DigA16 residues contribute approximately one-half of the digoxigenin interaction energy with DigA16 and, of these, the energetically most important are residues His35, Arg58, Ser87, Tyr88, and Phe114. Potential “sensor loops” were found for both receptors. These are the outlier loops between residues 114–121 and 63–67 for DigA16 and RBP, respectively, and they are located near the entrance of the ligand-binding cavity. Especially, the residues Glu119 (DigA16) and Leu64 (RBP) are critical for sensing. The ligand binding energies have been estimated based on the linear response approximation of binding affinity by using a previous parametrization for retinoids and RBP.     

Hierarchy of simulation models in predicting molecular recognition mechanisms from the binding energy landscapes: structural analysis of the peptide complexes with SH2 domains.

Computer simulations using the simplified energy function and simulated tempering dynamics have accurately determined the native structure of the pYVPML, SVLpYTAVQPNE, and SPGEpYVNIEF peptides in the complexes with SH2 domains. Structural and equilibrium aspects of the peptide binding with SH2 domains have been studied by generating temperature-dependent binding free energy landscapes. Once some native peptide-SH2 domain contacts are constrained, the underlying binding free energy profile has the funnel-like shape that leads to a rapid and consistent acquisition of the native structure. The dominant native topology of the peptide-SH2 domain complexes represents an extended peptide conformation with strong specific interactions in the phosphotyrosine pocket and hydrophobic interactions of the peptide residues C-terminal to the pTyr group. The topological features of the peptide-protein interface are primarily determined by the thermodynamically stable phosphotyrosyl group. A diversity of structurally different binding orientations has been observed for the amino-terminal residues to the phosphotyrosine. The dominant native topology for the peptide residues carboxy-terminal to the phosphotyrosine is tolerant to flexibility in this region of the peptide-SH2 domain interface observed in equilibrium simulations. The energy landscape analysis has revealed a broad, entropically favorable topology of the native binding mode for the bound peptides, which is robust to structural perturbations. This could provide an additional positive mechanism underlying tolerance of the SH2 domains to hydrophobic conservative substitutions in the peptide specificity region.     

Docking,molecular dynamics,binding energy-MM-PBSA studies of naphthofuran derivatives to identify potential dual inhibitors against BACE-1 and GSK-3β

BACE-1 and GSK-3β both are potential therapeutic drug targets for Alzheimer’s disease. Recently, both these targets received attention for designing dual inhibitors. Till now only two scaffolds (triazinone and curcumin) derivatives have been reported as BACE-1 and GSK-3β dual inhibitors. In our previous work, we have reported first in class dual inhibitor for BACE-1 and GSK-3β. In this study, we have explored other naphthofuran derivatives for their potential to inhibit BACE-1 and GSK-3β through docking, molecular dynamics, binding energy (MM-PBSA). These computational methods were performed to estimate the binding affinity of naphthofuran derivatives towards the BACE-1 and GSK-3β. In the docking results, two derivatives (NS7 and NS9) showed better binding affinity as compared to previously reported inhibitors. Hydrogen bond occupancy of NS7 and NS9 generated from MD trajectories showed good interaction with the flap residues Gln73, Thr72 of BACE-1 and Arg141, Thr138 residues of GSK-3β. MM-PBSA and energy decomposition per residue revealed different components of binding energy and relative importance of amino acid involved in binding. The results showed that the binding of inhibitors was majorly governed by the hydrophobic interactions and suggesting that hydrophobic interactions might be the key to design dual inhibitors for BACE1-1 and GSK-3β. Distance between important pair of amino acid residues indicated that BACE-1 and GSK-3β adopt closed conformation and become inactive after ligand binding. The results suggested that naphthofuran derivatives might act as dual inhibitor against BACE-1 and GSK-3β.     

胡椒碱分子键合人血清白蛋白位点的表征

利用荧光光谱法、紫外光谱法并结合计算机模拟技术在分子水平上研究了胡椒碱与人血清白蛋白(human serum albumin HSA)的键合作用.同步荧光及紫外光谱图表明,胡椒碱对HSA微环境有影响.位点竞争试验证明,胡椒碱分子键合在HSA的位点Ⅱ区.通过荧光光谱滴定数据求得不同温度下(300K 310K和318 K)药物与蛋白相互作用的结合常数及结合位点数.分子模拟的结果显示了胡椒碱与HSA的键合区域和键合模式,表明药物与蛋白有较强的键合作用;维持药物与蛋白质的相互作用力主要是疏水用,兼有氢键(位于氨基酸残基Arg 257,Arg 222及Arg218位).通过实验数据计算得到的热力学参数(ΔH0与ΔS0的值分别为原33.11 kJ·mol-1和原18.90 J·mol原1·K-1)确定了胡椒碱与HSA分子的相互作用力类型主要为氢键兼范德华力.     

Protein-protein recognition: an experimental and computational study of the R89K mutation in Raf and its effect on Ras binding

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras-Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [gamma-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9+/-1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein-protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 A away, namely residues 37-41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2-3 specific water molecules. This analysis will provide insights into the large experimental database of Ras-Raf mutations.     

Influence of cation-pi interactions on RNA-binding proteins

The energy contribution due to cation-π interactions has been computed for 37 RNA binding proteins. The contribution of these cation-π interacting residues in sequential separation, secondary structure involvement, solvent accessibility, and stabilization centers has been evaluated. Sequential separation of the cation-π involving residues show that, long range contacts predominates in all the proteins studied. Lys and Arg prefers to be in helical structures. Of the cation-π interacting residues, Arg and Lys were in the exposed regions and the aromatic residues (Phe, Tyr and Trp) were in the buried and partially buried regions in the protein structures. Stabilization centers for these proteins showed that all the five residues found in cation-π interactions are important in locating one or more of such centers. On the whole, the results presented in this work will be very useful for further investigations on the specificity and selectivity of RNA binding proteins and also for their structural studies.     

Structure of the central RNA recognition motif of human TIA-1 at 1.95A resolution

T-cell-restricted intracellular antigen-1 (TIA-1) regulates alternative pre-mRNA splicing in the nucleus, and mRNA translation in the cytoplasm, by recognizing uridine-rich sequences of RNAs. As a step towards understanding RNA recognition by this regulatory factor, the X-ray structure of the central RNA recognition motif (RRM2) of human TIA-1 is presented at 1.95 Å resolution. Comparison with structurally homologous RRM-RNA complexes identifies residues at the RNA interfaces that are conserved in TIA-1-RRM2. The versatile capability of RNP motifs to interact with either proteins or RNA is reinforced by symmetry-related protein-protein interactions mediated by the RNP motifs of TIA-1-RRM2. Importantly, the TIA-1-RRM2 structure reveals the locations of mutations responsible for inhibiting nuclear import. In contrast with previous assumptions, the mutated residues are buried within the hydrophobic interior of the domain, where they would be likely to destabilize the RRM fold rather than directly inhibit RNA binding.     

Identification of ligand effector binding sites in transmembrane regions of the human G protein-coupled C3a receptor

The human C3a anaphylatoxin receptor (C3aR) is a G protein-coupled receptor (GPCR) composed of seven transmembrane alpha-helices connected by hydrophilic loops. Previous studies of chimeric C3aR/C5aR and loop deletions in C3aR demonstrated that the large extracellular loop2 plays an important role in noneffector ligand binding; however, the effector binding site for C3a has not been identified. In this study, selected charged residues in the transmembrane regions of C3aR were replaced by Ala using site-directed mutagenesis, and mutant receptors were stably expressed in the RBL-2H3 cell line. Ligand binding studies demonstrated that R161A (helix IV), R340A (helix V), and D417A (helix VII) showed no binding activity, although full expression of these receptors was established by flow cytometric analysis. C3a induced very weak intracellular calcium flux in cells expressing these three mutant receptors. H81A (helix II) and K96A (helix III) showed decreased ligand binding activity. The calcium flux induced by C3a in H81A and K96A cells was also consistently reduced. These findings suggest that the charged transmembrane residues Arg161, Arg340, and Asp417 in C3aR are essential for ligand effector binding and/or signal coupling, and that residues His81 and Lys96 may contribute less directly to the overall free energy of ligand binding. These transmembrane residues in C3aR identify specific molecular contacts for ligand interactions that account for C3a-induced receptor activation.     

Human formyl peptide receptor function role of conserved and nonconserved charged residues.

In this study, we investigated the role of charged residues in ligand binding interactions of f-Met-Leu-Phe receptors (FPR). Charged residues of FPR, both conserved and nonconserved, which are located close to the membrane interface were mutated to alanine to determine their role in ligand binding. The mutated residues belonged to specific domains of FPR which have previously been implicated in FPR ligand binding interactions. We demonstrate that nonconserved charged residues such as Arg84, Lys85, Arg205 and Asp284 and conserved charge residue Arg163 seem to play a role in ligand binding. However, alteration of nonconserved charged residue Asp106 did not have any effect. In conclusion, specific charged residues of FPR, both conserved nonconserved, may contribute to FPR function either directly or indirectly.     

Energetics of hydrogen bond switch, residue burial and cavity analysis reveals molecular basis of improved heparin binding to antithrombin

Antithrombin III (ATIII) is the main inhibitor of the coagulation proteases like factor Xa and thrombin. Anticoagulant activity of ATIII is increased by several thousand folds when activated by vascular wall heparan sulfate proteoglycans (HSPGs) and pharmaceutical heparins. ATIII isoforms in human plasma, alpha-ATIII and beta-ATIII differ in the amount of glycosylation which is the basis of differences in their heparin binding affinity and function. Crystal structures and site directed mutagenesis studies have mapped the heparin binding site in ATIII, however the hydrogen bond switch and energetics of interaction during the course of heparin dependent conformational change remains largely unclear. An analysis of heparin bound conformational states of ATIII using PEARLS software showed that in heparin bound intermediate state, Arg 47 and Arg 13 residues make hydrogen bonds with heparin but in the activated conformation Lys 11 and Lys 114 have more hydrogen bond interactions. In the protease bound-antithrombin-pentasaccharide complex Lys 114, Pro 12 and Lys 125 form important hydrogen bonding interactions. The results showed that A-helix and N-terminal end residues are more important in the initial interactions but D-helix is more important during the latter stage of conformational activation and during the process of protease inhibition. We carried out the residue wise Accessible Surface Area (ASA) analysis of alpha and beta ATIII native states and the results indicated major differences in burial of residues from Ser 112 to Ser 116 towards the N-terminal end. This region is involved in the P-helix formation on account of heparin binding. A cavity analysis showed a progressively larger cavity formation during activation in the region just adjacent to the heparin binding site towards the C-terminal end. We hypothesize that during the process of conformational change after heparin binding beta form of antithrombin has low energy barrier to form D-helix extension toward N and C-terminal end as compared to alpha isoform.