Original Research ARTICLE
In silico Therapeutics for Neurogenic Hypertension and Vasovagal Syncope
- 1Laboratory of Radiobiology and Molecular Genetics-080, Institute of Nuclear Sciences Vinča, University of Belgrade, Belgrade, Serbia
- 2Center for Multidisciplinary Research-180, Institute of Nuclear Sciences Vinča, University of Belgrade, Belgrade, Serbia
Neurocardiovascular diseases (NCVD) are the leading cause of death in the developed world and will remain so till 2020. In these diseases the pathologically changed nervous control of cardiovascular system has the central role. The actual NCV syndromes are neurogenic hypertension, representing the sympathetically mediated disorder, and vasovagal syncope, which is the vagally mediated disorders. Vasovagal syncope, the disease far from its etiological treatment, could benefit from recruiting and application of antimuscarinic drugs used in other parasympathetic disorders. The informational spectrum method (ISM), a method widely applied for the characterization of protein-protein interactions in the field of immunology, endocrinology and anti HIV drug discovery, was applied for the first time in the analysis of neurogenic hypertension and vasovagal syncope therapeutic targets. In silico analysis revealed the potential involvement of apelin in neurogenic hypertension. Applying the EIIP/ISM bioinformatics concept in investigation of drugs for treatment of vasovagal syncope suggests that 78% of tested antimuscarinic drugs could have anti vasovagal syncope effect. The presented results confirm that ISM is a promissing method for investigation of molecular mechanisms underlying pathophysiological proceses of NCV syndromes and discovery of therapeutics targets for their treatment.
NCVD are the syndromes where autonomic nervous system (Zoccoli et al., 2001; Bojic, 2003) dysfunction plays a dominant etiological role (Goldstain, 2001; Bojic et al., 2012a,b). NCV disorders can be classified as sympathetically mediated disorders (i.e., neurogenic hypertension, NH) vs. vagally mediated disorders (i.e., vasovagal syncope, VVS), though in many disorders both systems are dysfunctional (Goldstain, 2001).
NH is characterized by an increased level of sympathetic nervous activity (SNA) (Fisher and Paton, 2011). SNA is in reciprocal interaction with the number of important systems for the pathophysiological profile of NH, like inflammation, angiotension II system and vascular dysfunction. Recently, due to its potential causal role in the genesis of NH, an emphasis was put on the role of oxidative stress in brain stem structures.
Reactive oxygen spices (ROS) in the brain can be generated in angiotesin II dependent and angiotensin II independent manners. The most studied mechanism for investigating the potential causal treatment of NH is Ang 1-7-MAS receptor-NO mechanism (Zimmerman, 2011). This mechanism counteracts prohypertensive actions of ROS and represents a good choice for investigation of therapeutic candidate by Informational Spectrum Method (ISM).
In the case of VVS, Nucleus Tractus Solitarii (NTS) in the brain stem is stimulated either directly (central VVS) or indirectly (peripheral VVS), provoking an enhancement of vagal tone and withdrawal of SNA tone. This dual response causes a continuum of cardiovascular phenotype responses. At the “vagal end” it leads to the cardioinhibitory type of VVS (bradicardia as a dominant cause of the VVS), while on the “sympathetic end” causes the vasodepressor type of VVS (hypotension as a dominant cause of the VVS). Mixed type of VVS is on the midway between these two extremes.
The treatment of VVS involves a layered approach with a combination of lifestyle changes, physical maneuvers, medications, and implantable devices. The vast majority of patients with VVS can be adequately controlled with non-pharmacological approaches and do not require pharmacological treatment (Raj and Coffin, 2013). There is, however, the minority of patients with refractory and recurrent VVS who can benefit from effective pharmacotherapy.
Since the majority of the VVS patients suffer from mixed and cardioinhibitory type of VVS, we took under consideration the modulation of vagal tone as the in silico strategy for the new anti VVS drugs.
The therapy of NH and VVS is far from efficient etiological treatment. The development of new drugs for treatment of NH and VVS is time and money consuming process, which can be accelerated by in silico screening of the molecular libraries for candidate novel drugs and by repurposing of approved drugs.
The aim of this study was to apply EIIP/ISM approach to identify candidate neuropeptides and small-molecules for the potential therapeutics of NH and VVS. The study may lead to new insights in the field of neurocardiovascular pharmacotherapy and pathophysiology.
The Long Range Molecular Interactions
Current concepts of intermolecular interactions in biological systems are based on the surface complementarity between interacting biomolecules and assumption that the first contact between interacting molecules is achieved accidentally by the thermal motions that cause molecular wander. If proteins are considered as spheres of 18 Å radius (typical of a small protein), and if spheres associate with every contact, without regard to orientation, the diffusion-limited association rate constant, calculated according to Smoluchowski's equation (Smoluchowski, 1916) is 7 × 109 M−1s−1. However, before chemical bond formation takes place, reacting molecular regions must be positioned close enough (at a distance of ~2 Å) and the appropriate atoms must be held in the correct orientation for the reaction that is to follow, because the attractive forces involved in the recognition and binding of molecules include all the weak non-covalent forces. It means that the protein's binding site is only a small fraction (~0.1%) of the surface area. Taking into account this limitation, the diffusion-limited association rate constant, predicted from a three-dimensional (3D) “random diffusion” model and calculated according Smoluchowski's equation is ~106 M−1s−1 for a protein-ligand and ~103 M−1s−1 for a protein-protein interaction. Northrup and Erickson have noted that protein-protein association generally occurs at rates that are 103–104 times faster than would be expected from simple considerations of collision frequencies and strict orientation effects which assume that productive binding occurs only when the molecules collide within 2 Å of their final binding site (Northrup and Erickson, 1992).
Electron-Ion Interaction Potential (EIIP)
In order to overcome the discrepancy between theoretically estimated values and real values of the associated rate constant for a intermolecular interactions in biological systems, the long-range intermolecular interactions (distances between 5 and 1000 A) between interacting molecules was proposed (Veljkovic, 1980). It has been showed that the EIIP and the average quasivalence numbers (AQVN) Z* represent essential molecular descriptors which determines the long-range properties of biological molecules (Veljkovic, 1980). These two molecular descriptors are defined by the following equations:
i, Type of the chemical element;
Z, Valence of the ith chemical element;
n, Number of the ith chemical element atoms in the compound;
m, Number of types of chemical elements in the compond;
N, total number of atoms.
The EIIP values calculated according to the Equation (2) are in Rydbergs (Ry = 13.6 eV).
Informational Spectrum Method (ISM)
The ISM a virtual spectroscopy method for calculation of the long-range properties of biological macromolecules, is based on a model that assigns to each amino acid a defined parameter describing a physico-chemical property involved in the biological activity of the protein and corresponding to electron-ion interaction potential (EIIP) (Veljkovic et al., 1985).
ISM method consists in three basic steps:
1. Transformation of alphabetic code of primary protein structure into a sequence of numbers representing EIIP of each component.
2. Conversion of numerical sequence by fast Fourier Transformation into information spectrum, which reveals dominant freqency peaks of the whole organic molecule.
3. Consensus Information Spectrum (CIS) analysis between information spectrums of two potentially similar or interactive molecules, which reveales functional locus of the interaction of two molecules.
Peak frequencies in CIS are common frequency components for the analyzed sequences. A measure of similarity for each peak is the signal-to-noise ratio (S/N), the ratio between the signal intensity at one particular IS frequency and the main value of the whole spectrum.
Schematic Presentation of ISM is given in Figure 1.
Figure 1. Basic steps of ISM method. Sequence → Transformation of primary protein sequence into sequence of numbers, by assigning of EIIP to each amino acid → Numerical presentation → Discrete Fourier transformation → Informational spectrum (IS).
The strong connection between EIIP and Z* molecular descriptors of small molecules and their biological activities (carcinogenisity, antibiotic activity, antiviral activity, toxicity, etc.) has been documented (for review see Veljkovic et al., 2011 and references therein).
The method has been successful in identification of functional protein domains representing candidate therapeutic targets for anti-HIV drugs (Veljkovic et al., 2007), anthrax (Doliana et al., 2008), and human influenza viruses (Veljkovic et al., 2009a,b; Perovic et al., 2013). It is recognized and used in more than 100 research centers worldwide (Veljkovic et al., 2011).
Sequences and Molecular Formulas of Compounds
Sequences were retrieved from UniProt Database with following accession numbers: (a) sequences of primates MAS1: Homo sapiens P04201, Pan troglodytes H2QU00, Macaca mulatta F7GJU7, Pongo abelii H2PL76GN, Chlorocebus sabaeus A0A0D9RJ10, Gorilla gorilla gorilla G3R4L5, Papio anubis A0A096NZT0. (b) sequences of primates Angiotensin 1-7: Homo sapiens P01019, Pan troglodytes H2Q1B7, Macaca mulatta G7MFR4, Gorilla gorilla gorilla Q9GLP6. (C) sequences of primates Apelin: Homo sapiens Q9ULZ1, Macaca mulatta F7GX01, Gorilla gorilla gorilla G3S9L8, Chlorocebus sabaeus A0A0D9R7Q2. Sequences are shown in Data Sheet 1.
Molecular formulas of compounds: muscarinic modulators from US Patent 7786308, US Patent 7378447US, and Patent Application 20100311746A1.
The primary structure of proteins encodes the information represented by the informational spectrum (IS) frequencies that correspond to the protein biological function. Mutually interacting proteins share common information represented by peaks in their cross-spectrum (Veljkovic et al., 2008). IS of human MAS receptor (MAS1) is presented on Figure 2A. It contains characteristic peak at the frequency F(0.240).
Figure 2. Informational spectrum (IS) of (A) MAS1, (B) angiotensin 1-7, (C) apelin 28, and consensus informational-spectrum (CIS) of (D) MAS1 and angiotensin 1-7 and (E) MAS1 and apelin.
MAS1 sequences are highly homologous among primates [see Supplementary Material: Data Sheet 1 and Exported Multiple Sequence Alignment (Edgar, 2004) of MAS1 as artwork using Jalview (Waterhouse et al., 2009; Image 1)]. Informational spectrum (IS) of MAS1 of Macaca mulatta (Figure 3B) and consensus informational spectrum (CIS) of MAS1 of primates are shown in Figure 3A. It can be seen that information represented by the CIS frequency F(0.240), is evolutionary conserved among proteins having the same biological function, as it was shown previously (Glisic et al., 2008).
Figure 3. (A) Consensus informational spectrum (CIS) of primates MAS1; (B) Informational spectrum (IS) of MAS1 of Macaca mulatta; (C) CS of Macaca mulatta MAS1 and angiotensin 1-7 and (D) CS of Macaca mulatta MAS1 and apelin 28, (E) CS of Macaca mulatta MAS1 and apelin 31 (F) CS of Macaca mulatta MAS1 and apelin 36.
The sequences of peptides Angiotensin 1-7 and Apelins of other primates are the same as the sequences of human peptides (Data Sheet 1) and consequently spectra are the same.
Figure 2B represents IS of the of human Angiotensin 1-7 (Ang 1-7). By performing cross-spectral analysis (CS) of MAS receptor and Ang 1-7 we have identified that these two molecules share common information corresponding to the IS frequency F(0.057) (Figure 2D). To confirm the reliability of our findings we have performed CS analysis of Macaca mulatta MAS1 and Ang 1-7 and identified the same frequency F(0.057) as dominant and evolutionary conserved (Figure 3C).
We further apply ISM analyses to identify new peptide interactors of human MAS1 representing potential candidate therapeutic agents. To identify peptides which share the common information represented by the frequency component F(0.057), 270 human peptides of the Human Neuropeptide sequence database (Kim et al., 2011) were screened by ISM and peptide apelin 28 (Figure 2C) was found to be among the peptides with the highest amplitudes and Signal/Noise values at the frequency F(0.057) in CIS with MAS1 (Figure 2E). Similar results were obtained for apelin-31 and apelin-36. In further analysis it was shown that the information represented by CS frequency for Macaca mulatta MAS receptor and apelins is the same (Figures 3D–F). According to the IS criterion these peptides are the potential candidate interactors of MAS1. Presented results indicate apelin as potential modulator of MAS1 receptor and novel candidate for treatment of NH.
With respect to the muscarinic antagonists applied in the therapy of VVS (atropine, propanteline bromide, and scopolamine), on the basis of AQVN/EIIP analysis we analyzed the distribution of other muscarinic antagonists which are not tested for their efficacy in the treatment of VVS (Figures 4, 5). The analysis revealed that 78% of tested drugs are in the active AQVN/EIIP domain, implying that there is consistent number of known antimuscarinic drugs that might have a therapeutic impact on VVS.
Figure 4. Schematic presentation of the EIIP/AQVN criterion for selection of candidate antimuscarinic drugs with potential therapeutic effect on VVS. Active domain (red): AQVN (2.25–2.65), EIIP (0.074–0.096). Chemical space (blue) AQVN (2.40–3.30) EIIP (0.000–0.116) EIIP/AQVN domain of homologous distribution of >90% compounds from PubChem Compound Database. Statistics: inside the active domain: 21 (78%), outside the active domain: 6 (22%). Asterisk (*)-position of atropine and propentheline bromide, known anti VVS drugs, within AQVN/EIIP active domain.
Figure 5. The list of antimuscarinic drugs tested for their AQVN/EIIP. Shaded fields- antimuscarinic drugs within the active domain (Figure 3), white fields-antimuscarinic drugs out of the active domain, bordered fields- antimuscarinic drugs tested for their therapeutical effect on VVS. Probability for the anti VVS activity of muscarinicdrugs, assessed by EIIP/AQVN criterion, is represented with the gray scale (black-the highest probability, white-the lowest probability).
We also analyzed muscarinic receptor modulators from three patents which are randomly selected among more than 30,000 patents from the patent database (http://www.freepatentsonline.com), in order test the suitability of the proposed EIIP/AQVN filter. Results of this analysis are presented in Table 1. As can be seen, all patented compounds are located within the active EIIP/AQVN domain (Figure 4).
Table 1. EIIP/AQVN values of patented muscarinic modulators: Muscarinic modulators (US Patent 7786308), Muscarinic agonists (US Patent 7378447), and Modulators of muscarinic receptors (US Patent Application 20100311746A1).
Brain oxidative stress/status defines the state of sympathoexcitation (increased level of ROS) or sympathoinhibition (the effect of NO), via angiotensinergic neurons (Zimmerman, 2011). Since oxidative injury contributes also to the progressive development of NH, targeting oxidative stress may represent one of promising therapeutic strategies for NH treatment.
The ISM analysis identified three forms of apelin as the potential candidates-interactors of MAS receptor on neurons in brainstem sympathetic centers and indicated apelin as potential modulator of MAS1 receptor and novel candidate for treatment of NH (Figure 6). There is lot of recent evidence of the role of apelin as therapeutic agent in protection of the brain against ischemic/reperfusion injury (Yang et al., 2015b), promotion of neurological function recovery after ischemic brain injury (Gu et al., 2013) and potential to cure acute and chronic neurological diseases (Cheng et al., 2012). Synthetic and biased agonists of apelin have been also developed and latter have shown in proof-of-concept studies clinical potential (Yang et al., 2015a). Apelin is adipokine secreted by peripheral tissues but also present in the hypothalamic neurons, having a major impact on the genesis and progression of diabetes (Drougard et al., 2014). This impact is obtained through ROS impact on sympathetic neuronal centers in the diencephalon which cause sympathoexcitation and consequently, the liver glycogenolysis and glusoneogenesis. Our data for the first time point on the potential role of apelin in the development of NH, again, through the action of ROS on sympathetic neuronal centers. Differently from other studies that investigate majorly apelin-APJ receptor signaling pathway (Chun et al., 2008; Yu et al., 2014) our analysis points on the novel mechanism, the apelin-MAS receptor dependent mechanism that could counterbalance sympathostimulating action of ROS. On the basis of the data in the literature (Chun et al., 2008; Yu et al., 2014) and on the basis of our results, there is high probability that apelin has sympathoinhibiting, anti ROS action also on the brainstem level. Further in vitro and in vivo studies that would elucidate the functional significance of apelin-MAS receptor interaction are needed. This finding could be of importance in designing a novel therapeutics for NH, but also for better understanding and treatment of the states like metabolic syndrome, where diabetes, hypertension, obesity, and obstructive sleep apnea occur in a cluster (Katsiki et al., 2014).
Figure 6. Proposed mechanism of antagonizing effect of apeline-MAS receptor and Angiotensin II-ATR1 molecular signaling pathways on preganglionic sympathetic neurons in the brainstem (adapted from Northrup and Erickson 1992.) Ang II, angiotensin II; ATR1, angiotensin 1 receptor; MAS R, MAS receptor; O2, superoxide; NO, nitric oxide; nNOS, neural nitric oxide synthase; SNA, sympathetic nervous activity; MAP, mean arterial pressure.
Muscarinic antagonists that were already tested for their efficacy in the treatment of VVS were:
1. atropine(i.v.), shown to be fully effective in the cardio-inhibitory form of tilt-induced vasovagal reflex, but with limited action in the vasodepressor form (Santini et al., 1999),
2. propantheline bromide (p.o.), shown to be highly effective in preventing VVS. In addition, propantheline bromide's effectiveness is present in the vasodepressor VVS and supports a role of direct cholinergic control of vascular tone (Yu and Sung, 1997),
3. scopolamine (t.d.), which showed no effect on VVS (Lee et al., 1996).
An AQVN/EIIP approach for new candidates in pharmacotherapy of VVS reveals that on the basis of their molecular properties, the majority of antimuscarinic drugs (78% of investigated substances) might have therapeutical potential for VVS (Figures 4, 5). More, the remaining 22% were out of AQVN/EIIP active domain, together with scopolamine, the drug that failed to show anti VVS therapeutical effect (Figure 4). Best candidates, identified by AQVN/EIIP approach, for future clinical studies are hydroxysine and cyclopentoate (Figure 5). The suitability of our method was further confirmed by the finding that all muscarinic modulators patented in three patents randomly selected among more than 30,000 patents from the patent database (http://www.freepatentsonline.com) have AQVN/EIIP values within AQVN/EIIP active domain (Table 1, Figure 4).
Computer-aided design techniques based on the long-range intermolecular interactions offer an insight into therapeutical and pathophysiological aspects of NH and VVS. These bioinformatics approaches also open a promising path for molecular investigation of other NCVD and could be of help in design of in vitro and in vivo animal and clinical studies of these diseases.
TB, VP and SG contributed equally to the conception of the work, acquisition, analysis and interpretation of data. TB, VP, and SG participated in drafting the manuscript, revisiting it critically and gave final approval of the version to be published. The authors reached the agreement to be accountable for all the aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was financed by the Ministry of Education, Science and Technological Development of the Republic of Serbia, project III 41028 and 173001.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnins.2015.00520
Data Sheet 1. Sequences of (a) primates MAS1 (b) sequences of primates Angiotensin 1-7 (c) sequences of primates Apelin 28, Apelin 31, and Apelin 36.
Bojic, T. (2003). Mechanisms of Cardiovascular Control and Effects of Acoustic Stimulation on Cardiovascular System during the Wake-Sleep Cycle. Ph.D. thesis, Alma Mater Università di Bologna, Bologna.
Bojic, T., Radak, D., Putnikovic, B., Alavantic, D., and Isenovic, E. R. (2012a). Methodology of monitoring cardiovascular regulation: baroreflex and central mechanisms of cardiovascular regulation. Vojnosanit. Pregl. 69, 1084–1090. doi: 10.2298/VSP110707019B
Bojic, T., Sudar, E., Mikhailidis, D. P., Alavantic, D., and Isenovic, E. R. (2012b). The role of G protein coupled receptor kinases in neurocardiovascular pathophysiology. Arch. Med. Sci. 8, 970–977. doi: 10.5114/aoms.2012.32466
Chun, H. J., Ali, Z. A., Kojima, Y., Kundu, R. K., Sheikh, A. Y., Agrawal, R., et al. (2008). Apelin signaling antagonizes Ang II effects in mouse models of atherosclerosis. J. Clin. Invest. 118, 3343–3354. doi: 10.1172/jci34871
Doliana, R., Veljkovic, V., Prljic, J., Veljkovic, N., De Lorenzo, E., Mongiat, M., et al. (2008). EMILINs interact with anthrax protective antigen and inhibit toxin action in vitro. Matrix Biol. 27, 96–106. doi: 10.1016/j.matbio.2007.09.008
Drougard, A., Duparc, T., Brenachot, X., Carneiro, L., Gouazé, A., Fournel, A., et al. (2014). Hypothalamic apelin/reactive oxygen species signaling controls hepatic glucose metabolism in the onset of diabetes. Antioxid. Redox Signal. 20, 557–573. doi: 10.1089/ars.2013.5182
Glisic, S., Arrigo, P., Alavantic, D., Perovic, V., Prljic, J., and Veljkovic, N. (2008). Lipoprotein lipase: a bioinformatics criterion for assessment of mutations as a risk factor for cardiovascular disease. Proteins 70, 855–862. doi: 10.1002/prot.21581
Gu, Q., Zhai, L., Feng, X., Chen, J., Miao, Z., Ren, L., et al. (2013). Apelin-36, a potent peptide, protects against ischemic brain injury by activating the PI3K/Akt pathway. Neurochem. Int. 63, 535–540. doi: 10.1016/j.neuint.2013.09.017
Katsiki, N., Athyros, V. G., Karagiannis, A., and Mikhailidis, D. P. (2014). Metabolic syndrome and non-cardiac vascular diseases: an update from human studies. Curr. Pharm. Des. 20, 4944–4952. doi: 10.2174/1381612819666131206100750
Northrup, S. H., and Erickson, H. P. (1992). Kinetics of protein-protein association explained by Brownian dynamics computer simulation. Proc. Natl. Acad. Sci. U.S.A. 89, 3338–3342. doi: 10.1073/pnas.89.8.3338
Perovic, V. R., Muller, C. P., Niman, H. L., Veljkovic, N., Dietrich, U., Tosic, D. D., et al. (2013). Novel phylogenetic algorithm to monitor human tropism in Egyptian H5N1-HPAIV reveals evolution toward efficient human-to-human transmission. PLoS ONE 8:e61572. doi: 10.1371/journal.pone.0061572
Raj, S. R., and Coffin, S. T. (2013). Medical therapy and physical maneuvers in the treatment of the vasovagal syncope and orthostatic hypotension. Prog. Cardiovasc. Dis. 55, 425–433. doi: 10.1016/j.pcad.2012.11.004
Santini, M., Ammirati, F., Colivicchi, F., Gentilucci, G., and Guido, V. (1999). The effect of atropine in vasovagal syncope induced by head-up tilt testing. Eur. Heart. J. 20, 1745–1751. doi: 10.1053/euhj.1999.1697
Veljkovic, N., Glisic, S., Perovic, V., and Veljkovic, V. (2011). The role of long-range intermolecular interactions in discovery of new drugs. Expert Opin. Drug. Discov. 6, 1263–1270. doi: 10.1517/17460441.2012.638280
Veljkovic, N., Glisic, S., Prljic, J., Perovic, V., Botta, M., and Veljkovic, V. (2008). Discovery of new terapeutic targets by the informational spectrum method. Curr. Protein Pept. Sci. 9, 493–562. doi: 10.2174/138920308785915245
Veljkovic, V., Cosic, I., Dimitrijevic, B., and Lalovic, D. (1985). Is it possible to analyze DNA and protein sequences by the methods of digital signal processing? IEEE Trans. Biomed. Eng. 32, 337–341.
Veljkovic, V., Niman, H. L., Glisic, S., Veljkovic, N., Perovic, V., and Muller, C. P. (2009a). Identification of hemagglutinin structural domain and polymorphisms which may modulate swine H1N1 interactions with human receptor. BMC Struct. Biol. 9:62. doi: 10.1186/1472-6807-9-62
Veljkovic, V., Veljkovic, N., Esté, J. A., Huther, A., and Dietrich, U. (2007). Application of the EIIP/ISM bioinformatics concept in development of new drugs. Curr. Med. Chem. 14, 441–453. doi: 10.2174/092986707779941014
Veljkovic, V., Veljkovic, N., Muller, C. P., Müller, S., Glisic, S., Perovic, V., et al. (2009b). Characterization of conserved properties of hemagglutinin of H5N1 and human influenza viruses: possible consequences for therapy and infection control. BMC Struct. Biol. 9:21. doi: 10.1186/1472-6807-9-21
Waterhouse, A. M., Procter, J. B., Martin, D. M. A., Clamp, M., and Barton, G. J. (2009). Jalview Version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics 25, 1189–1191. doi: 10.1093/bioinformatics/btp033
Yang, P., Maguire, J. J., and Davenport, A. P. (2015a). Apelin, Elabela/Toddler, and biased agonists as novel therapeutic agents in the cardiovascular system. Trends Pharmacol. Sci. 36, 560–567. doi: 10.1016/j.tips.2015.06.002
Yu, X. H., Tang, Z. B., Liu, L. J., Qian, H., Tang, S. L., Zhang, D. W., et al. (2014). Apelin and its receptor APJ in cardiovascular diseases. Clin. Chim. Acta 428, 1–8. doi: 10.1016/j.cca.2013.09.001
Keywords: in silico analysis, neurocardiovascular diseases, neurogenic hypertension, protein-protein interaction, vasovagal syncope
Citation: Bojić T, Perović VR and Glišić S (2016) In silico Therapeutics for Neurogenic Hypertension and Vasovagal Syncope. Front. Neurosci. 9:520. doi: 10.3389/fnins.2015.00520
Received: 03 September 2015; Accepted: 24 December 2015;
Published: 21 January 2016.
Edited by:James J. Galligan, Michigan State University, USA
Reviewed by:Lu Liu, University of New South Wales, Australia
Norbert Ikechukwu Nwankwo, University of Port Harcourt, Nigeria
Copyright © 2016 Bojić, Perović and Glišić. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.