Repertoire selection of variant single-chain Cro: toward directed DNA-binding specificity of helix-turn-helix proteins

A single-chain derivative of the lambda Cro repressor (scCro) has been randomly mutated in amino acid residues critical for specific DNA recognition to create libraries of protein variants. Utilizing phage display-afforded affinity selection, scCro variants have been isolated for binding to synthetic DNA ligands. Isolated scCro variants were analyzed functionally, both in fusion with phage particles and after expression of the corresponding free proteins. The binding properties with regard to specificity and affinity in binding to different DNA ligands were investigated by inhibition studies and determination of equilibrium dissociation constants for formed complexes. Variant proteins with altered DNA-sequence specificity were identified, which favored binding of targeted synthetic DNA sequences over a consensus operator sequence, bound with high affinity by wild-type Cro. The specificities were relatively modest (2-3-fold, as calculated from K(D) values), which can be attributed to the inherent properties in the design of the selection system; one half-site of the synthetic DNA sequences maintains the consensus operator sequence, and one "subunit" of the variant single-chain Cro dimers was conserved as wild-type sequence. The anticipated interaction between the wild-type subunit and the consensus DNA half-site of target DNA ligands is, hence, expected to contribute to the overlap in sequence discrimination. The binding affinity for the synthetic DNA sequences, however, was improved 10-30-fold in selected variant proteins as compared to "wild-type" scCro.     

Identification of a polymeric β-cyclodextrin-binding peptide from a phage-displayed peptide library

In this work, a phage-displayed peptide library was applied to identification of β-cyclodextrin (CD)-binding peptide tag, capable of being combined to target peptides or proteins in a homogeneous way by established methods such as peptide synthesis and the recombinant DNA technique. Four enriched sequences were obtained after five rounds of biopanning against polymeric β-CD beads. One of the sequences showed high binding affinity to β-CD beads with a dissociation constant of approximately 7 × 10^–6 M. The β-CD-binding sequence was used for immobilization of a hepatitis C virus (HCV) antigenic peptide on β-CD beads. The functionalized β-CD beads were successfully used for immunoassay of anti-HCV antibody with a detection limit of 1 ng. These results demonstrate that the identified peptide sequence has the potential of being used as an affinity tag to β-CD-containing surfaces.     

A size dependent predator-prey interaction: who pursues whom?

We investigate the properties of an (age, size) -structured model for a population of Daphnia that feeds on a dynamical algal food source. The stability of the internal equilibrium is studied in detail and combined with numerical studies on the dynamics of the model to obtain insight in the relation between individual behaviour and population dynamical phenomena. Particularly the change in the (age, size)-relation with a change in the food availability seems to be an important behavioural mechanism that strongly influences the dynamics. This influence is partly stabilizing and partly destabilizing and leads to the coexistence of a stable equilibrium and a stable limit cycle or even coexistence of two stable limit cycles for the same parameter values. The oscillations in this case are characterized by drastic changes in the size-structure of the population during a cycle. In addition the model exhibits the usual predator-prey oscillations that characterize Lotka-Volterra models.     

Thermodynamic basis of selectivity in guide-target-mismatched RNA interference

Silencing in RNAi is strongly affected by guide‐strand/target‐mRNA mismatches. Target nucleation is thought to occur at positions 2–8 of the guide (“seed region”); successful hybridization in this region is the primary determinant of target‐binding affinity and hence target cleavage. To define a molecular basis for the target sequence selectivity in RNAi, we studied all possible distinct single mismatches in seven positions of the seed region—a total of 21 substitutions. We report results from soft‐core thermodynamic integration simulations to determine changes in targeting binding‐free energies to Argonaute due to single mismatches in the guide strand, which arise during binding of an imperfectly matched target mRNA. In agreement with experiment, most mismatches impair target binding, consistent with a prominent role for binding affinity changes in RNAi sequence selectivity. Individual Argonaute residues located near the mismatched base pair are found to contribute significantly to binding affinity changes. We also use this methodology to analyze the mismatch‐dependent free energy changes for dissociation of a DNA?RNA hybrid from Argonaute, as a model for the escape of miRNAs from the silencing pathway. Several mismatched sequences of the miRNA have increased affinity to Argonaute, implying that some mismatches may reduce the probability for escape. Furthermore, calculations of base‐substitution‐dependent free energy changes for binding ssDNA reveal mild sequence sensitivity as expected for guide strand binding to Argonaute. Our findings give a thermodynamic basis for RNAi target sequence selectivity and suggest that miRNA mismatches may increase silencing effectiveness and thus could be evolutionarily advantageous. Proteins 2012; © 2011 Wiley Periodicals, Inc.     

