Design of Novel Series of Antimalarial PMX Inhibitors with Increased Half-Life via Molecular Property Optimization

Plasmepsin X (PMX) has been identified as a multistage antimalarial target. PMX is a malarial aspartyl protease essential for merozoite egress from infected red blood cells and invasion of the host erythrocytes. Previously, we reported the identification of PMX inhibitors by structure-based optimization of a cyclic guanidine core. Preclinical assessment of UCB7362, which displayed both in vitro and in vivo antimalarial activity, revealed a suboptimal dose paradigm (once daily dosing of 50 mg for 7 days for treatment of uncomplicated malaria) relative to current standard of care (three-dose regime). We report here the efforts toward extending the half-life (t1/2) by reducing metabolic clearance and increasing volume of distribution (Vss). Our efforts culminated in the identification of a biaryl series, with an expected longer t1/2 in human than UCB7362 while maintaining a similar in vitro off-target hit rate.

Water Networks in Complexes between Proteins and FDA-Approved Drugs

Water molecules at protein-ligand interfaces are often of significant pharmaceutical interest, owing in part to the entropy which can be released upon the displacement of an ordered water by a therapeutic compound. Protein structures may not, however, completely resolve all critical bound water molecules, or there may be no experimental data available. As such, predicting the location of water molecules in the absence of a crystal structure is important in the context of rational drug design. Grand canonical Monte Carlo (GCMC) is a computational technique that is gaining popularity for the simulation of buried water sites. In this work, we assess the ability of GCMC to accurately predict water binding locations, using a dataset that we have curated, containing 108 unique structures of complexes between proteins and Food and Drug Administration (FDA)-approved small-molecule drugs. We show that GCMC correctly predicts 81.4% of nonbulk crystallographic water sites to within 1.4 Å. However, our analysis demonstrates that the reported performance of water prediction methods is highly sensitive to the way in which the performance is measured. We also find that crystallographic water sites with more protein/ligand hydrogen bonds and stronger electron density are more reliably predicted by GCMC. An analysis of water networks revealed that more than half of the structures contain at least one ligand-contacting water network. In these cases, displacement of a water site by a ligand modification might yield unexpected results if the larger network is destabilized. Cooperative effects between waters should therefore be explicitly considered in structure-based drug design.

Open Binding Pose Metadynamics: An Effective Approach for the Ranking of Protein-Ligand Binding Poses

Predicting the correct pose of a ligand binding to a protein and its associated binding affinity is of great importance in computer-aided drug discovery. A number of approaches have been developed to these ends, ranging from the widely used fast molecular docking to the computationally expensive enhanced sampling molecular simulations. In this context, methods such as coarse-grained metadynamics and binding pose metadynamics (BPMD) use simulations with metadynamics biasing to probe the binding affinity without trying to fully converge the binding free energy landscape in order to decrease the computational cost. In BPMD, the metadynamics bias perturbs the ligand away from the initial pose. The resistance of the ligand to this bias is used to calculate a stability score. The method has been shown to be useful in reranking predicted binding poses from docking. Here, we present OpenBPMD, an open-source Python reimplementation and reinterpretation of BPMD. OpenBPMD is powered by the OpenMM simulation engine and uses a revised scoring function. The algorithm was validated by testing it on a wide range of targets and showing that it matches or exceeds the performance of the original BPMD. We also investigated the role of accurate water positioning on the performance of the algorithm and showed how the combination with a grand-canonical Monte Carlo algorithm improves the accuracy of the predictions.

Discovery and Characterization of Potent, Efficacious, Orally Available Antimalarial Plasmepsin X Inhibitors and Preclinical Safety Assessment of

Plasmepsin X (PMX) is an essential aspartyl protease controlling malaria parasite egress and invasion of erythrocytes, development of functional liver merozoites (prophylactic activity), and blocking transmission to mosquitoes, making it a potential multistage drug target. We report the optimization of an aspartyl protease binding scaffold and the discovery of potent, orally active PMX inhibitors with in vivo antimalarial efficacy. Incorporation of safety evaluation early in the characterization of PMX inhibitors precluded compounds with a long human half-life (t_1/2) to be developed. Optimization focused on improving the off-target safety profile led to the identification of UCB7362 that had an improved in vitro and in vivo safety profile but a shorter predicted human t_1/2. UCB7362 is estimated to achieve 9 log 10 unit reduction in asexual blood-stage parasites with once-daily dosing of 50 mg for 7 days. This work demonstrates the potential to deliver PMX inhibitors with in vivo efficacy to treat malaria.

Rings in Clinical Trials and Drugs: Present and Future

We present a comprehensive analysis of all ring systems (both heterocyclic and nonheterocyclic) in clinical trial compounds and FDA-approved drugs. We show 67% of small molecules in clinical trials comprise only ring systems found in marketed drugs, which mirrors previously published findings for newly approved drugs. We also show there are approximately 450 000 unique ring systems derived from 2.24 billion molecules currently available in synthesized chemical space, and molecules in clinical trials utilize only 0.1% of this available pool. Moreover, there are fewer ring systems in drugs compared with those in clinical trials, but this is balanced by the drug ring systems being reused more often. Furthermore, systematic changes of up to two atoms on existing drug and clinical trial ring systems give a set of 3902 future clinical trial ring systems, which are predicted to cover approximately 50% of the novel ring systems entering clinical trials.

