A genetic algorithm-based protocol for docking ensembles of small ligands using experimental restraints

Summary A genetic algorithm (GA) based method for docking ensembles of small, flexible ligands to receptor proteins using NMR-derived constraints is described. In this method, three translations and rotations of the ligand and the dihedral angles of the ligand are represented by binary strings and evolve under the genetic operators of cross-over, mutation, migration and selection. The fitness function for the selection process includes distance and dihedral restraints and a repulsive van der Waals term. The GA was applied to a three-atom model system as well as to the streptavidin-biotin complex using simulated intermolecular distance restraints. In both systems, the GA was able to obtain low-energy conformations when only a single binding site was simulated. Calculations were also performed using distance restraints from two distinct binding sites. In these simulations, the GA was able to obtain low-energy conformations corresponding to ligand molecules in each of the two sites. The inclusion of additional ligands in the ensemble did not result in an energetic benefit, confirming that only two ligand conformations were necessary to fulfill the distance restraints. This method allows for a direct investigation of the minimum number of ligand orientations necessary to fulfill experimental distance restraints, and simultaneously yields detailed structural information about each site.     

Stimulation of the labellar sugar receptor of the fleshfly by mono- and disaccharides

Responses of the labellar sugar receptor of the fleshfly, Boettcherisca peregrina, were studied over a wide range of concentrations of several sugars (sucrose, maltose, glucose, fructose, and mannose) in single solutions and in mixtures. The results suggest (a) that the receptor sites are not completely differentiated for glucose and for fructose combination, (b) that the receptor site is composed of two subunits. Such suggestions are based on the classical model, where the response is proportional to the number of the sites, two subunits of each site being simultaneously occupied with one molecule of disaccharides or two molecules of monosaccharides. It is shown, however, that an allosteric model gives a somewhat better interpretation of the experimental results.     

Co-operative regulation mechanism of muscle contraction: Inter-tropomyosin co-operation model

A new model is presented on the basis of our experimental data and the “tropomyosin-blocking theory” of muscle relaxation to explain the regulation of certain characteristics of muscle contraction, namely that the relation of contraction to pCa is co-operative while calcium-binding is essentially non-cooperative. Our experiments show that end-to-end interactions between adjacent tropomyosin molecules in the groove of the actin helix are essential for the co-operative regulation. The blocking theory says that the tropomyosin molecule in relaxed muscle sterically blocks the myosin attachment site on actin, whereas in contracting muscle it moves to a position away from the attachment site. In this model a concerted movement of tropomyosin molecules, brought about by their end-to-end interactions, is considered to be the essential mechanism of co-operative regulation, and it is assumed that the positional changes of tropomyosin occur primarily when the four calcium binding sites of troponin on the tropomyosin are saturated with calcium. Theoretical analysis of the model, based upon the two-state allosteric model, leads to a Michaelis-Menten equation for the Ca-binding function together with a co-operative equation for the state function, proportional to the contraction or ATPase activity. These two functions fit well the experimental data. With cardiac muscle the slope of the contraction versus pCa curve is slightly less steep than that obtained with skeletal muscle. This difference can be explained by the difference in the number of Ca-binding sites of troponins.     

Distance-restrained docking of rifampicin and rifamycin SV to RNA polymerase using systematic FRET measurements: developing benchmarks of model quality and reliability

We are developing distance-restrained docking strategies for modeling macromolecular complexes that combine available high-resolution structures of the components and intercomponent distance restraints derived from systematic fluorescence resonance energy transfer (FRET) measurements. In this article, we consider the problem of docking small-molecule ligands within macromolecular complexes. Using simulated FRET data, we have generated a series of benchmarks that permit estimation of model accuracy based on the quantity and quality of FRET-derived distance restraints, including the number, random error, systematic error, distance distribution, and radial distribution of FRET-derived distance restraints. We find that expected model accuracy is 10 A or better for models based on: i), > or =20 restraints with up to 15% random error and no systematic error, or ii), > or =20 restraints with up to 15% random error, up to 10% systematic error, and a symmetric radial distribution of restraints. Model accuracies can be improved to 5 A or better by increasing the number of restraints to > or =40 and/or by optimizing the distance distribution of restraints. Using experimental FRET data, we have defined the positions of the binding sites within bacterial RNA polymerase of the small-molecule inhibitors rifampicin (Rif) and rifamycin SV (Rif SV). The inferred binding sites for Rif and Rif SV were located with accuracies of, respectively, 7 and 10 A relative to the crystallographically defined binding site for Rif. These accuracies agree with expectations from the benchmark simulations and suffice to indicate that the binding sites for Rif and Rif SV are located within the RNA polymerase active-center cleft, overlapping the binding site for the RNA-DNA hybrid.     