Hybridization of Antisense Oligonucleotides with α‐Sarcin Loop Region of Escherichia coli 23S rRNA

Binding of complementary oligonucleotides (ONs) with α‐sarcin loop region (2638–2682) of Escherichia coli 23S rRNA was investigated. Four of the tested pentadecanucleotides efficiently bound to target sequences with association rate and equilibrium constants ～ 10³ M^? 1s^? 1 and 10⁷ M^? 1, respectively. ON S5 (CGAGAGGACCGGAGU) complementary to the sequence 2658–2672 displayed the highest affinity to the target. Activation energy for binding of ON S5 was measured to be 11 kcal/mol; this value corresponds to ～ 10% of the calculated enthalpy of the local RNA structure unfolding in the presence of this oligonucleotide. The activation energy value is evidence for the heteroduplex formation to occur via strand displacement pathway; the initiation of heteroduplex formation requires disruption of 1–2 base pairs in RNA hairpin.     

Population extinction and quasi-stationary behavior in stochastic density-dependent structured models

Density-independent and density-dependent, stochastic and deterministic, discrete-time, structured models are formulated, analysed and numerically simulated. A special case of the deterministic, density-independent, structured model is the well-known Leslie age-structured model. The stochastic, density-independent model is a multitype branching process. A review of linear, density-independent models is given first, then nonlinear, density-dependent models are discussed. In the linear, density-independent structured models, transitions between states are independent of time and state. Population extinction is determined by the dominant eigenvalue λ of the transition matrix. If λ ≤ 1, then extinction occurs with probability one in the stochastic and deterministic models. However, if λ > 1, then the deterministic model has exponential growth, but in the stochastic model there is a positive probability of extinction which depends on the fixed point of the system of probability generating functions. The linear, density-independent, stochastic model is generalized to a nonlinear, density-dependent one. The dependence on state is in terms of a weighted total population size. It is shown for small initial population sizes that the density-dependent, stochastic model can be approximated by the density-independent, stochastic model and thus, the extinction behavior exhibited by the linear model occurs in the nonlinear model. In the deterministic models there is a unique stable equilibrium. Given the population does not go extinct, it is shown that the stochastic model has a quasi-stationary distribution with mean close to the stable equilibrium, provided the population size is sufficiently large. For small values of the population size, complete extinction can be observed in the simulations. However, the persistence time increases rapidly with the population size. This author received partial support by the National Science Foundation grant # DMS-9626417.     

High affinity DNA-microtubule interactions: evidence for a conserved DNA-MAP interaction involving unusual high CsCl density repetitious DNA families

We have examined high affinity interactions of chick brain microtubule proteins with ³⁵S labelled tracer DNAs from chick, mouse and D. melanogaster under equilibrium conditions by the nitrocellulose filter binding technique. Ternary reaction mixtures of the above two components and a third component, an excess of unlabelled competitor DNA from either E. coli., mouse, D. melanogaster or chick, were used to measure small fractions of DNA in each case (1–4%) bound to microtubule protein under high stringency- large competitor DNA concentration and 0.5 M NaCl. As seen in part previously (Marx, K.A. and Denial, T. (1985) in The Molecular Basis of Cancer, 172B, 65–75 (Rein, ed), A. Liss, N.Y.) the measured order of competitor DNA strengths was identical for all three tracer DNAs. That is: chick > mouse > D. melanogaster > E. coli competitor DNA. Since the homologous interaction, chick competitor DNA with chick brain microtubule protein, is always the strongest interaction measured, we interpret this as evidence for a conserved protein-DNA sequence interaction. ³⁵S chick DNA tracer sequences, isolated from nitrocellulose filters following the stringent binding in the presence of 0.9 mM^–1 E. coli. competitor DNA, was used in driven reassociation reactions with total chick driver DNA. This fraction was found to be significantly enriched in repetitive chick DNA sequences. Since we have observed a similar phenomenon in mouse, we then compared the stringent binding mouse sequences and showed that the bulk of these sequences did not cross-hybridize with total chick DNA. Finally, all three ³⁵S tracer DNAs binding to nitrocellulose were isolated and sedimented to equilibrium on CsCl density gradients. The CsCl density distributions from all three DNAs showed significant (100-fold) enrichment in classical satellite DNAs as well as higher enrichment in two very unusual high CsCl density families of DNA (1.720–1.740 g/cm³; 1.750–1.765 g/cm³). These families are never observed as distinct bands in total DNA CsCl gradients, nor could we isolate them in purified tubulin control binding experiments. This apparently general phenomena may be identifying some of the sequence families involved in the high affinity microtubule interaction, which appears to be conserved in evolution.     