Higher Affinity Antibodies Bind With Lower Hydration and Flexibility in Large Scale Simulations

We have carried out a long-timescale simulation study on crystal structures of nine antibody-antigen pairs, in antigen-bound and antibody-only forms, using molecular dynamics with enhanced sampling and an explicit water model to explore interface conformation and hydration. By combining atomic level simulation and replica exchange to enable full protein flexibility, we find significant numbers of bridging water molecules at the antibody-antigen interface. Additionally, a higher proportion of interactions excluding bulk waters and a lower degree of antigen bound CDR conformational sampling are correlated with higher antibody affinity. The CDR sampling supports enthalpically driven antibody binding, as opposed to entropically driven, in that the difference between antigen bound and unbound conformations do not correlate with affinity. We thus propose that interactions with waters and CDR sampling are aspects of the interface that may moderate antibody-antigen binding, and that explicit hydration and CDR flexibility should be considered to improve antibody affinity prediction and computational design workflows.

The Chemical Synthesis of Knob Domain Antibody Fragments

Cysteine-rich knob domains found in the ultralong complementarity determining regions of a subset of bovine antibodies are capable of functioning autonomously as 3-6 kDa peptides. While they can be expressed recombinantly in cellular systems, in this paper we show that knob domains are also readily amenable to a chemical synthesis, with a co-crystal structure of a chemically synthesized knob domain in complex with an antigen showing structural equivalence to the biological product. For drug discovery, following the immunization of cattle, knob domain peptides can be synthesized directly from antibody sequence data, combining the power and diversity of the bovine immune repertoire with the ability to rapidly incorporate nonbiological modifications. We demonstrate that, through rational design with non-natural amino acids, a paratope diversity can be massively expanded, in this case improving the efficacy of an allosteric peptide. As a potential route to further improve stability, we also performed head-to-tail cyclizations, exploiting the proximity of the N and C termini to synthesize functional, fully cyclic antibody fragments. Lastly, we highlight the stability of knob domains in plasma and, through pharmacokinetic studies, use palmitoylation as a route to extend the plasma half-life of knob domains in vivo. This study presents an antibody-derived medicinal chemistry platform, with protocols for solid-phase synthesis of knob domains, together with the characterization of their molecular structures, in vitro pharmacology, and pharmacokinetics.

Water molecules at protein-drug interfaces: computational prediction and analysis methods

The fundamental importance of water molecules at drug-protein interfaces is now widely recognised and a significant feature in structure-based drug design. Experimental methods for analysing the role of water in drug binding have many challenges, including the accurate location of bound water molecules in crystal structures, and problems in resolving specific water contributions to binding thermodynamics. Computational analyses of binding site water molecules provide an alternative, and in principle complete, structural and thermodynamic picture, and their use is now commonplace in the pharmaceutical industry. In this review, we describe the computational methodologies that are available and discuss their strengths and weaknesses. Additionally, we provide a critical analysis of the experimental data used to validate the methods, regarding the type and quality of experimental structural data. We also discuss some of the fundamental difficulties of each method and suggest directions for future study.

Insight into small molecule binding to the neonatal Fc receptor by X-ray crystallography and 100 kHz magic-angle-spinning NMR

Aiming at the design of an allosteric modulator of the neonatal Fc receptor (FcRn)-Immunoglobulin G (IgG) interaction, we developed a new methodology including NMR fragment screening, X-ray crystallography, and magic-angle-spinning (MAS) NMR at 100 kHz after sedimentation, exploiting very fast spinning of the nondeuterated soluble 42 kDa receptor construct to obtain resolved proton-detected 2D and 3D NMR spectra. FcRn plays a crucial role in regulation of IgG and serum albumin catabolism. It is a clinically validated drug target for the treatment of autoimmune diseases caused by pathogenic antibodies via the inhibition of its interaction with IgG. We herein present the discovery of a small molecule that binds into a conserved cavity of the heterodimeric, extracellular domain composed of an α-chain and β2-microglobulin (β2m) (FcRnECD, 373 residues). X-ray crystallography was used alongside NMR at 100 kHz MAS with sedimented soluble protein to explore possibilities for refining the compound as an allosteric modulator. Proton-detected MAS NMR experiments on fully protonated [13C,15N]-labeled FcRnECD yielded ligand-induced chemical-shift perturbations (CSPs) for residues in the binding pocket and allosteric changes close to the interface of the two receptor heterodimers present in the asymmetric unit as well as potentially in the albumin interaction site. X-ray structures with and without ligand suggest the need for an optimized ligand to displace the α-chain with respect to β2m, both of which participate in the FcRnECD-IgG interaction site. Our investigation establishes a method to characterize structurally small molecule binding to nondeuterated large proteins by NMR, even in their glycosylated form, which may prove highly valuable for structure-based drug discovery campaigns.

Stem cell-derived neurons from autistic individuals with SHANK3 mutation show morphogenetic abnormalities during early development