Framework for strategic wind farm site prioritisation based on modelled wolf reproduction habitat in Croatia

In order to meet carbon reduction targets, many nations are greatly expanding their wind power capacity. However, wind farm infrastructure potentially harms wildlife, and we must therefore find ways to balance clean energy demands with the need to protect wildlife. Wide-ranging carnivores live at low density and are particularly susceptible to disturbance from infrastructure development, so are a particular concern in this respect. We focused on Croatia, which holds an important population of wolves and is currently planning to construct many new wind farms. Specifically, we sought to identify an optimal subset of planned wind farms that would meet energy targets while minimising potential impact on wolves. A suitability model for wolf breeding habitat was carried out using Maxent, based on six environmental variables and 31 reproduction site locations collected between 1997 and 2015. Wind farms were prioritised using Marxan to find the optimal trade-off between energy capacity and overlap with critical wolf reproduction habitat. The habitat suitability model predictions were consistent with the current knowledge: probability of wolf breeding site presence increased with distance to settlements, distance to farmland and distance to roads and decreased with distance to forest. Spatial optimisation showed that it would be possible to meet current energy targets with only 31% of currently proposed wind farms, selected in a way that reduces the potential ecological cost (overall predicted wolf breeding site presence within wind farm sites) by 91%. This is a highly efficient outcome, demonstrating the value of this approach for prioritising infrastructure development based on its potential impact on wide-ranging wildlife species.     

Structural map of the dicyclohexylcarbodiimide site of chloroplast coupling factor determined by resonance energy transfer

Fluorescence resonance energy-transfer measurements were made on the membrane-bound chloroplast coupling factor. The distances from the N,N'-dicyclohexylcarbodiimide-binding site on the membrane-bound portion of the enzyme (CF0) to the vesicle surface and to two sulfhydryl sites on the gamma-polypeptide were determined. The dicyclohexylcarbodiimide-binding site was labeled with the fluorescent species N-cyclohexyl-N'-pyrenylcarbodiimide. The vesicle surface was labeled with N-(7-nitro-2,1,3-benzoxadiazol-4-yl)phosphatidylethanolamine. Steady-state energy transfer between the fluorescent-labeled enzyme (energy donor) and varying concentrations of the ethanolamine derivative (energy acceptor) indicated that the distance of closest approach between the energy donor and the outer vesicle surface is 16-24 A. Two specific sites on the gamma-polypeptide were reacted with a coumarinylmaleimide derivative; one is a sulfhydryl that can be labeled only on the thylakoids under energized conditions (the "light" site), while the other is the disulfide site that regulates enzymatic activity. Energy-transfer measurements utilizing steady-state fluorescence and fluorescence lifetime methods indicated that the dicyclohexylcarbodiimide site is approximately 41 A from the light site and approximately 50 A from the gamma-disulfide site. These distances are used to extend the current structural model of the chloroplast coupling factor.     

A new assay method for hydrogenase based on an enzymic electrode reaction. The enzymic electric cell method.

A new assay method for hydrogenase [EC 1.12.2.1] based on the enzymic electrode reaction of H2-H+ equilibrium has been established. The method is based on the experimental fact that the short-circuit current of the electric cell composed of an electrode with hydrogenase and methylviologen as the mediator of H2-H+ equilibrium and a saturated calomel electrode as the counter electrode, is practically proportional to the amount of hydrogenase in the cell. The new method is referred to as the "enzymic electric cell method." This technique has applications not only to routine activity assay but also to the direct determination of the time course of enzyme denaturation, which has not previously been possible.     

A wearable solution for accurate step detection based on the direct measurement of the inter-foot distance

Accurate step detection is crucial for the estimation of gait spatio-temporal parameters. Although several step detection methods based on the use of inertial measurement units (IMUs) have been successfully proposed, they may not perform adequately when the foot is dragged while walking, when walking aids are used, or when walking at low speed. The aim of this study was to test an original step-detection method, the inter-foot distance step counter (IFOD), based on the direct measurement of the distance between feet. Gait data were recorded using a wearable prototype system (SWING^2DS), which integrates an IMU and two time-of-flight distance sensors (DSs). The system was attached to the medial side of the right foot with one DS positioned close to the forefoot (FORE_DS) and the other close to the rearfoot (REAR_DS). Sixteen healthy adults were asked to walk over ground for two minutes along a loop, including both rectilinear and curvilinear portions, during two experimental sessions. The accuracy of the IFOD step counter was assessed using a stereo-photogrammetric system as gold standard. The best performance was obtained for REAR_DS with an accuracy higher than 99.8% for the instrumented foot step and 88.8% for the non-instrumented foot step during both rectilinear and curvilinear walks. Key features of the IFOD step counter are that it is possible to detect both right and left steps by instrumenting one foot only and that it does not rely on foot impact dynamics. The IFOD step counter can be combined with existing IMU-based methods for increasing step-detection accuracy.     