Pair wise binding affinity: 3D QSAR studies on a set of triazolo [1, 5-a] quinoxalines as antagonists of AMPA and KA receptors

The glutamate receptor system is implicated in the development and maintenance of epileptic seizures and it has been reported that compounds showing high affinity for both AMPA and KA binding sites are more potent anticonvulsants than compounds having selective affinity toward AMPA or KA receptor. These outcomes make such inhibitors future potential antiepileptic drugs. So, the pair wise binding affinity for AMPA and KA receptors inhibition was proposed by using the addition between biological activities of ligands. This approach for evaluation of pair wise binding affinity was exemplified using set of triazolo [1,5-a] quinoxaline for AMPA and KA receptors. The biological activity towards AMPA and KA receptors (expressed as -log IC_5O) was taken as a dependent variable for building CoMFA and CoMSIA models. The resulting models show the ways of increasing binding affinity to both AMPA and KA receptors as potential target for epilepsy. The statistically significant results show that pair wise CoMFA and CoMSIA models are better then individual models. The resulting cross-validated r²_CV value 0.806 for CoMFA is greater then 0.780 for CoMSIA pair wise model. The non-cross validated run giving a coefficient of determination r² value of 0.946 and 0.908 for CoMFA and CoMSIA respectively, provided a good correlation between the observed and computed affinities of the compounds.     

Multi‐constraint computational design suggests that native sequences of germline antibody H3 loops are nearly optimal for conformational flexibility

The limited size of the germline antibody repertoire has to recognize a far larger number of potential antigens. The ability of a single antibody to bind multiple ligands due to conformational flexibility in the antigen‐binding site can significantly enlarge the repertoire. Among the six complementarity determining regions (CDRs) that generally comprise the binding site, the CDR H3 loop is particularly variable. Computational protein design studies showed that predicted low energy sequences compatible with a given backbone structure often have considerable similarity to the corresponding native sequences of naturally occurring proteins, indicating that native protein sequences are close to optimal for their structures. Here, we take a step forward to determine whether conformational flexibility, believed to play a key functional role in germline antibodies, is also central in shaping their native sequence. In particular, we use a multi‐constraint computational design strategy, along with the Rosetta scoring function, to propose that the native sequences of CDR H3 loops from germline antibodies are nearly optimal for conformational flexibility. Moreover, we find that antibody maturation may lead to sequences with a higher degree of optimization for a single conformation, while disfavoring sequences that are intrinsically flexible. In addition, this computational strategy allows us to predict mutations in the CDR H3 loop to stabilize the antigen‐bound conformation, a computational mimic of affinity maturation, that may increase antigen binding affinity by preorganizing the antigen binding loop. In vivo affinity maturation data are consistent with our predictions. The method described here can be useful to design antibodies with higher selectivity and affinity by reducing conformational diversity. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

The spread of a parasitic infection in a spatially distributed host population

Starting from a stochastic model for the spread of a parasitic infection in a spatially distributed host population we describe the way to a continuum formulation by a deterministic model in terms of a nonlinear partial differential equation and an integro-differential equation. The hosts are assumed to occupy fixed spatial positions, whereas the parasites are mobile, however can propagate only within the hosts. To perform the continuum limit we suppose that the size N _h of the host population, the size N _p of the parasite population, and the ratio N _p /N _h tend to infinity. Accordingly, the parameters determining the time evolution of the host and parasite populations are rescaled suitably.Parts of this work have been elaborated during a stay at the Institute of Applied Mathematics at the University of Zürich     

Analysis and Application of an Equilibrium Model forin vitroBioassay Systems with Three Components: Receptor, Hormone and Hormone-Binding-Protein

A simple theoretical framework is presented for bioassay studies using three componentin vitrosystems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an 11-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein toin vitrosystems. The algorithm is tested by application to a published data set from an experimental study in anin vitrosystem (Limet al., 1990,Endocrinology127,1287–1291). Predicted changes show good agreement (within 8%) with experimental observations.     

Subfamily specific conservation profiles for proteins based on n-gram patterns

Background  