Shank3 is a structural protein found predominantly at the postsynaptic density. Mutations in the SHANK3 gene have been associated with risk for autism spectrum disorder (ASD). We generated induced pluripotent stem cells (iPSCs) from control individuals and from human donors with ASD carrying microdeletions of SHANK3. In addition, we used Zinc finger nucleases to generate isogenic SHANK3 knockout human embryonic stem (ES) cell lines. We differentiated pluripotent cells into either cortical or olfactory placodal neurons. We show that patient-derived placodal neurons make fewer synapses than control cells. Moreover, patient-derived cells display a developmental phenotype: young postmitotic neurons have smaller cell bodies, more extensively branched neurites, and reduced motility compared with controls. These phenotypes were mimicked by SHANK3-edited ES cells and rescued by transduction with a Shank3 expression construct. This developmental phenotype is not observed in the same iPSC lines differentiated into cortical neurons. Therefore, we suggest that SHANK3 has a critical role in neuronal morphogenesis in placodal neurons and that early defects are associated with ASD-associated mutations.

Discovery of a junctional epitope antibody that stabilizes IL-6 and gp80 protein:protein interaction and modulates its downstream signaling

Protein:protein interactions are fundamental in living organism homeostasis. Here we introduce VHH6, a junctional epitope antibody capable of specifically recognizing a neo-epitope when two proteins interact, albeit transiently, to form a complex. Orthogonal biophysical techniques have been used to prove the "junctional epitope" nature of VHH6, a camelid single domain antibody recognizing the IL-6-gp80 complex but not the individual components alone. X-ray crystallography, HDX-MS and SPR analysis confirmed that the CDR regions of VHH6 interact simultaneously with IL-6 and gp80, locking the two proteins together. At the cellular level, VHH6 was able to alter the response of endothelial cells to exogenous IL-6, promoting a sustained STAT3 phosphorylation signal, an accumulation of IL-6 in vesicles and an overall pro-inflammatory phenotype supported further by transcriptomic analysis. Junctional epitope antibodies, like VHH6, not only offer new opportunities in screening and structure-aided drug discovery, but could also be exploited as therapeutics to modulate complex protein:protein interactions.

Computational design of an epitope-specific Keap1 binding antibody using hotspot residues grafting and CDR loop swapping

Therapeutic and diagnostic applications of monoclonal antibodies often require careful selection of binders that recognize specific epitopes on the target molecule to exert a desired modulation of biological function. Here we present a proof-of-concept application for the rational design of an epitope-specific antibody binding with the target protein Keap1, by grafting pre-defined structural interaction patterns from the native binding partner protein, Nrf2, onto geometrically matched positions of a set of antibody scaffolds. The designed antibodies bind to Keap1 and block the Keap1-Nrf2 interaction in an epitope-specific way. One resulting antibody is further optimised to achieve low-nanomolar binding affinity by in silico redesign of the CDRH3 sequences. An X-ray co-crystal structure of one resulting design reveals that the actual binding orientation and interface with Keap1 is very close to the design model, despite an unexpected CDRH3 tilt and V_H/V_L interface deviation, which indicates that the modelling precision may be improved by taking into account simultaneous CDR loops conformation and V_H/V_L orientation optimisation upon antibody sequence change. Our study confirms that, given a pre-existing crystal structure of the target protein-protein interaction, hotspots grafting with CDR loop swapping is an attractive route to the rational design of an antibody targeting a pre-selected epitope.

Efficient Generalized Born Models for Monte Carlo Simulations

The Generalized Born Surface Area theory (GBSA) has become a popular method to model the solvation of biomolecules. While efficient in the context of molecular dynamics simulations, GBSA calculations do not integrate well with Monte Carlo simulations because of the nonlocal nature of the Generalized Born energy. We present a method by which Monte Carlo Generalized Born simulations can be made seven to eight times faster on a protein-ligand binding free energy calculation with little or no loss of accuracy. The method can be employed in any type of Monte Carlo or Hybrid Monte Carlo-molecular dynamics simulation and should prove useful in numerous applications.

Structure-based virtual screening for novel ligands

Computational searches for novel ligands for a given protein binding site have become ubiquitous in the pharmaceutical industry, and are potentially equally useful in helping identify small-molecule tools for biology. Here we describe the steps needed to carry out a high-throughput docking (HTD) or three-dimensional (3D) pharmacophore virtual screen starting with a model of the target protein's structure. The advice given is, in most cases, software independent but some tips are provided which apply only to certain popular programs. Useful work can be carried out using free programs on a modest workstation. Of course, any resultant "hits" remain in the virtual world until they are experimentally tested.

Determination of iron sites and the amount of amorphization in radiation-damaged titanite (CaSiTiO5)

Iron is a ubiquitous impurity in metamict (radiation-damaged and partially amorphized) materials such as titanite (CaSiTiO(5)). Using (57)Fe Mössbauer spectroscopy we find that iron in metamict titanite is partitioned between amorphous and crystalline regions based on valence. Trivalent iron exists in the crystalline titanite matrix whereas divalent iron exists almost exclusively in radiation-amorphized regions. We find that the relative abundances of the oxidation states correlate with the volume fraction of amorphous and crystalline regions. Our data also show that oxidation of iron proceeds along with the recrystallization of the amorphized regions. Recrystallization is confirmed to occur over the range 700 °C < T < 925 °C, and no further structural changes are observed at higher temperatures. It is surprising that our Mössbauer measurements show divalent iron to be surrounded by titanite with a high degree of short-range structural order in the amorphized regions. This observation is fundamentally different from other metamict materials such as zircon (ZrSiO(4)), where amorphized regions show no short-range order.