Coarse-Grain Modeling of Shear-Induced Binding between von Willebrand Factor and Collagen

Von Willebrand factor (VWF) is a large multimeric protein that aids in blood clotting. Near injury sites, hydrodynamic force from increased blood flow elongates VWF, exposing binding sites for platelets and collagen. To investigate VWF binding to collagen that is exposed on injured arterial surfaces, Brownian dynamics simulations are performed with a coarse-grain molecular model. Accounting for hydrodynamic interactions in the presence of a stationary surface, shear flow conditions are modeled. Binding between beads in coarse-grain VWF and collagen sites on the surface is described via reversible ligand-receptor-type bond formation, which is governed via Bell model kinetics. For conditions in which binding is energetically favored, the model predicts a high probability for binding at low shear conditions; this is counter to experimental observations but in agreement with what prior modeling studies have revealed. To address this discrepancy, an additional binding criterion that depends on the conformation of a submonomer feature in the model local to a given VWF binding site is implemented. The modified model predicts shear-induced binding, in very good agreement with experimental observations; this is true even for conditions in which binding is significantly favored energetically. Biological implications of the model modification are discussed in terms of mechanisms of VWF activity.     

Transient Response in a Dendritic Neuron Model for Current Injected at One Branch

Mathematical expressions are obtained for the response function corresponding to an instantaneous pulse of current injected to a single dendritic branch in a branched dendritic neuron model. The theoretical model assumes passive membrane properties and the equivalent cylinder constraint on branch diameters. The response function when used in a convolution formula enables one to compute the voltage transient at any specified point in the dendritic tree for an arbitrary current injection at a given input location. A particular numerical example, for a brief current injection at a branch terminal, illustrates the attenuation and delay characteristics of the depolarization peak as it spreads throughout the neuron model. In contrast to the severe attenuation of voltage transients from branch input sites to the soma, the fraction of total input charge actually delivered to the soma and other trees is calculated to be about one-half. This fraction is independent of the input time course. Other numerical examples, which compare a branch terminal input site with a soma input site, demonstrate that, for a given transient current injection, the peak depolarization is not proportional to the input resistance at the injection site and, for a given synaptic conductance transient, the effective synaptic driving potential can be significantly reduced, resulting in less synaptic current flow and charge, for a branch input site. Also, for the synaptic case, the two inputs are compared on the basis of the excitatory post-synaptic potential (EPSP) seen at the soma and the total charge delivered to the soma.     

Food storing by marsh tits

Wild marsh tits (Parus palustris) were allowed to hoard radioactively labelled sunflower seeds, which were subsequently found using a portable scintillation counter. Seeds were stored singly, in various sites close to the ground. Different birds favoured different types of site, although this preference was changeable. Seed density decreased with increasing distance from the feeder, and there was a negative correlation between seed sequence number and the distance it was carried. Hoarded seeds disappeared more rapidly than control seeds in identical sites 100 cm away, suggesting that the birds remember their exact location. Seeds also tended to be stored and recovered in the same sequence. The adaptive significance of these results is discussed.     

Towards understanding the crosstalk between protein post‐translational modifications: Homo‐ and heterotypic PTM pair distances on protein surfaces are not random

Post‐translational modifications (PTMs) represent an important regulatory layer influencing the structure and function of proteins. With broader availability of experimental information on the occurrences of different PTM types, the investigation of a potential “crosstalk” between different PTM types and combinatorial effects have moved into the research focus. Hypothesizing that relevant interferences between different PTM types and sites may become apparent when investigating their mutual physical distances, we performed a systematic survey of pairwise homo‐ and heterotypic distances of seven frequent PTM types considering their sequence and spatial distances in resolved protein structures. We found that actual PTM site distance distributions differ from random distributions with most PTM type pairs exhibiting larger than expected distances with the exception of homotypic phosphorylation site distances and distances between phosphorylation and ubiquitination sites that were found to be closer than expected by chance. Random reference distributions considering canonical acceptor amino acid residues only were found to be shifted to larger distances compared to distances between any amino acid residue type indicating an underlying tendency of PTM‐amenable residue types to be further apart than randomly expected. Distance distributions based on sequence separations were found largely consistent with their spatial counterparts suggesting a primary role of sequence‐based pairwise PTM‐location encoding rather than folding‐mediated effects. Our analysis provides a systematic and comprehensive overview of the characteristics of pairwise PTM site distances on proteins and reveals that, predominantly, PTM sites tend to avoid close proximity with the potential implication that an independent attachment or removal of PTMs remains possible. Proteins 2016; 85:78–92. © 2016 Wiley Periodicals, Inc.     

