Minimal models of multi-site ligand-binding kinetics

Mathematical models based on the current understanding of co-operativity in ligand binding to the (macro) molecule and relating the dose-response (saturation) curve of the (macro) molecule ligation to intrinsic dissociation constants characterizing the affinities of ligand for binding sites of both unliganded and partly liganded (macro) molecule have been developed. The simplified models disregarding the structural properties and considerations concerning conformational changes of the (macro) molecule retain the ability to yield sigmoid curves of ligand binding and reflect the co-operativity. Model 1 contains only three parameters, parameter κ (a multiplier characterising the change in the affinity) reflects also the existence and type of co-operativity of ligand binding: κ<1 corresponds to positive co-operativity, κ>1 to the negative and κ=1 to the absence of any co-operativity. Model 2 contains an extra parameter, ω, equilibrium constant for the T₀↔R₀ transition but fails to produce dose-response, which would suggest negative co-operativity. For any fixed n>1, the deviation of the dose-response (saturation) curve from the Henri hyperbola depends either solely on parameter κ (Model 1) or also on parameter ω (Model 2). The (macro) molecule being a receptor, both models yield a diversity of dose-response curves due to possible variety of efficacies of the (macro) molecule. The models may be considered as extensions of the Henri model: in case the dissociation constants remain unchanged, the proposed models are reduced to the latter.     

Improved PCR method for the creation of saturation mutagenesis libraries in directed evolution: application to difficult-to-amplify templates

Saturation mutagenesis constitutes a powerful method in the directed evolution of enzymes. Traditional protocols of whole plasmid amplification such as Stratagene’s QuikChange™ sometimes fail when the templates are difficult to amplify. In order to overcome such restrictions, we have devised a simple two-primer, two-stage polymerase chain reaction (PCR) method which constitutes an improvement over existing protocols. In the first stage of the PCR, both the mutagenic primer and the antiprimer that are not complementary anneal to the template. In the second stage, the amplified sequence is used as a megaprimer. Sites composed of one or more residues can be randomized in a single PCR reaction, irrespective of their location in the gene sequence.The method has been applied to several enzymes successfully, including P450-BM3 from Bacillus megaterium, the lipases from Pseudomonas aeruginosa and Candida antarctica and the epoxide hydrolase from Aspergillus niger. Here, we show that megaprimer size as well as the direction and design of the antiprimer are determining factors in the amplification of the plasmid. Comparison of the results with the performances of previous protocols reveals the efficiency of the improved method. Joaquin Sanchis, Layla Fernández, and J. Daniel Carballeira contributed equally.     

Femtomolar Fab binding affinities to a protein target by alternative CDR residue co-optimization strategies without phage or cell surface display

In therapeutic or diagnostic antibody discovery, affinity maturation is frequently required to optimize binding properties. In some cases, achieving very high affinity is challenging using the display-based optimization technologies. Here we present an approach that begins with the creation and clonal, quantitative analysis of soluble Fab libraries with complete diversification in adjacent residue pairs encompassing every complementarity-determining region position. This was followed by alternative recombination approaches and high throughput screening to co-optimize large sets of the found improving mutations. We applied this approach to the affinity maturation of the anti-tumor necrosis factor antibody adalimumab and achieved ~500-fold affinity improvement, resulting in femtomolar binding. To our knowledge, this is the first report of the in vitro engineering of a femtomolar affinity antibody against a protein target without display screening. We compare our findings to a previous report that employed extensive mutagenesis and recombination libraries with yeast display screening. The present approach is widely applicable to the most challenging of affinity maturation efforts.     

Synthesis and Characterization of Time-resolved Fluorescence Probes for Evaluation of Competitive Binding to Melanocortin Receptors

