Insights on pH-dependent conformational changes of mosquito odorant binding proteins by molecular dynamics simulations

Chemical recognition plays an important role for the survival and reproduction of many insect species. Odorant binding proteins (OBPs) are the primary components of the insect olfactory mechanism and have been documented to play an important role in the host-seeking mechanism of mosquitoes. They are “transport proteins” believed to transport odorant molecules from the external environment to their respective membrane targets, the olfactory receptors. The mechanism by which this transport occurs in mosquitoes remains a conundrum in this field. Nevertheless, OBPs have proved to be amenable to conformational changes mediated by a pH change in other insect species. In this paper, the effect of pH on the conformational flexibility of mosquito OBPs is assessed computationally using molecular dynamics simulations of a mosquito OBP “CquiOBP1” bound to its pheromone 3OG (PDB ID: 3OGN). Conformational twist of a loop, driven by a set of well-characterized changes in intramolecular interactions of the loop, is demonstrated. The concomitant (i) closure of what is believed to be the entrance of the binding pocket, (ii) expansion of what could be an exit site, and (iii) migration of the ligand towards this putative exit site provide preliminary insights into the mechanism of ligand binding and release of these proteins in mosquitoes. The correlation of our results with previous experimental observations based on NMR studies help us provide a cardinal illustration on one of the probable dynamics and mechanism by which certain mosquito OBPs could deliver their ligand to their membrane-bound receptors at specific pH conditions.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

Thermodynamic dissection of the binding energetics of proline-rich peptides to the Abl-SH3 domain: implications for rational ligand design

The inhibition of the interactions between SH3 domains and their targets is emerging as a promising therapeutic strategy. To date, rational design of potent ligands for these domains has been hindered by the lack of understanding of the origins of the binding energy. We present here a complete thermodynamic analysis of the binding energetics of the p41 proline-rich decapeptide (APSYSPPPPP) to the SH3 domain of the c-Abl oncogene. Isothermal titration calorimetry experiments have revealed a thermodynamic signature for this interaction (very favourable enthalpic contributions opposed by an unfavourable binding entropy) inconsistent with the highly hydrophobic nature of the p41 ligand and the Abl-SH3 binding site. Our structural and thermodynamic analyses have led us to the conclusion, having once ruled out any possible ionization events or conformational changes coupled to the association, that the establishment of a complex hydrogen-bond network mediated by water molecules buried at the binding interface is responsible for the observed thermodynamic behaviour. The origin of the binding energetics for proline-rich ligands to the Abl-SH3 domain is further investigated by a comparative calorimetric analysis of a set of p41-related ligands. The striking effects upon the enthalpic and entropic contributions provoked by conservative substitutions at solvent-exposed positions in the ligand confirm the complexity of the interaction. The implications of these results for rational ligand design are discussed.     

Quantitative NMR studies of high molecular weight proteins: application to domain orientation and ligand binding in the 723 residue enzyme malate synthase G

A high-resolution multidimensional NMR study of ligand-binding to Escherichia coli malate synthase G (MSG), a 723-residue monomeric enzyme (81.4 kDa), is presented. MSG catalyzes the condensation of glyoxylate with an acetyl group of acetyl-CoA, producing malate, an intermediate in the citric-acid cycle. We show that despite the size of the protein, important structural and dynamic information about the molecule can be obtained on a per-residue basis. 15N-1HN residual dipolar couplings and carbonyl chemical shift changes upon alignment in Pf1 phage establish that there are no significant domain reorientations in the molecule upon ligand binding, in contrast to what was anticipated on the basis of both the X-ray structure of the glyoxylate-bound form of the enzyme and structural studies of a related set of proteins. The chemical shift changes of 1HN, 15N and 13CO nuclei upon binding of pyruvate, a glyoxylate-mimicking inhibitor, and acetyl-CoA have been mapped onto the three-dimensional structure of the molecule. Binding constants of pyruvate, glyoxylate, and acetyl-CoA (in the presence of pyruvate) have been measured, along with the kinetic parameters for glyoxylate and pyruvate binding. The on-rates of pyruvate and glyoxalate binding, approximately 1.2 x 10(6)M(-1)s(-1) and approximately 2.7 x 10(6)M(-1)s(-1), respectively, are significantly lower than what is anticipated from a simple diffusion-controlled process. Some structural implications of the chemical shift perturbations upon binding and the estimated ligand on-rates are discussed.     

Deciphering adaptor specificity in GFL-dependent RET-mediated proliferation and neurite outgrowth