Some Properties of Site-Specific and General Recombination Inferred from Int-Initiated Exchanges by Bacteriophage Lambda

The site-specific recombination at the attachment site for prophage integration might proceed by two general mechanisms: (1) a concerted reaction without a free intermediate; (2) a sequential mechanism differing from typical general recombination only by an inability of the cross-strand intermediate structure to migrate into the region of nonhomology adjacent to the attachment site. The blocked-migration model predicts frequent genetic exchange in the int xis region near the attachment site if Int-mediated recombination occurs between lambda phage with homologous attachment sites. We find such additional int xis exchanges, but only at very low frequency (1% of the Int-mediated recombination). We conclude that the resolution point only rarely moves away from the initial crossover point specified by Int and, therefore, that the Int reaction is mainly concerted. We interpret the rare additional int xis recombinants as indicative of occasional branch migration from an initial Int-mediated crossover. The frequency of the rare int xis recombinants is not simply related to distance from the attachment site to an int- or xis- mutation, suggesting that the heteroduplex distance is often at least a gene in length. The frequency of these additional exchanges is also not a strong function of distance between two mutations; from this we conclude that the resolution to the observed recombinant structure in the sequential cases occurs often by mismatch repair. We have found no marked effect of mutations in the bacterial recA, recB, recC, recF, or recL genes on the frequency of the int xis recombinants; this may indicate that none of these genes specifies a product uniquely required for resolution of a cross-strand intermediate.     

Relationship Between Aboveground Biomass and Multiple Measures of Biodiversity in Subtropical Forest of Puerto Rico

Anthropogenic activities have accelerated the rate of global loss of biodiversity, making it more important than ever to understand the structure of biodiversity hotspots. One current focus is the relationship between species richness and aboveground biomass (AGB) in a variety of ecosystems. Nonetheless, species diversity, evenness, rarity, or dominance represent other critical attributes of biodiversity and may have associations with AGB that are markedly different than that of species richness. Using data from large trees in four environmentally similar sites in the Luquillo Experimental Forest of Puerto Rico, we determined the shape and strength of relationships between each of five measures of biodiversity (i.e., species richness, Simpson's diversity, Simpson's evenness, rarity, and dominance) and AGB. We quantified these measures of biodiversity using either proportional biomass or proportional abundance as weighting factors. Three of the four sites had a unimodal relationship between species richness and AGB, with only the most mature site evincing a positive, linear relationship. The differences between the mature site and the other sites, as well as the differences between our richness–AGB relationships and those found at other forest sites, highlight the crucial role that prior land use and severe storms have on this forest community. Although the shape and strength of relationships differed greatly among measures of biodiversity and among sites, the strongest relationships within each site were always those involving richness or evenness.     

A possible mechanism for the dynamics of transition between polymerase and exonuclease sites in a high-fidelity DNA polymerase

The fidelity of DNA synthesis by DNA polymerase is significantly increased by a mechanism of proofreading that is performed at the exonuclease active site separate from the polymerase active site. Thus, the transition of DNA between the two active sites is an important activity of DNA polymerase. Here, based on our proposed model, the rates of DNA transition between the two active sites are theoretically studied. With the relevant parameters, which are determined from the available crystal structure and other experimental data, the calculated transfer rate of correctly base-paired DNA from the polymerase to exonuclease sites and the transfer rate after incorporation of a mismatched base are in good agreement with the available experimental data. The transfer rates in the presence of two and three mismatched bases are also consistent with the previous experimental data. In addition, the calculated transfer rate from the exonuclease to polymerase sites has a large value even with the high binding affinity of 3′-5′ ssDNA for the exonuclease site, which is also consistent with the available experimental value. Moreover, we also give some predictive results for the transfer rate of DNA containing only A:T base pairs and that of DNA containing only G:C base pairs.     

A Model for the Vertical Distribution of Flotation-Size Particles

Abstract

An attempt is made to isolate and model the effects of pedoturbation on the vertical distribution of flotation particles within a site matrix. The derived model states that the concentration of particles at a specific depth is directly proportional to the initial quantity deposited and inversely proportional to the square of the distance from the original deposition. The potential value of the model is exemplified through an analysis of differences among the vertical distributions of charcoal, burnt bone, and flakage at a site in central Minnesota.     