A new algorithm has been developed for generating conservation profiles that reflect the evolutionary history of the subfamily associated with a query sequence. It is based on n-gram patterns (NP{n,m}) which are sets of n residues and m wildcards in windows of size n+m. The generation of conservation profiles is treated as a signal-to-noise problem where the signal is the count of n-gram patterns in target sequences that are similar to the query sequence and the noise is the count over all target sequences. The signal is differentiated from the noise by applying singular value decomposition to sets of target sequences rank ordered by similarity with respect to the query.     

The number of segregating sites in a sample of DNA sequences from a geographically structured population

 The distribution of the number of segregating sites among randomly sampled DNA sequences from a geographically structured population is studied. We assume the infinitely-many-sites model of neutral genes and no recombination. Employing the genealogical process, we derive an equation for the generating function of the distribution of the number of segregating sites. First we study the strong-migration limit and prove that the distribution converges to that for a panmictic population. We also study the case of two sampled DNA sequences in the d-dimensional torus model with homogeneous migration. Received 13 July 1995; received in revised form 21 April 1997     

Multilocus dynamics under haploid selection

 A general haploid selection model with arbitrary number of multiallelic loci and arbitrary linkage distribution is considered. The population is supposed to be panmictic. A dynamically equivalent diploid selection model is introduced. There is a position effect in this model if the original haploid selection is not multiplicative. If haploid selection is additive then the fundamental theorem is established even with an estimate for the change in the mean fitness. On this basis exponential convergence to an equilibrium is proved. As rule, the limit states are single-gamete ones. If, moreover, linkage is tight, then the single-gamete state with maximal fitness attracts the population for almost all initial states. Received 27 November 1995; received in revised form 17 January 1996     

Discovery of entry inhibitors for HIV-1 via a new de novo protein design framework

A new (to our knowledge) de novo design framework with a ranking metric based on approximate binding affinity calculations is introduced and applied to the discovery of what we believe are novel HIV-1 entry inhibitors. The framework consists of two stages: a sequence selection stage and a validation stage. The sequence selection stage produces a rank-ordered list of amino-acid sequences by solving an integer programming sequence selection model. The validation stage consists of fold specificity and approximate binding affinity calculations. The designed peptidic inhibitors are 12-amino-acids-long and target the hydrophobic core of gp41. A number of the best-predicted sequences were synthesized and their inhibition of HIV-1 was tested in cell culture. All peptides examined showed inhibitory activity when compared with no drug present, and the novel peptide sequences outperformed the native template sequence used for the design. The best sequence showed micromolar inhibition, which is a 3-15-fold improvement over the native sequence, depending on the donor. In addition, the best sequence equally inhibited wild-type and Enfuvirtide-resistant virus strains.     

研究报告: 基于分子对接和步进序列群的蓖麻毒素核酸适配体序列优化

目的 有效结合分子对接预测和表面等离子体共振实验评价技术，获得亲和力更强、序列最短的最优适配体。方法 针对前期筛选出的靶向蓖麻毒素的3条80 nt单链DNA适配体（L14、P3、L7），在明确各自二维随机区茎环序列与靶蛋白结合能力的基础上，以H-DOCK分子对接为指导，分别确定蓖麻毒素适配体随机区的最短结合单元，从而构建两端延长步进序列群，以表面等离子体共振技术测定序列群序列的亲和力和动力学参数，明确适配体的结合关键结构，从而筛选得到最优适配体。结果 3条全长适配体的随机区适配体L14r、P3r、L7r均可形成一定的茎环结构，其中L14r较L14的亲和力增强9倍、L7r增强2倍、P3r基本不变。对随机区适配体和蓖麻毒素进行分子对接，结果显示，L14r、P3r、L7r的对接分数值皆优于阴性序列40T，结合关键氨基酸个数分别为11、8、9个，存在距离小于5 ?的预测结合位点分别为20、12、15个，具有良好的与蓖麻毒素的结合能力。进一步明确了蓖麻毒素活性口袋所容纳的适配体最短结合单元L14rm、P3rm、L7rm的序列构成，在此基础上构建出两端延长步进序列群。针对该步进群，基于结合关键氨基酸个数、结合位点个数、对接得分等参数的变化和表面等离子体共振测定结果筛选出最优适配体。所获得的最优适配体L14rm、L7rm-2亲和力继续增强了1~2倍。结论 随机区适配体能有效地与蓖麻毒素结合，较之全长适配体亲和力更强，分子对接结合步进序列群设计，仅使用17条序列，便有效获得了3条最优适配体并明确其结合作用。3条结合蓖麻毒素的最优适配体——L14rm、P3r、L7rm-2的K_D值分别为（64±30）、（167±19）、（120±1）nmol/L，亲和力提高到全长适配体的14、1、4倍。