Dog obesity: owner attitudes and behaviour

Animal (dog) factors that contribute to obesity are classified into three broad categories: genetic pre-disposition, reproductive management and dietary/exercise (human) management. This paper examined the latter-dietary/exercise (human) management. A quantitative analysis of questionnaire responses from dog owners and veterinarians was used to determine the routine care and obesity management strategies for dogs. A total of 550 questionnaires were distributed to dog owners in Victoria, Australia. Owners were asked to score the body condition of their animal by comparison with photographic images of animals with condition score ranging from 2 to 5. The management routines of 219 dog owners were received, including data on 302 dogs. There were 168 households with normal weight animals (BCS 2 and 3) and 51 with obese animals (BCS 4 and 5). The mean number of people per household (normally involved with caring for the animal(s)) with normal weight dogs was significantly lower than that of households with dogs categorised as overweight or obese (Kruskal-Wallis, Chi; chi(2)=6.28; 2.2 (s=0.79) vs. 2.5 (s=1.66); d.f.=2, P<0.05). Dog owners identified a preference for main meal feeding of 'twice a day' (60%), followed by 'once daily' (33%), 'greater than or equal to three times daily' (2%), and 'always feed available' (5%). There was a significant difference (Chi; chi(2)=6.30; d.f.=1; P<0.05) in the frequency of main meal feeding between households. Normal weight animals had food divided into two portions, whereas obese animals or animals from mixed households were more often fed their meal in either one or three-plus portions. Almost all owners fed treats (99%) in the daily feed. Households with normal weight dogs gave treats significantly less frequently than households with obese or mixed weight dogs (Chi; chi(2)=31.81; d.f.=6; P<0.001). The frequency of exercise differed between households (Chi; chi(2)=9.9; d.f.=1; P<0.01), with normal weight dogs being exercised daily compared to weekly for overweight or mixed weight dogs. Owners who indicated that their dog was confined to a yard as its exercise regime rather than walked were also significantly more likely to be an obese or mixed weight household (Chi; chi(2)=18.4; d.f.=1; P<0.001). In conclusion, there is clear evidence that obesity in dogs is affected by the interrelationships between food management, exercise and social factors.

Novel potent BRAF inhibitors: toward 1 nM compounds through optimization of the central phenyl ring

BRAF, a serine/threonine specific protein kinase that is part of the MAPK pathway and acts as a downstream effector of RAS, is a potential therapeutic target in melanoma. We have developed a series of small-molecule BRAF inhibitors based on a 1H-imidazo[4,5-b]pyridine-2(3H)-one scaffold (ring A) as the hinge binding moiety and a number of substituted phenyl rings C that interact with the allosteric binding site. The introduction of various groups on the central phenyl ring B combined with appropriate A- and C-ring modifications afford very potent compounds that inhibit (V600E)BRAF kinase activity in vitro and oncogenic BRAF signaling in melanoma cells. Substitution on the central phenyl ring of a 3-fluoro, a naphthyl, or a 3-thiomethyl group improves activity to yield compounds with an IC(50) of 1 nM for purified (V600E)BRAF and nanomolar activity in cells.

Pyridoimidazolones as novel potent inhibitors of v-Raf murine sarcoma viral oncogene homologue B1 (BRAF)

BRAF is a serine/threonine kinase that is mutated in a range of cancers, including 50-70% of melanomas, and has been validated as a therapeutic target. We have designed and synthesized mutant BRAF inhibitors containing pyridoimidazolone as a new hinge-binding scaffold. Compounds have been obtained which have low nanomolar potency for mutant BRAF (12 nM for compound 5i) and low micromolar cellular potency against a mutant BRAF melanoma cell line, WM266.4. The series benefits from very low metabolism, and pharmacokinetics (PK) that can be modulated by methylation of the NH groups of the imidazolone, resulting in compounds with fewer H-donors and a better PK profile. These compounds have great potential in the treatment of mutant BRAF melanomas.

Prediction of binding modes for ligands in the cytochromes P450 and other heme-containing proteins

The cytochromes P450 (P450s) are a family of heme-containing monooxygenase enzymes involved in a variety of functions, including the metabolism of endogenous and exogenous substances in the human body. During lead optimization, and in drug development, many potential drug candidates are rejected because of the affinity they display for drug-metabolising P450s. Recently, crystal structures of human enzymes involved in drug metabolism have been determined, significantly augmenting the prospect of using structure-based design to modulate the binding and metabolizing properties of compounds against P450 proteins. An important step in the application of structure-based metabolic optimization is the accurate prediction of docking modes in heme binding proteins. In this paper we assess the performance of the docking program GOLD at predicting the binding mode of 45 heme-containing complexes. We achieved success rates of 64% and 57% for Chemscore and Goldscore respectively; these success rates are significantly lower than the value of 79% observed with both scoring functions for the full GOLD validation set. Re-parameterization of metal-acceptor interactions and lipophilicity of planar nitrogen atoms in the scoring functions resulted in a significant increase in the percentage of successful dockings against the heme binding proteins (Chemscore 73%, Goldscore 65%). The modified scoring functions will be useful in docking applications on P450 enzymes and other heme binding proteins.

