Hepatocyte determinants of susceptibility to pre-erythrocytic Plasmodium infection

Thesis (Ph.D.)--University of Washington, 2016-08

DATA DRIVEN APPROACHES TO IDENTIFY DETERMINANTS OF HEART DISEASES AND CANCER RESISTANCE

Cancer and cardio-vascular diseases are the leading causes of death world-wide. Caused by systemic genetic and molecular disruptions in cells, these disorders are the manifestation of profound disturbance of normal cellular homeostasis. People suffering or at high risk for these disorders need early diagnosis and personalized therapeutic intervention. Successful implementation of such clinical measures can significantly improve global health. However, development of effective therapies is hindered by the challenges in identifying genetic and molecular determinants of the onset of diseases; and in cases where therapies already exist, the main challenge is to identify molecular determinants that drive resistance to the therapies. Due to the progress in sequencing technologies, the access to a large genome-wide biological data is now extended far beyond few experimental labs to the global research community. The unprecedented availability of the data has revolutionized the capabilities of computational researchers, enabling them to collaboratively address the long standing problems from many different perspectives. Likewise, this thesis tackles the two main public health related challenges using data driven approaches. Numerous association studies have been proposed to identify genomic variants that determine disease. However, their clinical utility remains limited due to their inability to distinguish causal variants from associated variants. In the presented thesis, we first propose a simple scheme that improves association studies in supervised fashion and has shown its applicability in identifying genomic regulatory variants associated with hypertension. Next, we propose a coupled Bayesian regression approach -- eQTeL, which leverages epigenetic data to estimate regulatory and gene interaction potential, and identifies combinations of regulatory genomic variants that explain the gene expression variance. On human heart data, eQTeL not only explains a significantly greater proportion of expression variance in samples, but also predicts gene expression more accurately than other methods. We demonstrate that eQTeL accurately detects causal regulatory SNPs by simulation, particularly those with small effect sizes. Using various functional data, we show that SNPs detected by eQTeL are enriched for allele-specific protein binding and histone modifications, which potentially disrupt binding of core cardiac transcription factors and are spatially proximal to their target. eQTeL SNPs capture a substantial proportion of genetic determinants of expression variance and we estimate that 58% of these SNPs are putatively causal. The challenge of identifying molecular determinants of cancer resistance so far could only be dealt with labor intensive and costly experimental studies, and in case of experimental drugs such studies are infeasible. Here we take a fundamentally different data driven approach to understand the evolving landscape of emerging resistance. We introduce a novel class of genetic interactions termed synthetic rescues (SR) in cancer, which denotes a functional interaction between two genes where a change in the activity of one vulnerable gene (which may be a target of a cancer drug) is lethal, but subsequently altered activity of its partner rescuer gene restores cell viability. Next we describe a comprehensive computational framework --termed INCISOR-- for identifying SR underlying cancer resistance. Applying INCISOR to mine The Cancer Genome Atlas (TCGA), a large collection of cancer patient data, we identified the first pan-cancer SR networks, composed of interactions common to many cancer types. We experimentally test and validate a subset of these interactions involving the master regulator gene mTOR. We find that rescuer genes become increasingly activated as breast cancer progresses, testifying to pervasive ongoing rescue processes. We show that SRs can be utilized to successfully predict patients' survival and response to the majority of current cancer drugs, and importantly, for predicting the emergence of drug resistance from the initial tumor biopsy. Our analysis suggests a potential new strategy for enhancing the effectiveness of existing cancer therapies by targeting their rescuer genes to counteract resistance. The thesis provides statistical frameworks that can harness ever increasing high throughput genomic data to address challenges in determining the molecular underpinnings of hypertension, cardiovascular disease and cancer resistance. We discover novel molecular mechanistic insights that will advance the progress in early disease prevention and personalized therapeutics. Our analyses sheds light on the fundamental biological understanding of gene regulation and interaction, and opens up exciting avenues of translational applications in risk prediction and therapeutics.

Coarse Graining to Invesitigate Peptide-Lipid Interactions