Site-specific OH attack to the sugar moiety of DNA: a comparison of experimental data and computational simulation.

Little computational or experimental information is available on site-specific hydroxyl attack probabilities to DNA. In this study, an atomistic stochastic model of OH radical reactions with DNA was developed to compute relative OH attack probabilities at individual deoxyribose hydrogen atoms. A model of the self-complementary decamer duplex d(CCAACGTTGG) was created including Na(+) counter ions and the water molecules of the first hydration layer. Additionally, a method for accounting for steric hindrance from nonreacting atoms was implemented. The model was then used to calculate OH attack probabilities at the various C-H sites of the sugar moiety. Results from this computational model show that OH radicals exhibit preferential attack at different deoxyribose hydrogens, as suggested by their corresponding percentage solvent-accessible surface areas. The percentage OH attack probabilities for the deoxyribose hydrogens [1H(5')+2H(5'), H(4'), H(3'), 1H(2')+2H(2'), H(1')] were calculated as approximately 54.6%, 20.6%, 15.0%, 8.5% and 1.3%, respectively, averaged across the sequence. These results are in good agreement with the latest experimental site-specific DNA strand break data of Balasubramanian et al. [Proc. Natl. Acad. Sci. USA 95, 9738-9742 (1998)]. The data from this stochastic model suggest that steric hindrance from nonreacting atoms significantly influences site-specific hydroxyl radical attack probabilities in DNA. A number of previous DNA damage models have been based on the assumption that C(4') is the preferred site, or perhaps the only site, for OH-mediated DNA damage. However, the results of the present study are in good agreement the experimental results of Balasubramanian et al. in which OH radicals exhibit preferential initial attack at sugar hydrogen atoms in the order 1H(5')+2H(5') > H(4') > H(3') > 1H(2')+2H(2') > H(1').     

Calibrating and testing tissue equivalent proportional counters with37Ar

A method for testing and calibrating tissue equivalent proportional counters with³⁷Ar is described.³⁷Ar is produced by exposure of argon in its normal isotope composition to thermal neutrons. It is shown that - up to volume ratios of 0.01 of argon to the tissue equivalent gas - there is no appreciable effect of the argon admixture on the function of the proportional counter. Conventional calibration methods with characteristic x-rays or with -particles require modifications of the detectors, and they test only small sub-volumes in the counters. In contrast, argon permits calibrations and tests of the resolution that are representative for the entire counter volume and that do not require changes in detector construction. The method is equally applicable to multi-element proportional counters; it is here exemplified by its application to a long cylindrical counter of simplified design that is part of such a multi-element configuration.     

Dispersal of the emerald ash borer, Agrilus planipennis, in newly-colonized sites

• 1 Emerald ash borer Agrilus planipennis Fairmaire (Coleoptera: Buprestidae) is an invasive forest insect pest threatening more than 8 billion ash (Fraxinus spp.) trees in North America. Development of effective survey methods and strategies to slow the spread of A. planipennis requires an understanding of dispersal, particularly in recently established satellite populations.
• 2 We assessed the dispersal of A. planipennis beetles over a single generation at two sites by intensively sampling ash trees at known distances from infested ash logs, the point source of the infestations. Larval density was recorded from more than 100 trees at each site.
• 3 Density of A. planipennis larvae by distance for one site was fit to the Ricker function, inverse power function, and the negative exponential function using a maximum likelihood approach. The prediction of the best model, a negative exponential function, was compared with the results from both sites.
• 4 The present study demonstrates that larval densities rapidly declined with distance, and that most larvae (88.9 and 90.3%) were on trees within 100 m of the emergence point of the adults at each site. The larval distribution pattern observed at both sites was adequately described by the negative exponential function.

     

One parameter to describe the mechanism of splice sites competition

The choice of a splice site is not only related to its own intrinsic strength, but also is influenced by its flanking competitors. Splice site competition is an important mechanism for splice site prediction, especially, it is a new insight for alternative splice site prediction. In this paper, the position weight matrix scoring function is used to represent splice site strength, and the mechanism of splice site competition is described by only one parameter: scoring function subtraction. While applying on the alternative splice site prediction, based on the only one parameter, 68.22% of donor sites and 70.86% of acceptor sites are correctly classified into alternative and constitutive. The prediction abilities are approximately equal to the recent method which is based on the mechanism of splice site competition. The results reveal that the scoring function subtraction is the best parameter to describe the mechanism of splice sites competition.