Origin of the Verwey transition in magnetite

Comprehensive x-ray powder diffraction studies were carried out in magnetite in the 80-150 K and 0-12 GPa ranges with a membrane-driven diamond anvil cell and helium as a pressure medium. Careful data analyses have shown that a reversible, cubic to a distorted-cubic, structural transition takes place with increasing pressure, within the (P,T) regime below the Verwey temperature TV(P). The experimental documentation that TV(P)=Tdist(P) implies that the pressure-temperature-driven metal-insulator Verwey transition is caused by a gap opening in the electronic band structure due to the crystal-structural transformation to a lower-symmetry phase. The distorted-cubic insulating phase comprises a relatively small pressure-temperature range of the stability field of the cubic metallic phase that extends to 25 GPa.

Principal atrial fibrillation discharges by the new ACC/AHA/ESC classification

BACKGROUND

The American College of Cardiology, American Heart Association, and European Society of Cardiology Board (ACC/AHA/ESC) 2001 guidelines for management of patients with atrial fibrillation (AF) include a new classification system consisting of 4 categories: first-detected episode; recurrent paroxysmal (self-terminating); recurrent persistent (requiring cardioversion); and permanent. The frequency of hospital discharges within these categories has not been reported.

METHODS

The new classification system was applied to 135 consecutive hospital discharges with a principal diagnosis of AF.

RESULTS

Classification of AF in these discharged patients included 74 (55%) with first-detected episode; 28 (21%) with recurrent paroxysmal AF; 17 (13%) with recurrent persistent AF; and 16 (12%) with permanent AF. Hypertension (n = 48; 35%) was the most common primary cause of AF, followed by alcohol related (n = 23; 17%), coronary artery disease (n = 20; 15%), and valvular heart disease (n = 17; 12%). For the 102 patients with first-detected and recurrent paroxysmal AF, 71 (69%) converted spontaneously to normal sinus rhythm within 48 hours of admission. Of the 48 patients with a discharge diagnosis of AF, 32 (67%) were receiving anticoagulation therapy.

CONCLUSIONS

Most hospital discharges with a principal diagnosis of AF represent the first-detected episode. Diverse causes contribute to AF, and to examine them would help direct therapy. Importantly, in our analysis, 69% of those patients with first-detected or recurrent paroxysmal AF converted spontaneously to normal sinus rhythm within 48 hours of admission.

Modeling Water Molecules in Protein−Ligand Docking Using GOLD

We implemented a novel approach to score water mediation and displacement in the protein-ligand docking program GOLD. The method allows water molecules to switch on and off and to rotate around their three principal axes. A constant penalty, sigma(p), representing the loss of rigid-body entropy, is added for water molecules that are switched on, hence rewarding water displacement. We tested the methodology in an extensive validation study. First, sigma(p) is optimized against a training set of 58 protein-ligand complexes. For this training set, our algorithm correctly predicts water mediation/displacement in approximately 92% of the cases. We observed small improvements in the quality of the predicted binding modes for water-mediated complexes. In the second part of this work, an entirely independent set of 225 complexes is used. For this test set, our algorithm correctly predicts water mediation/displacement in approximately 93% of the cases. Improvements in binding mode quality were observed for individual water-mediated complexes.

Comparing protein-ligand docking programs is difficult

There is currently great interest in comparing protein-ligand docking programs. A review of recent comparisons shows that it is difficult to draw conclusions of general applicability. Statistical hypothesis testing is required to ensure that differences in pose-prediction success rates and enrichment rates are significant. Numerical measures such as root-mean-square deviation need careful interpretation and may profitably be supplemented by interaction-based measures and visual inspection of dockings. Test sets must be of appropriate diversity and of good experimental reliability. The effects of crystal-packing interactions may be important. The method used for generating starting ligand geometries and positions may have an appreciable effect on docking results. For fair comparison, programs must be given search problems of equal complexity (e.g. binding-site regions of the same size) and approximately equal time in which to solve them. Comparisons based on rescoring require local optimization of the ligand in the space of the new objective function. Re-implementations of published scoring functions may give significantly different results from the originals. Ostensibly minor details in methodology may have a profound influence on headline success rates.

The parameterization and validation of generalized born models using the pairwise descreening approximation

Generalized Born Surface Area (GBSA) models for water using the Pairwise Descreening Approximation (PDA) have been parameterized by two different methods. The first method, similar to that used in previously reported parameterizations, optimizes all parameters against the experimental free energies of hydration of organic molecules. The second method optimizes the PDA parameters to compensate only for systematic errors of the PDA. The best models are compared to Poisson-Boltzmann calculations and applied to the computation of potentials of mean force (PMFs) for the association of various molecules. PMFs present a more rigorous test of the ability of a solvation model to correctly reproduce the screening of intermolecular interactions by the solvent, than its accuracy at predicting free energies of hydration of small molecules. Models derived with the first method are sometimes shown to fail to compute accurate potentials of mean force because of large errors in the computation of Born radii, while no such difficulties are observed with the second method. Furthermore, accurate computation of the Born radii appears to be more important than good agreement with experimental free energies of solvation. We discuss the source of errors in the potentials of mean force and suggest means to reduce them. Our findings suggest that Generalized Born models that use the Pairwise Descreening Approximation and that are derived solely by unconstrained optimization of parameters against free energies of hydration should be applied to the modeling of intermolecular interactions with caution.