Glial cell derived neurotrophic factor (GDNF)-dependent receptor tyrosine kinase RET activity is required for proper development of the nervous system and genitourinary tract. Loss-of-function mutations in RET are associated with enteric nervous system abnormalities (Hirschsprung disease) and renal deficits (Potter's syndrome), whereas activating mutations lead to hereditary cancer syndromes (multiple endocrine neoplasia type 2A and type 2B). RET activation is crucial for the proper regulation of a variety of cellular processes including cell migration, proliferation and neurite outgrowth. By analyzing a series of RET mutants we found that Y1062 was critical for stimulating GDNF-mediated proliferation as well as proliferation stimulated by GDNF-independent oncogenic forms of RET. Studies using small interfering RNA driven by lentivirus to knock-down expression of particular adaptor proteins that interact with RET phospho-Y1062, demonstrated that only Src-homology 2 and growth factor receptor binding protein 2 were necessary for RET-mediated proliferation by wild type and oncogenic forms of RET. Interestingly, we discovered that Y1062 was also required for GDNF-stimulated neurite outgrowth. However, small interfering RNAs to either Src-homology 2 or growth factor receptor binding protein 2 or a panel of other adaptor proteins known to interact with RET Y1062 were incapable of blocking GDNF-stimulated neurite formation, indicating that differential use of intracellular adaptors is responsible for regulating alternative RET-stimulated cellular events such as proliferation versus a differentiation response like neurite outgrowth.     

Interactions in native binding sites cause a large change in protein dynamics

Cellular functions are regulated by molecules that interact with proteins and alter their activities. To enable such control, protein activity, and therefore protein conformational distributions, must be susceptible to alteration by molecular interactions at functional sites. Here we investigate whether interactions at functional sites cause a large change in the protein conformational distribution. We apply a computational method, called dynamics perturbation analysis (DPA), to identify sites at which interactions have a large allosteric potential D(x), which is the Kullback-Leibler divergence between protein conformational distributions with and without an interaction. In DPA, a protein is decorated with surface points that interact with neighboring protein atoms, and D(x) is calculated for each of the points in a coarse-grained model of protein vibrations. We use DPA to examine hundreds of protein structures from a standard small-molecule docking test set, and find that ligand-binding sites have elevated values of D(x): for 95% of proteins, the probability of randomly obtaining values as high as those in the binding site is 10(-3) or smaller. We then use DPA to develop a computational method to predict functional sites in proteins, and find that the method accurately predicts ligand-binding-site residues for proteins in the test set. The performance of this method compares favorably with that of a cleft analysis method. The results confirm that interactions at small-molecule binding sites cause a large change in the protein conformational distribution, and motivate using DPA for large-scale prediction of functional sites in proteins. They also suggest that natural selection favors proteins whose activities are capable of being regulated by molecular interactions.     

C2H2型锌指蛋白的研究进展

锌指基因家族是人体中最大的基因家族,它参与细胞分化、胚胎发育,并与许多疾病的发生相关．根据半胱氨酸(c)和组氨酸(H)的数目和位置可将锌指蛋白分为c2H2、c2Hc2、c2c2 CHCC2C2、C2C2C2C2等亚类．c2H2型锌指是最普遍的类型,它们作为重要的转录调控因子参与许多的生理过程．c2H2型锌指蛋白包含的锌指数目从1个到30多个不等．依据锌指的数量以及在蛋白中的分布情况,大多数c2H2型锌指蛋白属于下列3类之一：1)含3个c：H：锌指的蛋白(tC2H2);2)含多个锌指的c2H2型锌指结构蛋白(mac2H2);3)锌指成对间隔排列的c2H2型锌指蛋白(spC2H2)、一些c2H2型锌指蛋白能识别并结合特异性RNA或DNA片段．另一些则只能与RNA结合．通常锌指蛋白含锌指数目越多。它选择结合的能力就越强．     

Molecular dynamics characterization of the active cavity of carboxypeptidase A and some of its inhibitor adducts.

Molecular dynamics (MD) calculations have been performed on carboxypeptidase A and on its adducts with inhibitors, such as d-phenylalanine (dPhe) and acetate. The catalytically essential zinc ion present in the protein was explicitly included in all the simulations. The simulation was carried out over a sphere of 15 A centered on the zinc ion. The crystallographic water molecules were explicitly taken into account; then the protein was solvated with a 18 A sphere of water molecules. MD calculations were carried out for 45-60 ps. There is no large deviation from the available X-ray structures of native and the dPhe adduct for the MD structures. Average MD structures were calculated starting from the X-ray structure of the dPhe adduct, and, from a structure obtained by docking the inhibitor in the native structure. Comparison between these two structures and with that of the native protein shows that some of the key variations produced by inhibitor binding are reproduced by MD calculations. Addition of acetate induces structural changes relevant for the understanding of the interaction network in the active cavity. The structural variations induced by different inhibitors are examined. The effects of these interactions on the catalytic mechanism and on the binding of substrate are discussed.     

Structural activation pathways from dynamic olfactory receptor-odorant interactions

We have simulated an odor ligand's dynamic behavior in the binding region of an olfactory receptor (OR). Our short timescale computational studies (up to 200 ps) have helped identify unprecedented postdocking ligand behavior of ligands. From in vacuo molecular dynamics simulations of interactions between models of rat OR I7 and 10 aldehyde ligands, we have identified a dissociative pathway along which the ligand exits and enters the OR-binding pocket--a transit event. The ligand's transit through the receptor's binding region may mark the beginning of a signal transduction cascade leading to odor recognition. We have graphically traced the rotameric changes in key OR amino acid side chains during the transit. Our results have helped substantiate or refute previously held notions of amino acid contribution to ligand stability in the binding pocket. Our observations of ligand activity when compared to those of experimental (electroolfactogram response) OR-activation studies provide a view to predicting the stability of ligands in the binding pocket as a precursor to OR activation by the ligand.