SUMMARY: Network medicine leverages the quantification of information flow within sub-cellular networks to elucidate disease etiology and comorbidity, as well as to predict drug efficacy and identify potential therapeutic targets. However, current Network Medicine toolsets often lack computationally efficient data processing pipelines that support diverse scoring functions, network distance metrics, and null models. These limitations hamper their application in large-scale molecular screening, hypothesis testing, and ensemble modeling. To address these challenges, we introduce NetMedPy, a highly efficient and versatile computational package designed for comprehensive Network Medicine analyses. AVAILABILITY AND IMPLEMENTATION: NetMedPy is an open-source Python package under an MIT license. Source code, documentation, and installation instructions can be downloaded from https://github.com/menicgiulia/NetMedPy and https://pypi.org/project/NetMedPy. The package can run on any standard desktop computer or computing cluster.

MOTIVATION: Disease gene prioritization methods assign scores to genes or proteins according to their likely relevance for a given disease based on a provided set of seed genes. This scoring can be used to find new biologically relevant genes or proteins for many diseases. Although methods based on classical random walks have proven to yield competitive results, quantum walk methods have not been explored to this end. RESULTS: We propose a new algorithm for disease gene prioritization based on continuous-time quantum walks using the adjacency matrix of a protein-protein interaction (PPI) network. We demonstrate the success of our proposed quantum walk method by comparing it to several well-known gene prioritization methods on three disease sets, across seven different PPI networks. In order to compare these methods, we use cross-validation and examine the mean reciprocal ranks of recall and average precision values. We further validate our method by performing an enrichment analysis of the predicted genes for coronary artery disease. AVAILABILITY AND IMPLEMENTATION: The data and code for the methods can be accessed at https://github.com/markgolds/qdgp.

The physiochemical determinants of drug-target interactions in the microenvironment of the cell are complex and generally not defined by simple diffusion and intrinsic chemical reactivity. Non-specific interactions of drugs and macromolecules in cells are rarely considered formally in assessing pharmacodynamics. Here, we demonstrate that non-specific interactions lead to very slow incorporation kinetics of DNA binding drugs. We observe a rate of drug incorporation in cell nuclei three orders of magnitude slower than in vitro due to anomalous drug diffusion within cells. This slow diffusion, however, has an advantageous consequence: it leads to virtually irreversible binding of the drug to specific DNA targets in cells. We show that non-specific interactions drive slow drug diffusion manifesting as slow reaction front propagation. We study the effect of non-specific interactions in different cellular compartments by permeabilization of plasma and nuclear membranes in order to pinpoint differential compartment effects on variability in intracellular drug kinetics. These results provide the basis for a comprehensive model of the determinants of intracellular diffusion of small-molecule drugs, their target-seeking trajectories, and the consequences of these processes on the apparent kinetics of drug-target interactions.

Impaired nitric oxide (NO˙)-cyclic guanosine 3', 5'-monophosphate (cGMP) signaling has been observed in many cardiovascular disorders, including heart failure and pulmonary arterial hypertension. There are several enzymatic determinants of cGMP levels in this pathway, including soluble guanylyl cyclase (sGC) itself, the NO˙-activated form of sGC, and phosphodiesterase(s) (PDE). Therapies for some of these disorders with PDE inhibitors have been successful at increasing cGMP levels in both cardiac and vascular tissues. However, at the systems level, it is not clear whether perturbation of PDE alone, under oxidative stress, is the best approach for increasing cGMP levels as compared with perturbation of other potential pathway targets, either alone or in combination. Here, we develop a model-based approach to perturbing this pathway, focusing on single reactions, pairs of reactions, or trios of reactions as targets, then monitoring the theoretical effects of these interventions on cGMP levels. Single perturbations of all reaction steps within this pathway demonstrated that three reaction steps, including the oxidation of sGC, NO˙ dissociation from sGC, and cGMP degradation by PDE, exerted a dominant influence on cGMP accumulation relative to other reaction steps. Furthermore, among all possible single, paired, and triple perturbations of this pathway, the combined perturbations of these three reaction steps had the greatest impact on cGMP accumulation. These computational findings were confirmed in cell-based experiments. We conclude that a combined perturbation of the oxidatively-impaired NO˙-cGMP signaling pathway is a better approach to the restoration of cGMP levels as compared with corresponding individual perturbations. This approach may also yield improved therapeutic responses in other complex pharmacologically amenable pathways.