Virtual Screening Using Protein−Ligand Docking: Avoiding Artificial Enrichment

This study addresses a number of topical issues around the use of protein-ligand docking in virtual screening. We show that, for the validation of such methods, it is key to use focused libraries (containing compounds with one-dimensional properties, similar to the actives), rather than "random" or "drug-like" libraries to test the actives against. We also show that, to obtain good enrichments, the docking program needs to produce reliable binding modes. We demonstrate how pharmacophores can be used to guide the dockings and improve enrichments, and we compare the performance of three consensus-ranking protocols against ranking based on individual scoring functions. Finally, we show that protein-ligand docking can be an effective aid in the screening for weak, fragment-like binders, which has rapidly become a popular strategy for hit identification. All results presented are based on carefully constructed virtual screening experiments against four targets, using the protein-ligand docking program GOLD.

The influence of whole grain inclusion in pelleted broiler diets on proventricular dilatation and ascites mortality

1. Pelleted diets, incorporating whole or ground wheat or barley in the pellets, were fed to broiler chickens and proventricular dilatation, digesta characteristics and mortality were recorded. 2. The incidence of proventricular dilatation was reduced when whole grain was incorporated into pelleted diets at 200 g/kg. 3. Enzyme application did not influence proventricular dilatation. 4. Digesta contents of intestinal sections were greater and pH lower in anterior intestinal sections when birds had dilated proventriculi. 5. Consideration of overall ascites mortality suggests that the incidence of the condition may be curtailed by whole grain inclusion in pelleted diets. 6. Production efficiency may be enhanced with whole grain inclusion in pellets when feed processing costs and flock health responses are considered.

The incorporation of whole grain into pelleted broiler chicken diets. II. Gastrointestinal and digesta characteristics

1. Pelleted diets, incorporating whole or ground wheat or barley in the pellets, were fed to broiler chickens and the performance, gastrointestinal development and digesta characteristics of those chickens recorded. 2. Body weight was similar with whole grain incorporation or enzyme application. Food conversion efficiency was improved by addition of enzyme to a full wheat diet. 3. Whole grain inclusion in pellets reduced proventriculus proportional mass and increased gizzard proportional mass with no apparent effects of exogenous feed enzyme addition. 4. Relative ileal mass was reduced by enzyme inclusion in a wheat diet and by inclusion of whole barley in the pellets. 5. Digesta viscosity was increased by whole wheat inclusion but reduced by enzyme inclusion irrespective of grain processing. Inclusion of 200 g/kg of barley did not alter viscosity of digesta beyond the duodenum. 6. Fresh excreta pH was higher with both cereals in the grower phase when whole grain was fed. On the barley diet, this was affected by enzyme addition, which created higher pH from a point earlier in the grower phase than when no enzyme was employed. 7. Evidence of a complex interaction between higher viscosity and pH being involved in differences in ileal relative mass was found through significant relationships being produced by enzyme use on the wheat diet alone.

Pressure induced self-oxidation of Fe(OH)2

Mössbauer spectroscopy, x-ray diffraction, and electrical resistance [R(P,T)] studies in Fe(OH)(2) to 40 GPa revealed an unforeseen process by which a gradual Fe2+ oxidation takes place, starting at approximately 8 GPa reaching 70% Fe3+ abundance at 40 GPa. The nonreversible process Fe2+-->Fe3++e(-) occurs with no structural transition. The "ejected" electrons form a deep band within the high-pressure electronic manifold becoming weakly localized at P>50 GPa. This process is attributed to an effective ionization potential created by the pressure induced orientationally deformed (OH) dipoles and the unusual small binding energy of the valence electron in Fe2+(OH)(2).

FDS: flexible ligand and receptor docking with a continuum solvent model and soft-core energy function

The docking of flexible small molecule ligands to large flexible protein targets is addressed in this article using a two-stage simulation-based method. The methodology presented is a hybrid approach where the first component is a dock of the ligand to the protein binding site, based on deriving sets of simultaneously satisfied intermolecular hydrogen bonds using graph theory and a recursive distance geometry algorithm. The output structures are reduced in number by cluster analysis based on distance similarities. These structures are submitted to a modified Monte Carlo algorithm using the AMBER-AA molecular mechanics force field with the Generalized Born/Surface Area (GB/SA) continuum model. This solvent model is not only less expensive than an explicit representation, but also yields increased sampling. Sampling is also increased using a rotamer library to direct some of the protein side-chain movements along with large dihedral moves. Finally, a softening function for the nonbonded force field terms is used, enabling the potential energy function to be slowly turned on throughout the course of the simulation. The docking procedure is optimized, and the results are presented for a single complex of the arabinose binding protein. It was found that for a rigid receptor model, the X-ray binding geometry was reproduced and uniquely identified based on the associated potential energy. However, when side-chain flexibility was included, although the X-ray structure was identified, it was one of three possible binding geometries that were energetically indistinguishable. These results suggest that on relaxing the constraint on receptor flexibility, the docking energy hypersurface changes from being funnel-like to rugged. A further 14 complexes were then examined using the optimized protocol. For each complex the docking methodology was tested for a fully flexible ligand, both with and without protein side-chain flexibility. For the rigid protein docking, 13 out of the 15 test cases were able to find the experimental binding mode; this number was reduced to 11 for the flexible protein docking. However, of these 11, in the majority of cases the experimental binding mode was not uniquely identified, but was present in a cluster of low energy structures that were energetically indistinguishable. These results not only support the presence of a rugged docking energy hypersurface, but also suggest that it may be necessary to consider the possibility of more than one binding conformation during ligand optimization.