Probes for use in time-resolved fluorescence competitive binding assays at melanocortin receptors based on the parental ligands MSH(4), MSH(7), and NDP-α-MSH were prepared by solid phase synthesis methods, purified, and characterized. The saturation binding of these probes was studied using HEK-293 cells engineered to overexpress the human melanocortin 4 receptor (hMC4R) as well as the human cholecystokinin 2 receptor (hCCK2R). The ratios of non-specific binding to total binding approached unity at high concentrations for each probe. At low probe concentrations, receptor-mediated binding and uptake was discernable, and so probe concentrations were kept as low as possible in determining K_d values. The Eu-DTPA-PEGO-MSH(4) probe exhibited low specific binding relative to non-specific binding, even at low nanomolar concentrations, and was deemed unsuitable for use in competition binding assays. The Eu-DTPA-PEGO probes based on MSH(7) and NDP-α-MSH exhibited K_d values of 27 ± 3.9 nM and 4.2 ± 0.48 nM, respectively, for binding with hMC4R. These probes were employed in competitive binding assays to characterize the interactions of hMC4R with monovalent and divalent MSH(4), MSH(7), and NDP-α-MSH constructs derived from squalene. Results from assays with both probes reflected only statistical enhancements, suggesting improper ligand spacing on the squalene scaffold for the divalent constructs. The K_i values from competitive binding assays that employed the MSH(7)-based probe were generally lower than the K_i values obtained when the probe based on NDP-α-MSH was employed, which is consistent with the greater potency of the latter probe. The probe based on MSH(7) was also competed with monovalent, divalent, and trivalent MSH(4) constructs that previously demonstrated multivalent binding in competitive binding assays against a variant of the probe based on NDP-α-MSH. Results from these assays confirm multivalent binding, but suggest a more modest increase in avidity for these MSH(4) constructs than was previously reported.     

On a population pathogen model incorporating species dispersal with temporal variation in dispersal rate

In the present paper, we consider a mathematical model of ecosystem population interaction where the population suffers from a susceptible–infectious–susceptible disease. Dispersal of both the susceptible and the infective is incorporated using reaction–diffusion equations. We first study the stability criteria of the basic (non-spatial) model around the disease-free and the infected steady states. We find that the loss rate of the infective species controls disease prevalence. Also without predation pressure, the disease will continue to exist among the population. Then we analyze the spatial model with species dispersal in constant as well as in time-varying form. It is observed that though constant dispersal is unable to generate diffusion-driven instability, dispersal with sinusoidal variation in dispersion rate can generate diffusive instability when the wave number of the perturbation lies within a given range. Numerical simulations are performed to illustrate analytical studies.     

Simulation and estimation of gene number in a biological pathway using almost complete saturation mutagenesis screening of haploid mouse cells

Background

Genome-wide saturation mutagenesis and subsequent phenotype-driven screening has been central to a comprehensive understanding of complex biological processes in classical model organisms such as flies, nematodes, and plants. The degree of “saturation” (i.e., the fraction of possible target genes identified) has been shown to be a critical parameter in determining all relevant genes involved in a biological function, without prior knowledge of their products. In mammalian model systems, however, the relatively large scale and labor intensity of experiments have hampered the achievement of actual saturation mutagenesis, especially for recessive traits that require biallelic mutations to manifest detectable phenotypes.

Results

By exploiting the recently established haploid mouse embryonic stem cells (ESCs), we present an implementation of almost complete saturation mutagenesis in a mammalian system. The haploid ESCs were mutagenized with the chemical mutagen N-ethyl-N-nitrosourea (ENU) and processed for the screening of mutants defective in various steps of the glycosylphosphatidylinositol-anchor biosynthetic pathway. The resulting 114 independent mutant clones were characterized by a functional complementation assay, and were shown to be defective in any of 20 genes among all 22 known genes essential for this well-characterized pathway. Ten mutants were further validated by whole-exome sequencing. The predominant generation of single-nucleotide substitutions by ENU resulted in a gene mutation rate proportional to the length of the coding sequence, which facilitated the experimental design of saturation mutagenesis screening with the aid of computational simulation.

Conclusions

Our study enables mammalian saturation mutagenesis to become a realistic proposition. Computational simulation, combined with a pilot mutagenesis experiment, could serve as a tool for the estimation of the number of genes essential for biological processes such as drug target pathways when a positive selection of mutants is available.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1016) contains supplementary material, which is available to authorized users.     

NMR spectroscopic and molecular modeling studies of protein-carbohydrate and protein-peptide interactions