Experimental characterization of molecular details is challenging, and although single molecule experiments have gained prominence, oligomer characterization remains largely unexplored. The ability to monitor the time evolution of individual molecules while they self assemble is essential in providing mechanistic insights about biological events. Molecular dynamics (MD) simulations can fill the gap in knowledge between single molecule experiments and ensemble studies like NMR, and are increasingly used to gain a better understanding of microscopic properties. Coarse-grained (CG) models aid in both exploring longer length and time scale molecular phenomena, and narrowing down the key interactions responsible for significant system characteristics. Over the past decade, CG techniques have made a significant impact in understanding physicochemical processes. However, the realm of peptide-lipid interfacial interactions, primarily binding, partitioning and folding of amphipathic peptides, remains largely unexplored compared to peptide folding in solution. The main drawback of existing CG models is the inability to capture environmentally sensitive changes in dipolar interactions, which are indigenous to protein folding, and lipid dynamics. We have used the Drude oscillator approach to incorporate structural polarization and dipolar interactions in CG beads to develop a minimalistic peptide model, WEPPROM (Water Explicit Polarizable PROtein Model), and a lipid model WEPMEM (Water Explicit Polarizable MEmbrane Model). The addition of backbone dipolar interactions in a CG model for peptides enabled us to achieve alpha-beta secondary structure content de novo, without any added bias. As a prelude to studying amphipathic peptide-lipid membrane interactions, the balance between hydrophobicity and backbone dipolar interactions in driving ordered peptide aggregation in water and at a hydrophobic-hydrophilic interface, was explored. We found that backbone dipole interactions play a crucial role in driving ordered peptide aggregation, both in water and at hydrophobic-hydrophilic interfaces; while hydrophobicity is more relevant for aggregation in water. A zwitterionic (POPC: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and an anionic lipid (POPS: 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine) are used as model lipids for WEPMEM. The addition of head group dipolar interactions in lipids significantly improved structural, dynamic and dielectric properties of the model bilayer. Using WEPMEM and WEPPROM, we studied membrane-induced peptide folding of a cationic antimicrobial peptide with anticancer activity, SVS-1. We found that membrane-induced peptide folding is driven by both (a) cooperativity in peptide self interaction and (b) cooperativity in membrane-peptide interactions. The dipolar interactions between the peptide and the lipid head-groups contribute to stabilizing folded conformations. The role of monovalent ion size and peptide concentration in driving lipid domain formation in anionic/zwitterionic lipid mixtures was also investigated. Our study suggest monovalent ion size to be a crucial determinant of interaction with lipid head groups, and hence domain formation in lipid mixtures. This study reinforces the role of dipole interactions in protein folding, lipid membrane properties, membrane induced peptide folding and lipid domain formation. Therefore, the models developed in this thesis can be used to explore a multitude of biomolecular processes, both at longer time-scales and larger system sizes.

THE ROLE OF THE MECHANICAL ENVIRONMENT ON CANCER CELL TRANSMIGRATION AND MRNA LOCALIZATION

Most cancer-related deaths are due to metastasis formation, the ability of cancer cells to break away from the primary tumor site, transmigrate through the endothelium, and form secondary tumors in distant areas. Many studies have identified links between the mechanical properties of the cellular microenvironment and the behavior of cancer cells. Cells may experience heterogeneous microenvironments of varying stiffness during tumor progression, transmigration, and invasion into the basement membrane. In addition to mechanical factors, the localization of RNAs to lamellipodial regions has been proposed to play an important part in metastasis. This dissertation provides a quantitative evaluation of the biophysical effects on cancer cell transmigration and RNA localization. In the first part of this dissertation, we sought to compare cancer cell and leukocyte transmigration and investigate the impact of matrix stiffness on transmigration process. We found that cancer cell transmigration includes an additional step, ‘incorporation’, into the endothelial cell (EC) monolayer. During this phase, cancer cells physically displace ECs and spread into the monolayer. Furthermore, the effects of subendothelial matrix stiffness and endothelial activation on cancer cell incorporation are cell-specific, a notable difference from the process by which leukocytes transmigrate. Collectively, our results provide mechanistic insights into tumor cell extravasation and demonstrate that incorporation into the endothelium is one of the earliest steps. In the next part of this work, we investigated how matrix stiffness impacts RNA localization and its relevance to cancer metastasis. In migrating cells, the tumor suppressor protein, adenomatous polyposis coli (APC) targets RNAs to cellular protrusions. We observed that increasing stiffness promotes the peripheral localization of these APC-dependent RNAs and that cellular contractility plays a role in regulating this pathway. We next investigated the mechanism underlying the effect of substrate stiffness and cellular contractility. We found that contractility drives localization of RNAs to protrusions through modulation of detyrosinated microtubules, a network of modified microtubules that associate with, and are required for localization of APC-dependent RNAs. These results raise the possibility that as the matrix environment becomes stiffer during tumor progression, it promotes the localization of RNAs and ultimately induces a metastatic phenotype.

Plant-derived compounds stimulate the decomposition of organic matter in arctic permafrost soils

Arctic ecosystems are warming rapidly, which is expected to promote soil organic matter (SOM) decomposition. In addition to the direct warming effect, decomposition can also be indirectly stimulated via increased plant productivity and plant-soil C allocation, and this so called "priming effect" might significantly alter the ecosystem C balance. In this study, we provide first mechanistic insights into the susceptibility of SOM decomposition in arctic permafrost soils to priming. By comparing 119 soils from four locations across the Siberian Arctic that cover all horizons of active layer and upper permafrost, we found that an increased availability of plant-derived organic C particularly stimulated decomposition in subsoil horizons where most of the arctic soil carbon is located. Considering the 1,035 Pg of arctic soil carbon, such an additional stimulation of decomposition beyond the direct temperature effect can accelerate net ecosystem C losses, and amplify the positive feedback to global warming.

Mechanistic features of Salmonella typhimurium propionate kinase (TdcD): Insights from kinetic and crystallographic studies

Short-chain fatty acids (SCFAs) play a major role in carbon cycle and can be utilized as a source of carbon and energy by bacteria. Salmonella typhimurium propionate kinase (StTdcD) catalyzes reversible transfer of the gamma-phosphate of ATP to propionate during L-threonine degradation to propionate. Kinetic analysis revealed that StTdcD possesses broad ligand specificity and could be activated by various SCFAs (propionate > acetate approximate to butyrate), nucleotides (ATP approximate to GTP > CTP approximate to TTP; dATP > dGTP > dCTP) and metal ions (Mg2+ approximate to Mn2+ > Co2+). Inhibition of StTdcD by tricarboxylic acid (TCA) cycle intermediates such as citrate, succinate, alpha-ketoglutarate and malate suggests that the enzyme could be under plausible feedback regulation. Crystal structures of StTdcD bound to PO4 (phosphate), AMP, ATP, Ap4 (adenosine tetraphosphate), GMP, GDP, GTP, CMP and CTP revealed that binding of nucleotide mainly involves hydrophobic interactions with the base moiety and could account for the broad biochemical specificity observed between the enzyme and nucleotides. Modeling and site-directed mutagenesis studies suggest Ala88 to be an important residue involved in determining the rate of catalysis with SCFA substrates. Molecular dynamics simulations on monomeric and dimeric forms of StTdcD revealed plausible open and closed states, and also suggested role for dimerization in stabilizing segment 235-290 involved in interfacial interactions and ligand binding. Observation of an ethylene glycol molecule bound sufficiently close to the gamma-phosphate in StTdcD complexes with triphosphate nucleotides supports direct in-line phosphoryl transfer. (C) 2013 Elsevier B.V. All rights reserved.

New insights on nevirapine use: A mechanistic perspective of its toxic events

Introdução: A nevirapina (NVP) é um fármaco amplamente utilizado para o tratamento da infecção pelo vírus da imunodeficiência humana de tipo 1 (VIH-1), no entanto, a sua utilização na terapêutica crónica tem sido associada à toxicidade hepática e cutânea. O sexo feminino é um factor de risco para o desenvolvimento destes eventos tóxicos, mas as razões para essa diferença entre o sexo feminino e masculino não estão completamente esclarecidas. Diferenças na biotransformação da NVP e na formação de metabolitos tóxicos podem ser as causas subjacentes. O presente trabalho teve como objectivo explorar as diferenças entre homens e mulheres na biotransformação da NVP, como um potencial factor de toxicidade induzida por este fármaco anti-retroviral. Materiais e Métodos: Todos os indivíduos incluídos no presente estudo eram adultos com infecção por VIH-1 confirmada, tratados com 400 mg de NVP uma vez ao dia, durante pelo menos 1 mês. Foram colhidas amostras de sangue e os níveis de NVP e dos metabolitos de fase I foram determinados por cromatografia líquida de alta performance. Os dados antropométricos e clínicos e os perfis de metabolitos foram avaliados de forma a averiguar possíveis diferenças relacionadas com o sexo dos indivíduos. Resultados: Foram incluídos 52 doentes (63% do sexo masculino). O peso corporal foi inferior nas mulheres (p = 0.028) e o sexo feminino foi associado a maiores níveis de fosfatase alcalina (p = 0.036) e lactato desidrogenase (p = 0.037). Os níveis plasmáticos de NVP (p = 0.030) e 3-hidroxi-NVP (p = 0.035), assim como as proporções de 12-hidroxi-NVP (p = 0.037) e 3-hidroxi-NVP (p = 0.001) foram maiores nas mulheres, quando ajustados pelo peso corporal dos indivíduos. Discussão: Existem diferenças na biotransformação da NVP entre homens e mulheres, particularmente na formação de 12-hidroxi-NVP e 3-hidroxi-NVP. Estes resultados apontam para uma formação de metabolitos reactivos, que é dependente do sexo e que pode contribuir para o perfil de dimorfismo sexual associado às reacções tóxicas induzidas pela NVP.

RNA-based gene duplication: mechanistic and evolutionary insights.

Gene copies that stem from the mRNAs of parental source genes have long been viewed as evolutionary dead-ends with little biological relevance. Here we review a range of recent studies that have unveiled a significant number of functional retroposed gene copies in both mammalian and some non-mammalian genomes. These studies have not only revealed previously unknown mechanisms for the emergence of new genes and their functions but have also provided fascinating general insights into molecular and evolutionary processes that have shaped genomes. For example, analyses of chromosomal gene movement patterns via RNA-based gene duplication have shed fresh light on the evolutionary origin and biology of our sex chromosomes.

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background: Bacteria such as Escherichia coli and Salmonella typhimurium can utilize acetate as the sole source of carbon and energy. Acetate kinase (AckA) and phosphotransacetylase (Pta), key enzymes of acetate utilization pathway, regulate flux of metabolites in glycolysis, gluconeogenesis, TCA cycle, glyoxylate bypass and fatty acid metabolism. Results: Here we report kinetic characterization of S. typhimurium AckA (StAckA) and structures of its unliganded (Form-I, 2.70 angstrom resolution) and citrate-bound (Form-II, 1.90 angstrom resolution) forms. The enzyme showed broad substrate specificity with k(cat)/K-m in the order of acetate > propionate > formate. Further, the K-m for acetyl-phosphate was significantly lower than for acetate and the enzyme could catalyze the reverse reaction (i.e. ATP synthesis) more efficiently. ATP and Mg2+ could be substituted by other nucleoside 5'-triphosphates (GTP, UTP and CTP) and divalent cations (Mn2+ and Co2+), respectively. Form-I StAckA represents the first structural report of an unliganded AckA. StAckA protomer consists of two domains with characteristic beta beta beta alpha beta alpha beta alpha topology of ASKHA superfamily of proteins. These domains adopt an intermediate conformation compared to that of open and closed forms of ligand-bound Methanosarcina thermophila AckA (MtAckA). Spectroscopic and structural analyses of StAckA further suggested occurrence of inter-domain motion upon ligand-binding. Unexpectedly, Form-II StAckA structure showed a drastic change in the conformation of residues 230-300 compared to that of Form-I. Further investigation revealed electron density corresponding to a citrate molecule in a pocket located at the dimeric interface of Form-II StAckA. Interestingly, a similar dimeric interface pocket lined with largely conserved residues could be identified in Form-I StAckA as well as in other enzymes homologous to AckA suggesting that ligand binding at this pocket may influence the function of these enzymes. Conclusions: The biochemical and structural characterization of StAckA reported here provides insights into the biochemical specificity, overall fold, thermal stability, molecular basis of ligand binding and inter-domain motion in AckA family of enzymes. Dramatic conformational differences observed between unliganded and citrate-bound forms of StAckA led to identification of a putative ligand-binding pocket at the dimeric interface of StAckA with implications for enzymatic function.

New Insights into Electronic and Geometric Effects in the Enhanced Photoelectrooxidation of Ethanol Using ZnO Nanorod/Ultrathin Au Nanowire Hybrids

Oxidation of small organic molecules in a fuel cell is a viable method for energy production. However, the key issue is the development of suitable catalysts that exhibit high efficiencies and remain stable during operation. Here, we demonstrate that amine-modified ZnO nanorods on which ultrathin Au nanowires are grown act as an excellent catalyst for the oxidation of ethanol. We show that the modification of the ZnO nanorods with oleylamine not only modifies the electronic structure favorably but also serves to anchor the Au nanowires on the nanorods. The adsorption of OH- species on the Au nanowires that is essential for ethanol oxidation is facilitated at much lower potentials as compared to bare Au nanowires leading to high activity. While ZnO shows negligible electrocatalytic activity under normal conditions, there is significant enhancement in the activity under light irradiation. We demonstrate a synergistic enhancement in the photoelectrocatalytic activity of the ZnO/Au nanowire hybrid and provide mechanistic explanation for this enhancement based on both electronic as well as geometric effects. The principles developed are applicable for tuning the properties of other metal/semiconductor hybrids with potentially interesting applications beyond the fuel cell application demonstrated here.

Insights into the Structural Dynamics of Nucleocytoplasmic Transport of tRNA by Exportin-t

Exportin-t (Xpot) transports mature 5'- and 3'-end processed tRNA from the nucleus to the cytoplasm by associating with a small G-protein Ran (RAs-related nuclear protein), in the nucleus. The release of tRNA in cytoplasm involves RanGTP hydrolysis. Despite the availability of crystal structures of nuclear and cytosolic forms of Xpot, the molecular details regarding the sequential events leading to tRNA release and subsequent conformational changes occurring in Xpot remain unknown. We have performed a combination of classical all-atom and accelerated molecular dynamics simulations on a set of complexes involving Xpot to study a range of features including conformational flexibility of free and cargo-bound Xpot and functionally critical contacts between Xpot and its cargo. The systems investigated include free Xpot and its different complexes, bound either to Ran (GTP/GDP) or tRNA or both. This approach provided a statistically reliable estimate of structural dynamics of Xpot after cargo release. The mechanistic basis for Xpot opening after cargo release has been explained in terms of dynamic structural hinges, about which neighboring region could be displaced to facilitate the nuclear to cytosolic state transition. Post-RanGTP hydrolysis, a cascade of events including local conformational change in RanGTP and loss of critical contacts at Xpot/tRNA interface suggest factors responsible for eventual release of tRNA. The level of flexibility in different Xpot complexes varied depending on the arrangement of individual HEAT repeats. Current study provides one of the most comprehensive and robust analysis carried out on this protein using molecular dynamics schemes.

Muscle signaling in exercise intolerance: Insights from the McArdle mouse model

We recently generated a knock-in mouse model (PYGM p.R50X/p.R50X) of McArdle disease (myophosphorylase deficiency). One mechanistic approach to unveil the molecular alterations caused by myophosphorylase deficiency, which is arguably the paradigm of 'exercise intolerance', is to compare the skeletal-muscle tissue of McArdle, heterozygous, and healthy (wild type (wt)) mice. We analyzed in quadriceps muscle of p.R50X/p.R50X (n=4), p.R50X/wt (n=6) and wt/wt mice (n=5) (all male, 8 wk-old) molecular markers of energy-sensing pathways, oxidative phosphorylation (OXPHOS) and autophagy/proteasome systems, oxidative damage and sarcoplamic reticulum (SR) Ca handling. We found a significant group effect for total AMPK (tAMPK) and ratio of phosphorylated (pAMPK)/tAMPK (P=0.012 and 0.033), with higher mean values in p.R50X/p.R50X mice vs. the other two groups. The absence of massive accumulation of ubiquitinated proteins, autophagosomes or lysosomes in p.R50X/p.R50X mice suggested no major alterations in autophagy/proteasome systems. Citrate synthase activity was lower in p.R50X/p.R50X mice vs. the other two groups (P=0.036) but no statistical effect existed for respiratory chain complexes. We found higher levels of 4-hydroxy-2-nonenal-modified proteins in p.R50X/p.R50X and p.R50X/wt mice compared with the wt/wt group (P=0.011). Sarco(endo)plasmic reticulum ATPase 1 (SERCA1) levels detected at 110kDa tended to be higher in p.R50X/p.R50X and p.R50X/wt mice compared with wt/wt animals (P=0.076), but their enzyme activity was normal. We also found an accumulation of phosphorylated SERCA1 in p.R50X/p.R50X animals. Myophosphorylase deficiency causes alterations in sensory energetic pathways together with some evidence of oxidative damage and alterations in Ca handling but with no major alterations in OXPHOS capacity or autophagy/ubiquitination pathways, which suggests that the muscle tissue of patients is likely to adapt overall favorably to exercise training interventions.

Synthesis and enzymatic evaluation of the guanosine analogue 2-amino-6-mercapto-7-methylpurine ribonucleoside (MESG): insights into the phosphorolysis reaction mechanism based on the blueprint transition state: SN1 or SN2?

A modified experimental procedure for the synthesis of MESG (2-amino-6-mercapto-7-methylpurine ribonucleoside) 1 has been successfully performed and its full characterization is presented. High resolution ESI(+)-MSMS indicates both the nucleoside bond cleavage as the main fragmentation in the gas phase and a possible SN1 mechanism. Ab initio transition state calculations based on the blue print transition state support this mechanistic rationale and discard an alternative SN2 mechanism. Assays using purine nucleoside phosphorylase (PNP) enzyme (human and M. tuberculosis sources) indicate its efficiency in the phosphorolysis of MESG and allow the quantitative determination of inorganic phosphate in real time assay.

Structural insights into the allosteric activation mechanism of cGMP-dependent protein kinase I

Cyclic GMP-dependent protein kinase (PKG) is a key transducer in the NO-cGMP signaling pathway. In this line, PKG has been considered an important drug target for treating hypertensive cardiovascular and pulmonary diseases. However, the investigation of PKG’s allosteric activation mechanism has been hampered by a lack of structural information. One of the fundamental questions on the cGMP-dependent activation of PKG is how the enzyme can distinguish cGMP over cAMP and selectively respond to cGMP. To ensure proper signaling, PKG must have developed unique features to ensure its activation upon the right activation signal. In this thesis, the cGMP-selective activation mechanism of PKG was studied through determining crystal structures of three truncated constructs of the regulatory domain [CNB-A (92-227), CNB-B (271-369), and CNB-A/B (92-351)] of PKG Iβ in the absence or presence of cyclic nucleotides. Herein, two individual CNB domain structures with biochemical data revealed that the C-terminal CNB domain (CNB-B) is responsible for cGMP selectivity, while the N-terminal CNB-domain (CNB-A) has a higher binding affinity for both cGMP and cAMP without showing any selectivity. Based on these crystal structures, mutagenesis studies were performed in which the critical residues for cyclic nucleotide selectivity and activation were identified. Furthermore, we discovered that the conformational changes of the C-terminal helix of the CNB-B that bridges between the regulatory and catalytic domains including the hydrophobic capping interaction are crucial for PKG activation. In addition, to observe the global conformation of the activated R-domain, I solved a co-crystal structure of the CNB-A/B with cGMP. Although a monomeric construct was crystallized, the structure displays a dimer. Strikingly, the CNB-A domain and its bound cGMP provide a key interface for this dimeric interaction. Using small angle X-ray scattering (SAXS), the existence of the cGMP-mediated dimeric interface within the CNB domains was confirmed. Furthermore, measuring cGMP-binding affinities (EC50) of the dimeric interface mutants as well as determining activation constants (Ka) revealed that the interface formation is important for PKG activation. To conclude, this thesis study provides a new mechanistic insight in PKG activation along with a newly found interface that can be targeted for designing PKG-specific activity modulators.

Effects of agricultural management practices on earthworm populations and crop yield: validation and application of a mechanistic modelling approach

There is little consensus on how agriculture will meet future food demands sustainably. Soils and their biota play a crucial role by mediating ecosystem services that support agricultural productivity. However, a multitude of site-specific environmental factors and management practices interact to affect the ability of soil biota to perform vital functions, confounding the interpretation of results from experimental approaches. Insights can be gained through models, which integrate the physiological, biological and ecological mechanisms underpinning soil functions. We present a powerful modelling approach for predicting how agricultural management practices (pesticide applications and tillage) affect soil functioning through earthworm populations. By combining energy budgets and individual-based simulation models, and integrating key behavioural and ecological drivers, we accurately predict population responses to pesticide applications in different climatic conditions. We use the model to analyse the ecological consequences of different weed management practices. Our results demonstrate that an important link between agricultural management (herbicide applications and zero, reduced and conventional tillage) and earthworms is the maintenance of soil organic matter (SOM). We show how zero and reduced tillage practices can increase crop yields while preserving natural ecosystem functions. This demonstrates how management practices which aim to sustain agricultural productivity should account for their effects on earthworm populations, as their proliferation stimulates agricultural productivity. Synthesis and applications. Our results indicate that conventional tillage practices have longer term effects on soil biota than pesticide control, if the pesticide has a short dissipation time. The risk of earthworm populations becoming exposed to toxic pesticides will be reduced under dry soil conditions. Similarly, an increase in soil organic matter could increase the recovery rate of earthworm populations. However, effects are not necessarily additive and the impact of different management practices on earthworms depends on their timing and the prevailing environmental conditions. Our model can be used to determine which combinations of crop management practices and climatic conditions pose least overall risk to earthworm populations. Linking our model mechanistically to crop yield models would aid the optimization of crop management systems by exploring the trade-off between different ecosystem services.