Improved protein-ligand docking using GOLD

The Chemscore function was implemented as a scoring function for the protein-ligand docking program GOLD, and its performance compared to the original Goldscore function and two consensus docking protocols, "Goldscore-CS" and "Chemscore-GS," in terms of docking accuracy, prediction of binding affinities, and speed. In the "Goldscore-CS" protocol, dockings produced with the Goldscore function are scored and ranked with the Chemscore function; in the "Chemscore-GS" protocol, dockings produced with the Chemscore function are scored and ranked with the Goldscore function. Comparisons were made for a "clean" set of 224 protein-ligand complexes, and for two subsets of this set, one for which the ligands are "drug-like," the other for which they are "fragment-like." For "drug-like" and "fragment-like" ligands, the docking accuracies obtained with Chemscore and Goldscore functions are similar. For larger ligands, Goldscore gives superior results. Docking with the Chemscore function is up to three times faster than docking with the Goldscore function. Both combined docking protocols give significant improvements in docking accuracy over the use of the Goldscore or Chemscore function alone. "Goldscore-CS" gives success rates of up to 81% (top-ranked GOLD solution within 2.0 A of the experimental binding mode) for the "clean list," but at the cost of long search times. For most virtual screening applications, "Chemscore-GS" seems optimal; search settings that give docking speeds of around 0.25-1.3 min/compound have success rates of about 78% for "drug-like" compounds and 85% for "fragment-like" compounds. In terms of producing binding energy estimates, the Goldscore function appears to perform better than the Chemscore function and the two consensus protocols, particularly for faster search settings. Even at docking speeds of around 1-2 min/compound, the Goldscore function predicts binding energies with a standard deviation of approximately 10.5 kJ/mol.

A review of protein-small molecule docking methods

The binding of small molecule ligands to large protein targets is central to numerous biological processes. The accurate prediction of the binding modes between the ligand and protein, (the docking problem) is of fundamental importance in modern structure-based drug design. An overview of current docking techniques is presented with a description of applications including single docking experiments and the virtual screening of databases.

TANK2, a new TRF1-associated poly(ADP-ribose) polymerase, causes rapid induction of cell death upon overexpression

Tankyrase (TANK1) is a human telomere-associated poly(ADP-ribose) polymerase (PARP) that binds the telomere-binding protein TRF1 and increases telomere length when overexpressed. Here we report characterization of a second human tankyrase, tankyrase 2 (TANK2), which can also interact with TRF1 but has properties distinct from those of TANK1. TANK2 is encoded by a 66-kilobase pair gene (TNKS2) containing 28 exons, which express a 6.7-kilobase pair mRNA and a 1166-amino acid protein. The protein shares 85% amino acid identity with TANK1 in the ankyrin repeat, sterile alpha-motif, and PARP catalytic domains but has a unique N-terminal domain, which is conserved in the murine TNKS2 gene. TANK2 interacted with TRF1 in yeast and in vitro and localized predominantly to a perinuclear region, similar to the properties of TANK1. In contrast to TANK1, however, TANK2 caused rapid cell death when highly overexpressed. TANK2-induced death featured loss of mitochondrial membrane potential, but not PARP1 cleavage, suggesting that TANK2 kills cells by necrosis. The cell death was prevented by the PARP inhibitor 3-aminobenzamide. In vivo, TANK2 may differ from TANK1 in its intrinsic or regulated PARP activity or its substrate specificity.

Structural requirements of Tom40 for assembly into preexisting TOM complexes of mitochondria

Tom40 is the major subunit of the translocase of the outer mitochondrial membrane (the TOM complex). To study the assembly pathway of Tom40, we have followed the integration of the protein into the TOM complex in vitro and in vivo using wild-type and altered versions of the Neurospora crassa Tom40 protein. Upon import into isolated mitochondria, Tom40 precursor proteins lacking the first 20 or the first 40 amino acid residues were assembled as the wild-type protein. In contrast, a Tom40 precursor lacking residues 41 to 60, which contains a highly conserved region of the protein, was arrested at an intermediate stage of assembly. We constructed mutant versions of Tom40 affecting this region and transformed the genes into a sheltered heterokaryon containing a tom40 null nucleus. Homokaryotic strains expressing the mutant Tom40 proteins had growth rate defects and were deficient in their ability to form conidia. Analysis of the TOM complex in these strains by blue native gel electrophoresis revealed alterations in electrophoretic mobility and a tendency to lose Tom40 subunits from the complex. Thus, both in vitro and in vivo studies implicate residues 41 to 60 as containing a sequence required for proper assembly/stability of Tom40 into the TOM complex. Finally, we found that TOM complexes in the mitochondrial outer membrane were capable of exchanging subunits in vitro. A model is proposed for the integration of Tom40 subunits into the TOM complex.