Investigations of the conformations of carbohydrates, their analogues and their molecular mimics are described, with emphasis on structural and functional information that can be gained by NMR spectroscopic techniques in combination with molecular modeling. The transferred nuclear Overhauser effect (trNOE) has been employed to determine the bound conformations of carbohydrates and other bioactive molecules in complex with protein receptors. The corresponding experiments in the rotating frame (trROE) and selective editing experiments (e.g., QUIET-NOESY) are used to eliminate indirect cross-relaxation pathways (spin diffusion), thereby minimizing errors in the data used for calculation of conformations. Saturation transfer difference NMR experiments reveal detailed information about intermolecular contacts between ligand and protein. Computational techniques are integrated with NMR-derived information to construct structural models of these bioactive molecules and of their complexes with proteins. Recent investigations into the nature of molecular mimicry with regard to protein-ligand interactions are described, along with applications in determining the mode of action of enzyme inhibitors. The results are relevant for the design of the next generation of drug and vaccine candidates.     

Saturation labeling with cysteine-reactive cyanine fluorescent dyes provides increased sensitivity for protein expression profiling of laser-microdissected clinical specimens

Laser capture microdissection (LCM) provides the capability to isolate and analyze small numbers of cells from a specific area of a histologic section. LCM has particular value for analysis of early stage tumors, which are often small and intermixed with non-tumor tissue. It has previously been shown that a new generation of cysteine-reactive cyanine dyes can, in principle, provide increased sensitivity for two-dimensional fluorescence difference gel electrophoresis (2-D DIGE) profiling when sample quantitities are limiting. However, the comparative advantage of the new dyes in a clinical setting has not been established. Here, we report that cysteine-reactive dyes allowed the identification of more features than established, lysine-reactive dyes with a given number of cells. This was true both with extracts prepared from human papillomavirus E6 and E7-transduced human keratinocytes, a model for early-stage cervical cancer, and with LCM samples. In an experiment comparing LCM clinical samples of gastric adenocarcinoma versus precancerous, spasmolytic polypeptide expressing metaplasia (SPEM) from the same patient, cysteine labeling allowed the identification of more than 1000 discrete protein spots in samples containing 5000 cells. This is a 5- to 50-fold smaller sample than used in previous studies. Both labeling methods had a comparable success rate for protein identification by mass spectrometry (MS). The proteins associated with more than 40 differentially abundant spots in the clinical samples were identified by MS. In this exploratory analysis, changes in expression levels of cytoskeletal proteins, molecular chaperones, and cell-signaling proteins were seen. The identification of a number of proteins that are potentially relevant to tumor progression suggests that the method holds promise for biomarker discovery.     

Water quality and soil properties of ponds as influenced by Makhana (Euryale ferox) cultivation

Summary Analysis of regularly collected samples of soil and water from Makhana (Euryale ferox) growing ponds indicated that pH of soil and water were established near neutrality and C/N of soil near ten. Bicarbonate dominated soil extract while chloride dominated standing water. During grand growth period (March to June), Ca²⁺+Mg²⁺, NH ₄ ⁺ K⁺ increased while Na⁺ increased during decay of the plant (July to September) both in soil and plants. The soil had high contents of micro-nutrients, K and P. S.A.R. of soil was 3 to 4 times that of water. Free energy of exchange indicated stability and spontaneity.     

A 31P-NMR saturation transfer study of the regulation of creatine kinase in the rat heart

(1) ³¹P nuclear magnetic resonance was used to measure the creatine kinase-catalysed fluxes in Langendorff-perfused rat hearts consuming oxygen at different rates and using either of two exogenous substrates (11 mM glucose or 5 mM acetate). (2) Fluxes in the direction of ATP synthesis were between 3.5–12-times the steady-state rates of ATP utilization (estimated from rates of O₂-consumption), demonstrating that the reaction is sufficiently rapid to maintain the cytosolic reactants near their equilibrium concentrations. (3) Under all conditions studied, the cytosolic free [ADP] was primarily responsible for regulating the creatine kinase fluxes. The enzyme displayed a K_m for cytosolic ADP of 35 μM and an apparent V_max of 5.5 mM/s in the intact tissue. (4) Although the reaction is maintained in an overall steady-state, the measured ratio of the forward flux (ATP synthesis) to the reverse flux (phosphocreatine synthesis) was significantly greater than unity under some conditions. It is proposed that this discrepancy may be a consequence of participation of ATP in reactions other than the PCr /ag ATP or ATP /ag ADP + P_i interconversions specifically considered in the analysis. (5) The results support the view that creatine kinase functions primarily to maintain low cytosolic concentrations of ADP during transient periods in which energy utilization exceeds production.