The incorporation of whole grain into pelleted broiler chicken diets: production and physiological responses

1. Pelleted diets, incorporating whole or ground triticale or wheat in the pellets, were fed to broiler chickens and the performance, gastro-intestinal development and health of those chickens recorded. 2. The use of whole triticale in the pelleted food produced similar bodyweight responses to when ground triticale was incorporated in the food. Food conversion efficiency was enhanced when whole triticale was used compared to ground triticale and was similar to when an exogenous enzyme was added to the ground triticale diet. The incorporation of whole wheat into the pelleted food produced similar production responses to the use of ground wheat. 3. The gastro-intestinal development of the birds was greater when they were fed the whole grain diet as gizzard weights were increased by the incorporation of whole grain into the diet. Gastro-intestinal development was thought to be related to the presence of large fibre particles in the food. 4. Proventricular dilatation and mortality due to ascites were reduced by feeding pelleted diets containing whole grain, which was possibly related to the early development of the gastro-intestinal tract relative to total body development. 5. It is suggested that the current use of exogenous enzyme additions to broiler diets may be reduced by incorporating whole grain into pelleted diets for broiler chickens.

Proneness to eating disorders: weightlifters compared to exercisers

Previous findings were confirmed suggesting that weightlifters may be more prone to eating disorders than are casual exercisers. 25 weightlifters' mean eating disorder score was significantly higher than that of 25 casual exercisers when contrasted for practical difference.

Influence of the method of calcium provision on Marek's disease losses in compound and choice fed layers

1. Two strains of layers were fed compound or choice forms of a standard diet and given calcium as either ground limestone or coarse limestone grit daily or every 2nd day. Mortality caused by Marek's disease was monitored during the rearing and laying phases. 2. A dramatic difference in Marek's disease mortality was found between the strains. Heavy Marek's disease mortality in the imported strain A in the rearing phase was not markedly affected by dietary or calcium treatment. 3. Method of feeding and calcium provision interacted to influence significantly the risk of death due to Marek's disease in the laying phase in layer strain A. 4. Choice feeding may adversely influence Marek's disease losses in layers.

U-richness is a defining feature of plant introns and may function as an intron recognition signal in maize

Using a large set of plant gene sequences we compared individual introns to their flanking exons. Both Zea mays and Arabidopsis thaliana introns are U-rich but display no apparent bias for A. We identified fifteen 11-mer U-rich motifs as frequent elements of maize introns, and these are virtually absent from exons. By mutagenesis, we show that the single U-rich motif in the Bronze2 intron of maize plays a key role in intron processing in vivo.

Reconstitution of human telomerase with the template RNA component hTR and the catalytic protein subunit hTRT

The maintenance of chromosome termini, or telomeres, requires the action of the enzyme telomerase, as conventional DNA polymerases cannot fully replicate the ends of linear molecules. Telomerase is expressed and telomere length is maintained in human germ cells and the great majority of primary human tumours. However, telomerase is not detectable in most normal somatic cells; this corresponds to the gradual telomere loss observed with each cell division. It has been proposed that telomere erosion eventually signals entry into senescence or cell crisis and that activation of telomerase is usually required for immortal cell proliferation. In addition to the human telomerase RNA component (hTR; ref. 11), TR1/TLP1 (refs 12, 13), a protein that is homologous to the p80 protein associated with the Tetrahymena enzyme, has been identified in humans. More recently, the human telomerase reverse transcriptase (hTRT; refs 15, 16), which is homologous to the reverse transcriptase (RT)-like proteins associated with the Euplotes aediculatus (Ea_p123), Saccharomyces cerevisiae (Est2p) and Schizosaccharomyces pombe (5pTrt1) telomerases, has been reported to be a telomerase protein subunit. A catalytic function has been demonstrated for Est2p in the RT-like class but not for p80 or its homologues. We now report that in vitro transcription and translation of hTRT when co-synthesized or mixed with hTR reconstitutes telomerase activity that exhibits enzymatic properties like those of the native enzyme. Single amino-acid changes in conserved telomerase-specific and RT motifs reduce or abolish activity, providing direct evidence that hTRT is the catalytic protein component of telomerase. Normal human diploid cells transiently expressing hTRT possessed telomerase activity, demonstrating that hTRT is the limiting component necessary for restoration of telomerase activity in these cells. The ability to reconstitute telomerase permits further analysis of its biochemical and biological roles in cell aging and carcinogenesis.

Kinship support and maternal and adolescent well-being in economically disadvantaged African-American families

This study tested a conceptual model developed to explain the link between kinship support and the psychological well-being of economically disadvantaged African-American adolescents. The relation of kinship support with maternal and adolescent well-being and mothers' child-rearing practices was assessed in 51 African-American families whose incomes placed them at or below the poverty threshold. Findings revealed that kinship social support to mothers/female guardians was positively associated with adolescent psychological well-being, maternal well-being, and more adequate maternal parenting practices (acceptance, firm control and monitoring of behavior, autonomy granting). Maternal well-being and more adequate maternal parenting practices were positively related to adolescent well-being. Evidence of the mediational role of maternal well-being and parenting practices was revealed. When the effects of maternal well-being and maternal parenting practices were controlled, significant relations between kinship support and adolescent well-being were no longer apparent.