Next Article in Journal
Design and Synthesis of Coumarin Derivatives as Cytotoxic Agents through PI3K/AKT Signaling Pathway Inhibition in HL60 and HepG2 Cancer Cells
Next Article in Special Issue
Antibodies-Abzymes with Antioxidant Activities in Two Th and 2D2 Experimental Autoimmune Encephalomyelitis Mice during the Development of EAE Pathology
Previous Article in Journal
Diaza-1,3-butadienes as Useful Intermediate in Heterocycles Synthesis
Previous Article in Special Issue
Effects of Chemical Structures Interacting with Amine Oxidases on Glucose, Lipid and Hydrogen Peroxide Handling by Human Adipocytes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Insights into β-Carboline Inhibition of Monoamine Oxidase A

1
National Institute of Chemistry, SI-1000 Ljubljana, Slovenia
2
Jožef Stefan Institute, SI-1000 Ljubljana, Slovenia
3
Faculty of Pharmacy, University of Ljubljana, SI-1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(19), 6711; https://doi.org/10.3390/molecules27196711
Submission received: 16 September 2022 / Revised: 3 October 2022 / Accepted: 6 October 2022 / Published: 9 October 2022

Abstract

:
Monoamine oxidases (MAOs) are an important group of enzymes involved in the degradation of neurotransmitters and their imbalanced mode of action may lead to the development of various neuropsychiatric or neurodegenerative disorders. In this work, we report the results of an in-depth computational study in which we performed a static and a dynamic analysis of a series of substituted β-carboline natural products, found mainly in roasted coffee and tobacco smoke, that bind to the active site of the MAO-A isoform. By applying molecular docking in conjunction with structure-based pharmacophores and molecular dynamics simulations coupled with dynamic pharmacophores, we extensively investigated the geometric aspects of MAO-A binding. To gain insight into the energetics of binding, we used the linear interaction energy (LIE) method and determined the key anchors that allow productive β-carboline binding to MAO-A. The results presented herein could be applied in the rational structure-based design and optimization of β-carbolines towards preclinical candidates that would target the MAO-A enzyme and would be applicable especially in the treatment of mental disorders such as depression.

1. Introduction

Monoamine oxidases (MAOs) are an important group of enzymes responsible for the oxidative deamination of amine neurotransmitters, primarily in the central nervous system, thereby regulating their levels. When the function of MAOs is increased, this can lead to the development of mental disorders (e.g., autism and depression), and even their normal function is associated with the production of hydrogen peroxide, which can lead to neuronal damage and the development of neurodegenerative disorders (Parkinson’s or Alzheimer’s disease).
MAOs exist in two distinct isoforms, MAO-A and MAO-B, which share about 70% homology in their sequence but differ in some parts of their structure, tissue distribution, and selectivity—MAO-A preferentially metabolizes serotonin and norepinephrine, whereas MAO-B predominantly degrades dopamine [1]. Both isoenzymes contain in their active site the FAD (flavin adenine dinucleotide) co-factor bound to one of the cysteine residues. The overall catalytic cycle of MAOs begins with the conversion of the amine substrate to the corresponding imine, which is followed by the regeneration of the reduced flavin cofactor (FADH2) with molecular oxygen. The latter is the part of the catalytic cycle in which hydrogen peroxide is formed. The rate-limiting step of this cycle is the cleavage of the C-H bond, vicinal to the amino group and the subsequent transfer of the hydrogen atom to the flavin cofactor. There are several proposals for the mechanism of the rate-limiting step [2,3,4,5,6,7,8], with most evidence pointing to the hydride transfer mechanism [9,10,11,12,13].
Because of their involvement in the development of neurodegenerative and neuropsychiatric diseases, MAOs are intriguing targets for drug development, several of which are already on the market. The best-known drugs are the MAO-B-specific rasagiline and selegiline (both used mainly to treat Parkinson’s disease) and the MAO-A-specific moclobemide (used mainly to treat depression). In addition, studies have shown that smoking can also affect the activity of MAO. Results interestingly suggest that tobacco smokers have up to 28% lower brain MAO-A activity and 40% lower brain MAO-B activity compared with nonsmokers [14,15]. Considering the involvement of MAOs in the development of neurodegenerative diseases, this may be related to the lowest rate of Parkinson’s disease in smokers [16]. Several studies suggest that β-carbolines which are found primarily in tobacco smoke, coffee, dark chocolate, salmon, raisins, certain spices, cooked foods, and alcoholic beverages, among others [17], are responsible for the inhibition of MAO [18].
β-carbolines (or harmala alkaloids) are naturally occurring, biologically active alkaloids that are derivatives of indole [19] and whose basic structural element is shown in Figure 1. They were first isolated from the plant Peganum harmala, a succulent in the Nitrariaceae family [19], but are also found in several other plant families [20]. Some of the harmala alkaloids are the reason for the hallucinogenic properties of ayahuasca, a psychoactive drink used in spiritual ceremonies in South America and Africa [21,22]. They not only inhibit MAO enzymes but are also known to produce a variety of different effects in the body, with their mechanism of action being associated with imidazole [23,24], serotonin [25], benzodiazepine [26] and dopamine receptors [27,28], as well as affecting cerebral neurotransmitter concentrations [29]. They have also been shown to intercalate into DNA [30,31,32], interact with DNA topoisomerases [31], and have a protective role against oxidative stress in human tissues when accumulated [33]. Other studies focused on their antiplatelet [34], antimalarial [35] and antidepressant effects [36]. The latter, in particular, was the main topic of many studies and some specific β-carbolines have been analyzed in detail. Harman has been shown to reduce depressive behavior in animals, presumably through several different modes of action [37,38,39], the most important of which is MAO-A inhibition [40]. Similarly, preclinical studies have shown that harmine has a potential antidepressant effect in various depression models [39,41,42]. Currently WHO estimates that depression affects about 280 million people worldwide. This gives rise to the urgent need to develop new antidepressants that can effectively treat depression and improve the quality of life of those affected [43].
We can divide the β-carbolines into three groups based on the saturation of one of their rings: fully aromatic β-carbolines, dihydro-β-carbolines which are partially saturated, and fully saturated tetrahydro-β-carbolines [44]. All usually contain different substituents on one of their three rings. Depending on the surrounding solvent and pH, the compounds can exist in one of four states: as a cation, in neutral form, as a zwitterion, or as an anion [45]. Studies focusing on the selectivity of some β-carbolines have found that a methyl group on the C1 atom appears to be critical for the compounds to be selective towards MAO-A rather than MAO-B—norharman (without substituents) is 10–20 times less active as a MAO-A inhibitor, but two orders of magnitude more active as a MAO-B inhibitor than harman or harmaline (both have a methyl group at position 1) [46]. A crystal structure of one of the fully aromatic β-carbolines (harmine, methyl group at position 1 and a methoxy group at position 7) bound to the active site of MAO-A was determined [47] (Figure 1) and harmine was characterized as a reversible MAO-A inhibitor. This crystal structure presents an excellent starting point for our computational studies.
The aim of the present work was to study the binding properties and molecular recognition between β-carbolines and the MAO-A isoform at the atomistic level by applying a variety of computational approaches. The study is based on a series of differently substituted β-carbolines [48], and we sought to first gain a static insight into their binding and later extended this with a more dynamic approach. Within the latter, we performed molecular dynamics (MD) simulations of the considered β-carbolines and analyzed their geometry (RMSD), interactions with the MAO-A active site (dynophores) and calculated their predicted binding free energy (with the linear interaction energy method). All the data obtained could be used to develop new MAO-A selective inhibitors that could potentially be used in the effective treatment of mental disorders, especially depression.

2. Results

2.1. Static Insight into the β-Carboline–MAO-A Molecular Recognition

We began our computational analysis of a selected series of β-carbolines 1–11, which bind non-covalently into the active site of the MAO-A isoform, with a static view of our system. A crystal structure of MAO-A with harmine (compound 6), a β-carboline that acts as a reversible MAO-A inhibitor, had already been determined (PDB code 2Z5X), providing an excellent starting point for our study [47]. Thus, we were able to first redock harmine into the MAO-A active site and deduce whether our docking software is able to reproduce the observed X-ray binding pose. For this purpose, we used three different scoring functions, available in GOLD [49], namely: GoldScore, ChemScore and ChemPLP. All three scoring functions were able to successfully reproduce the X-ray binding pose (see Supplementary Materials, Figure S14). Consequently, we were able to proceed with docking the remainder of the compounds studied.
Using the LigandScout program package [50], we then visualized the key intermolecular interactions by generating a 3D structure-based pharmacophore for harmine (compound 6, Figure 2A), 6-methoxyharmalan (compound 10, Figure 2B), and other β-carbolines based on the binding modes determined using the GoldScore scoring function. In the X-ray-determined binding model of harmine in the MAO-A active site, researchers previously concluded that harmine interacts with Tyr69, Asn181, Phe208, Val210, Gln215, Cys323, Ile325, Ile335, Leu337, Phe352, Tyr407, Tyr444, and the cofactor flavin adenine dinucleotide (FAD), with seven water molecules occupying the space between harmine and these amino acids [47]. Our pharmacophore model detects interactions of harmine with six of these groups (Tyr69, Phe208, Phe352, Ile335, Tyr407, and Tyr444) as well as an additional interaction with Ile180 and hydrogen bonds with the conserved water molecules (Figure 2A).
We can roughly divide the studied set of β-carbolines into two substructural groups: those that are neutral (compounds 1, 3, 6, 9, and 10) and those that have a positive charge on the pyridine nitrogen atom (2, 4, 5, 7, 8, and 11), with some important differences in the pharmacophores of the two groups. All neutral compounds (except inhibitor 9) exhibit hydrophobic interactions of all three aromatic β-carboline rings with MAO-A active site residues, namely tyrosine (Tyr407), isoleucines (Ile180, Ile335), leucine (Leu337), phenylalanines (Phe208, Phe352), and methionine (Met350). The addition of a methyl group to the inhibitor structure leads to a further hydrophobic interaction with the active site tyrosines (Tyr407 and Tyr444). In all neutral β-carbolines, the pyridine nitrogen atom acts as a hydrogen bond acceptor (interaction with tyrosine residues), while the pyrrole nitrogen acts as a hydrogen bond donor (interaction with the conserved water molecules). The pharmacophore models predict no additional interactions as a result of the addition of the methoxy group (compounds 6, 9, and 10).
For the β-carbolines with a positive charge on the pyridine nitrogen atom (compounds 2, 4, 5, 7, 8, and 11), pharmacophore models generally predict fewer interactions than for the neutral β-carbolines. The focus in this group is on the hydrophobic interactions formed between the active site residues and the benzene ring portion of the β-carboline structure. In contrast to what we found for the neutral β-carboline group, in this case, we see interactions occurring with the oxygen atom of the methoxy group. The latter acts as a hydrogen bond acceptor with water or with one of the MAO-A active site tyrosines. In addition, for compounds 8 and 11, the model predicts an interaction between the positively charged nitrogen atom and a phenylalanine residue.
Furthermore, we wanted to determine if there was a correlation between the experimentally determined binding free energies of the β-carbolines and the values predicted by all three scoring functions available in GOLD: GoldScore, ChemScore, and ChemPLP. That is, we wanted to determine whether a more negative free energy of binding (i.e., a smaller inhibition constant Ki, meaning that the compound is a better inhibitor) corresponds to a higher scoring function value and vice versa. Our results show that there is little to no correlation between the values of the scoring functions and the experimentally determined binding free energy (see Figure 2C). This is not surprising since docking is usually a good predictor of the binding geometry but generally fails in exactly predicting the binding energy [51]. To obtain this information, we need to use more thorough dynamical methods [52], so we have extended our study to include MD simulations.

2.2. Binding Free Energies of β-Carbolines

Following the established linear interaction energy (LIE) protocol, we first calculated the interaction energies between compounds 1–11 in the aqueous environment (free state, see Figure 3A, left) and in the protein active site (bound state, see Figure 3A, right). For all compounds considered, the interaction energy in the bound state is more favorable than in the free state (by about 3–10 kcal mol−1, depending on the compound), which is in line with our expectations and indicates that the binding of the compounds is energetically favored. Interestingly, splitting the interaction energy into the van der Waals and electrostatic components shows that the bound state is energetically more favorable with respect to the van der Waals component, while the opposite is true for the electrostatic component. This is not surprising since the active site of MAO-A is highly hydrophobic (as confirmed by our additional detailed active site analysis), so it is to be expected that van der Waals interactions are the predominant driving force for binding these β-carbolines. The compounds may form more electrostatic interactions in water than in the very hydrophobic active site.
In general, the observed fluctuations (i.e., standard deviations) for the electrostatic contributions of charged compounds 1, 3, 6, 9, and 10 (average standard deviation of 1.9 kcal mol−1) are higher than for neutral compounds 2, 4, 5, 7, 8, and 11 (average standard deviation of 1.0 kcal mol−1), which is consistent with our expectations because the electrostatics of charged moieties are expected to fluctuate more than those of their neutral counterparts. The fluctuations for the van der Waals contributions are generally smaller, with no significant difference between the charged and uncharged β-carbolines (0.7 and 0.6 kcal mol−1, respectively) (Table 1).
Figure 3B shows the values of the van der Waals and electrostatic components of the interaction energy for β-carboline 2 over the course of the MD simulation in the bound and free states. The van der Waals component of the interaction energy varies between −10 and −25 kcal mol−1 in the free state and between −30 and −40 kcal mol−1 in the bound state. There is a clear separation between the two states, with a clear preference for the bound state. On the other hand, the electrostatic component varies between −70 and −100 kcal mol−1 in the bound state and between −70 and −130 kcal mol−1 in the free state, with a clear overlap between the two states. In this case, the fluctuations are slightly larger than for the van der Waals component, and we can infer a slight preference for the free state, but the distinction between the two states is much less pronounced.
The experimentally determined inhibition constants (Ki) for the examined compounds were taken from ref. [48]. The binding free energy was calculated from these values by using the standard formula:
G b i n d i n g = R T ln K i
where R is the gas constant and T is the temperature. The binding free energy was calculated according to the LIE formalism, as described in the computational methods section. The values determined for the coefficients α, β and γ were 0.54, 0.26 and −2.65, respectively, resulting in the final LIE equation:
G b i n d i n g = 0.54 ( V L P v d W V L W v d W ) + 0.26 ( V L P e l V L W e l ) 2.65
Coefficient α is, therefore, parameterized to be larger than coefficient β, giving greater weight to the van der Waals component of the interaction energy, which in turn is consistent with the hydrophobicity of the active site and the compounds themselves. It should be noted that similar α and β coefficients have been frequently determined in LIE studies of various bio-macromolecular systems [53,54,55]. The experimentally derived and the calculated binding free energies for all investigated β-carbolines are listed in Table 2, along with the average RMSD values of the ligand for the duration of the MD simulation.
The binding free energy values calculated from our MD simulations are in broad agreement with the available experimental data. Figure 3C further illustrates the correlation between the experimentally determined and the calculated free energies of binding together with the standard deviations. Looking at the deviations, most of the values fall on the diagonal line, which is a perfect match between the experimental and calculated free energies, with some inhibitors deviating to a larger extent such as compounds 2, 6, 8, and 10.

2.3. Dynamic Insight into the β-Carboline–MAO-A Molecular Recognition

The obtained MD trajectories of all ligands 1–11 in the bound state provide a considerable amount of structural data for further analysis to allow a more dynamic picture of the molecular recognition between β-carbolines and the active site of the MAO-A enzyme. To complement the binding free energy data obtained by LIE calculations, we investigated the interaction between β-carbolines and MAO-A using dynamic pharmacophore models. The dynamic pharmacophore (dynophore) approach allows us to combine the features of MD simulations and static 3D structure-based pharmacophore models, as it provides us with the statistical characterization of the pharmacophore features of the ligand throughout the simulation. This information is then distilled into so-called “superfeatures” represented by the clouds of a given pharmacophore model [56,57,58]. Furthermore, the dynamic picture of molecular recognition was complemented by RMSD calculations to evaluate the stability of the proposed ligand-binding pose.
Compound 1 (i.e., norharman) is the simplest of the β-carbolines studied, having no additional methyl or methoxy groups. It is also the least potent inhibitor, as indicated by the experimentally determined inhibition constant (Ki = 3.34 ± 0.10 µM), which is converted to a binding free energy value of −7.6 kcal mol−1. This was well reflected in our LIE calculations, with compound 1 having the lowest calculated binding free energy of −7.9 ± 0.5 kcal mol−1. This observation could be related to the highest RMSD value in the group of simulated β-carbolines, i.e., compound 1 undergoes a significant geometric shift at about the 30 ns mark of the simulation (see Figure 4C). Before the shift, this compound is seen in a binding pose similar to the X-ray-determined pose of β-carboline harmine (compound 6), whereas after the shift, it occupies a different part of the active site. It no longer interacts with the aromatic cage tyrosines (Tyr407 and Tyr444), which have been shown to be very important for binding [59]. This pose is, therefore, not so well-suited for MAO-A inhibition, which is clearly reflected in our results, especially in the dynophore analysis (Figure 4A,B). The latter shows hydrophobic interactions between the aromatic rings of inhibitor 1 and several active site amino acids (Ile180, Ile335, Leu337, Met350, Phe352, Tyr407, and Tyr444) as well as the pyridine and pyrrole nitrogen atoms, which act as hydrogen bond acceptor and donor, respectively. Frame-by-frame analysis clearly shows that the interactions with both tyrosines are present in the first 30 ns of the simulation. After the molecule moves away from the aromatic cage, new interactions appear with Leu337, Met350, and Ile335.
β-carbolines 2–5 represent variations in compound 1 with additional methyl groups at different positions. In all compounds, most of the interactions are hydrophobic and originate mainly from the two aromatic rings. In compounds 2, 4, and 5, there is an additional methyl group on the pyridine nitrogen, which makes it positively charged and unable to act as a hydrogen bond acceptor, as is the case in compounds 1 and 3. An additional methyl group on the C1 atom (next to the pyridine nitrogen, see Figure 1) in compounds 3, 4, and 5 is able to form additional stabilizing hydrophobic interactions with the active site tyrosines (Tyr197, Tyr407, andTyr444) and isoleucines (Ile180 and Ile207). Interestingly, in compound 3, the pyrrole nitrogen can act as a hydrogen bond donor and interact with Asn181.
Figure 4D shows the distance between the pyrrole nitrogen and Asn181 for compounds 1–5 during the MD simulation. For inhibitors 2 and 3, this distance is relatively stable, with the average distance for compound 3 being slightly shorter, presumably resulting in the interaction being reflected in the dynophore analysis. For the other three inhibitors, the distance is less stable and for inhibitor 1, we can see a clear correlation between the distance and the interaction pattern in the dynophore analysis. Namely, the dynophore analysis reflects the interaction between Asn181 and the pyrrole nitrogen when the distance is shorter than about 2.9 Å.
All four compounds in this group have relatively small RMSD values (see Table 2), implying that their binding modes are stable throughout the MD simulation. Moreover, the calculated LIE values for the free energy of binding of inhibitors 3–5 agree remarkably well with the experimentally determined values for their free energies of binding.
According to experimental data, β-carboline 6 (i.e., harmine) is the strongest MAO-A inhibitor in the series (Ki = 0.005 ± 0.0002 µM), and the binding free energy value we calculated also places it among the strongest inhibitors (ΔGcalc = −10.0 ± 0.4). Structurally, it is most similar to inhibitor 3, with an additional methoxy group present at the C7 atom. Therefore, the predicted interactions (see Figure 5A) are similar to those of compound 3: hydrophobic interactions of the two aromatic rings and the methyl group, and the pyridine and pyrrole nitrogen as hydrogen bond acceptor and donor, respectively. Additional interactions are possible with the methoxy group.
Dynophore analysis predicts interactions with Leu337, Ile180, Ile335, Phe208, Val210, Tyr407, Tyr444, Phe352, Ser209, Cys323, Tyr197, Asn181, and Ile207. This covers 9 of the 12 interactions determined in the X-ray structure (compared to 6 of 12 that we determined using static pharmacophore analysis). The RMSD for compound 6 is relatively large, and (similar to inhibitor 1) we can see a clear shift in the molecule at about the 40 ns mark of the simulation (Figure 5C). In contrast to compound 1, the molecule moves closer to the aromatic tyrosine cage. Figure 5B clearly shows the appearance of interactions with Tyr407 and Tyr444 at this time point and the disappearance (or weakening) of interactions with Leu337, Ile180, and Phe352. At the same time point, we also note an abrupt shortening of the distance between the pyrrole nitrogen and residue Asn181 and a subsequent stabilization of this distance (Figure 5D).
Similar to inhibitor 6, compounds 7–9 also have a methoxy group at the C7 atom, with additional methyl groups at various positions. The dynophore analysis predicts hydrophobic interactions of aromatic rings and the methyl group at the C1 atom, as well as additional interactions of the methoxy oxygen atom with Cys323 and water molecules. Analogously to compound 6, some interactions with the pyrrole or pyridine nitrogen atom can also be formed. The calculated RMSD values for inhibitors 7–9 are relatively small, while the agreement between the experimentally determined free energies of binding and their calculated counterparts is not optimal.
Finally, compounds 10 and 11 both have a methoxy group on the C6 atom instead of the C7 atom. The dynophore analysis of both inhibitors (Figure 6A,C) again shows mainly the presence of hydrophobic interactions with both rings and the methyl group at C1. The prevalence of hydrophobic interactions with the pyridine ring is significantly lower for inhibitor 10, probably due to there being one less double bond. The average calculated RMSD value is quite high for both inhibitors (3.7 Å and 3.5 Å, respectively). Compound 10 clearly moves away from the aromatic cage at 10–20 ns of the MD simulation (Figure 6B). On the other hand, compound 11 undergoes a larger geometric shift in the first 5 ns of the simulation and remains relatively stable thereafter (Figure 6D). The calculated binding free energies for compounds 10 and 11 are −10.1 and −8.5 kcal mol−1, respectively, the latter being in better agreement with the experimentally determined values than the former.

2.4. Guidelines for Further Optimization of β-Carbolines as Non-Covalent MAO-A Inhibitors

The low molecular weight of the investigated β-carboline compounds, as well as several substitution possibilities provided by their tricyclic core scaffold, suggests that these compounds could be further optimized to obtain even more active and pharmacokinetically more suitable preclinical candidates targeting the MAO-A enzyme. Furthermore, these β-carbolines can be classified as natural products which are known to be the most productive source of lead compounds to develop new drug molecules [60,61].
Previous docking studies of MAO-A and MAO-B enzymes with β-carboline derivatives have shown that the addition of lipophilic and bulky groups at C7 increases the inhibitory potency against MAO-A more than against MAO-B enzyme [62]. This is due to the structurally unfavorable overlap of these derivatives as well as some of the β-carbolines in this study substituted at the C7 position with Tyr326 residue of MAO-B, as compared to the structurally smaller Ile335 found in the MAO-A isoform [47]. The mentioned residues comprise one of the crucial differences between the two MAO isoforms and play an important role in determining the compound’s binding selectivity. The advantage of the β-carbolines substituted at C7 was further substantiated by our analysis of the MAO-A active site using molecular interaction fields (MIFs). The contours of the obtained MIF using a hydrophobic probe identified an empty lipophilic pocket and the possibility of hydrophobic interactions with residues Ile180, Phe208, Ile325, and Ile335. The other amino acid residues identified by our Apo Site Grid analysis (see Figure 7A) as capable of forming hydrophobic interactions were Tyr69, Ile207, Leu337, Met350, Phe352, and the aromatic cage tyrosines Tyr407 and Tyr444. For all these residues, our pharmacophore and dynophore analysis also revealed favorable hydrophobic interactions.
In addition, MIFs further revealed some MAO-A active site residues as favorable to form hydrogen bonds, namely Asn181, Gln215, Ser209, and Thr336—all of which, except for the latter, we also found to form interactions with the pyrrole or pyridine nitrogen using the dynophore analysis. However, since the active site is predominantly hydrophobic, hydrogen bonding is possible mainly with the water molecules present in the MAO-A binding cavity (for the average number of water molecules present in the binding site cavity during the MD simulations, see Supplementary Materials, Table S12). Based on all the collected data from the Apo Site Grid analysis as well as the previous structure and energy data from the MD simulations, we present some basic ideas for optimizing the β-carboline structure to fully utilize all the possible interactions within the MAO-A active site (see Figure 7B). The main possible points of optimization are carbon atoms C1 and C7, where the alkyl chain could potentially be extended to allow for more hydrophobic interactions with the active site. The C6 atom should be left unsubstituted, since the addition of functional groups on this atom does not contribute positively to the MAO-A inhibitory activity.

3. Discussion

With this computational study, we aimed to provide additional insight into MAO-A molecular recognition by examining the binding modes of a series of substituted β-carboline natural products. MAOs are enzymes with a variety of important functions in the central nervous system, so the development of drugs targeting these enzymes could be useful for several therapeutic applications. For instance, they show potential for development of novel drugs for treatment of neurodegeneration, while drugs targeting the MAO-A isoform are particularly useful in the treatment of mental disorders such as depression.
After docking a set of 11 β-carbolines into the MAO-A active site, the derived structure-based pharmacophores revealed key interactions with most of the amino acids identified as important in the X-ray-determined binding model of harmine 6, with an emphasis on the occurrence of hydrophobic interactions with the active site. There were minor differences in the binding modes of the neutral β-carbolines compared to their charged counterparts, with the charged β-carbolines exhibiting fewer interactions. While docking is invaluable for evaluating the geometric aspects of β-carboline binding, it is less adequate for predicting binding energies. Thus, we used the linear interaction energy (LIE) approach and generated MD simulations for all inhibitors as well as explored the molecular recognition feature with dynamic pharmacophores (i.e., dynophores). Similar to the obtained pharmacophores, the generated dynophores predicted mainly hydrophobic interactions with the active site residues along with some hydrogen bonding of the pyrrole and pyridine nitrogen with the water molecules present in the active site and with some amino acid residues. In addition, during the performed MD simulations, we also examined the positional changes in the MAO-A active site residues compared to their original positions in the X-ray structure. We found that no significant shifts in the most important residues (e.g., Tyr69, Ile180, Phe208, Leu337, Phe352, Tyr407, and Tyr444) occurred, suggesting that no substantial reorganization of the MAO-A active site is required for the non-covalent binding of β-carbolines.
Decomposition of the interaction energy into electrostatic and van der Waals contributions within the LIE formalism showed that van der Waals interactions appear to be the predominant driving force for β-carboline binding, consistent with the observed high hydrophobicity of the MAO-A binding site. The empirical LIE coefficients α, β, and γ were determined to be 0.54, 0.26, and −2.65, respectively, and the calculated free energies of binding were in good agreement with the experimentally determined values.
Dynamic pharmacophore analysis confirms that the relative orientation of the β-carbolines is important for their ability to form interactions with the targeted active site. Most importantly, stabilizing interactions with the tyrosines of the so-called aromatic cage (Tyr407 and Tyr444) are possible only when the β-carboline is in a parallel position such as the conformation obtained by X-ray crystallography. Other sites that appear to be important for the formation of interactions are the C1 atom, where the methyl group at this position can form hydrophobic interactions, and the C7 atom, where the methoxy group at this position forms additional hydrogen bonds with active site water molecules and various amino acid residues such as Ser209 and Cys323.
The results presented in this study thus provide a static and dynamic atomistic picture of molecular recognition occurring between substituted β-carbolines and MAO-A. Furthermore, the data obtained can serve as guidelines in further rational optimization and development of β-carbolines towards active compounds that would potentially be useful as MAO-A-specific inhibitors with therapeutic application in the treatment of depression.

4. Materials and Methods

4.1. Molecular Docking and Binding Site Analysis

Molecular docking studies were performed using the available co-crystal structure of the β-carboline harmine (compound 6) with the MAO-A enzyme isoform (PDB code 2Z5X) and the GOLD (v. 5.3.0) docking program [49]. The binding site included all atoms within a radius of 20 Å around the center, which was fixed to the N5 atom of the FAD moiety. Prior to docking, all water molecules were removed except for seven water molecules occupying the binding cavity, which were identified as important for the binding modes and were therefore retained, and all were used in the docking procedure [47]. The ligand obtained from the crystal structure (harmine, i.e., compound 6) was set as the reference ligand for the docking procedure and the conformations of the remaining ligands were generated as described in the next subsection. Standard protonation states of amino acid residues at neutral pH were determined. A set of 25 algorithm runs were performed for each ligand. To achieve optimal accuracy in the docking process, the following genetic algorithm settings were used: population size 100, selection pressure 1.1, number of operations 100,000, 5 islands, and niche size 2. The scoring functions used to evaluate the ligand positions obtained from the docking studies were GoldScore, ChemScore and ChemPLP. Docking results were visualized in LigandScout [50], which was also used to determine the 3D structure-based pharmacophores of the binding modes considered.
LigandScout was also used for the Apo Site Grid analysis with the default settings used. During this task, molecular probes, such as hydrogen bond acceptor, hydrogen bond donor, positive ionizable, negative ionizable, hydrophobic probes, and others, scanned the β-carboline binding site of the MAO-A enzyme, deriving contours of the corresponding molecular interaction fields (MIFs).

4.2. Preparation of the Initial Structures for Molecular Dynamics Simulations

Our MD studies were also all based on the experimentally solved crystal structure of the MAO-A enzyme with the bound reversible inhibitor harmine (compound 6) obtained from the Protein Data Bank (PDB code 2Z5X). The studied structures of β-carbolines 1–11 were prepared manually by modifying the experimentally determined compound 6 with PyMOL [63]. To obtain structures 9–11, the corresponding pyridine bond had to be hydrogenated, resulting in a partially unsaturated ring. Molecules 2, 4, 7, and 11 were modified by methylation of the pyridine ring nitrogen and molecules 5 and 8 by additional methylation of the pyrrole ring nitrogen resulting in six charged β-carbolines. The methoxy group of compound 6 had to be removed to give compounds 1–5. In addition, deletion of the methyl group yielded compound 1.
Before starting the MD simulations, the force field parameters had to be developed for all the compounds we studied. The simulations were performed using the OPLS-AA force field [64]. After all the compounds were extracted from their initial complexes, they were fully optimized at the Hartree-Fock (HF) level of theory using the 6–31G(d) basis set encoded in the Gaussian16 program package [65]. For all compounds, atomic charges were determined by fitting to the HF/6–31G(d) calculated electrostatic potential according to the RESP scheme as implemented in AmberTools18 [66]. Parameters were determined using the ffld_server utility assisted by the Maestro v. 11.7 graphical interface [67]. The MAO-A enzyme was parametrized following the procedure described in our previous work [68]. Parameters and atom types for all simulated ligands 1–11 are available in the Supporting Information (see Supplementary Materials, Tables S1–S11).
All crystal waters were removed from the original protein complex and a spherical cell of TIP3P waters with a radius of 30 Å centered on the crystal coordinates of the central N5 atom was constructed. Similarly, for the free-state simulations, a spherical cell of TIP3P waters with the same radius of 30 Å centered on the N5 atom of the ligands was prepared. For the positively charged ligands, a chloride ion was included in both the bound and free-state simulations to ensure electroneutrality. All topology and coordinate files needed to initiate the MD simulations were created using the qprep5 utility of the MD package Q v. 5.06 [69].

4.3. Molecular Dynamics Simulations

The calculation of binding free energies for our compounds was facilitated by generating MD trajectories of the studied compounds in the bound and free states using the qdyn5 utility of the MD package Q v. 5.06 [69].
In the free state, a position constraint of 100 kcal mol−1 Å−2 was imposed on the central nitrogen atom of all studied compounds to prevent diffusion of the compound to the edge of the simulation sphere. In addition, a flat-bottom harmonic potential of 75 kcal mol−1 Å−2 was placed on the chloride ion for the charged compounds in both states and set to zero for distances less than 26.5 Å from the center of the simulation sphere.
Protein atoms beyond the 30 Å spherical cell were restrained to their original crystal structure coordinates using harmonic constraints. Non-bonding interactions of these atoms were turned off. The local reaction field (LRF) method was used to evaluate long-range electrostatics for distances beyond the 10 Å cut-off. For distances below 10 Å, the non-bonding interactions were explicitly evaluated.
All systems were first relaxed in 14 separate steps, by slowly increasing the temperature and the step size. Initially, two separate simulations were performed with 30,000 and 20,000 steps at 5 K and a step size of 0.005 fs and 0.01 fs, respectively. For the bound state, all protein and ligand atoms were constrained to their original coordinates, allowing only the movement of water molecules. Then, the system was further relaxed in four additional stages of 20,000 steps and one of 50,000 steps, all at 5 K, with step sizes of 0.01 fs, 0.04 fs, 0.1 fs, 0.3 fs, and 1.0 fs, respectively. Then, the system was heated in five stages of 40,000 steps with a step size of 1.0 fs in 50 K increments starting at 50 K and ending at 250 K. Finally, two simulations of 10,0000 steps were performed at 298 K, with a step size of 1.0 and 2.0 fs, respectively. Data for the calculation of binding free energies were collected from 10 consecutive simulations with 5,000,000 steps and a step size of 2.0 fs at 298 K (NVT canonical ensemble), giving a total of 0.1 μs MD simulation time for each investigated compound in both the free and bound states.
Visualization and analysis of the generated MD trajectories were facilitated by the VMD program package [70]. RMSD analysis was performed for the entire trajectory, and the calculated RMSD values were calculated referring to the initial structure of the protein and ligand, respectively. All calculations were performed using the computational facilities of the Ažman High-Performance Computing Center of the National Institute of Chemistry.

4.4. Generation of Dynamical Pharmacophore Models

To obtain a more detailed description of our system in terms of dynamic pharmacophore (dynophore) models, 1000 exported MD equidistant frames were used in dynamic pharmacophore analysis with the DynophoreApp developed at the Molecular Design Lab of Freie Universität (FU) Berlin [58]. These dynophore calculations were performed on computers at the Molecular Design Lab, Berlin, Germany, and subsequently visualized in LigandScout [50].

4.5. Calculation of Binding Free Energies Using the Linear Interaction Energy Method

The linear interaction energy (LIE) method [63,71] was developed by Åqvist to facilitate the calculation of binding free energies of ligands that bind to their target in a non-covalent way (Gbind). We have successfully used it in our previous study, where we calculated the binding free energies of N-sulphonyl-glutamic acid inhibitors of MurD ligase [72]. LIE is based on a modified linear response approximation for electrostatic interactions and on an empirical term treating the non-polar interactions. The main advantage of LIE compared to other methods that can be used for the same purpose (e.g., free energy perturbation and thermodynamic integration) is that we only need to perform two MD simulations—one in the free state (i.e., solvated, in water) and one in the bound state (i.e., solvated protein). This means that the sampling of any intermediate states is not necessary which consequently drastically reduces the associated computational costs. A thermodynamic cycle can be drawn as outlined in Figure 8, where the top row represents the two physical states mentioned above and the bottom row represents two unphysical intermediate states.
An equation can be derived because the free energy is a state function and, therefore, the sum of the changes in the thermodynamic cycle must be zero:
G b i n d i n g + G d e c o u p l i n g P G d e c o u p l i n g W G b i n d i n g d e c o u p l e d = 0 .
Since, for a decoupled ligand, all non-bonding interactions are switched off, the G b i n d i n g d e c o u p l e d term is zero, and the equation can be rearranged as:
G b i n d i n g = G d e c o u p l i n g W G d e c o u p l i n g P .
The two terms on the right-hand side of the equation essentially represent the difference in free energy for the transfer of the ligand from its environment (W represents water or the free state, while P represents the protein or the bound state) to the gas phase. In the context of LIE, these can be written as:
G s o l v a t i o n Q = G d e c o u p l i n g Q = α V L Q v d W + β V L Q e l + γ ,
where the 〈⬚〉 brackets denote the MD mean values of the non-bonded van der Waals (vdW) and electrostatic (el) interactions between the ligand L and its environment Q (i.e., either the solvated protein binding site P or the aqueous solution W). Consequently, the final LIE equation has the following form:
G b i n d i n g = α ( V L P v d W V L W v d W ) + β ( V L P e l V L W e l ) + γ .
The empirical parameters of this equation are the coefficients α and β for the nonpolar and polar binding energy contributions, respectively, and an additional (optional) constant γ , which is used mainly when we are dealing with very hydrophobic active sites [73]. In our study, we determined the empirical coefficients α , β , and γ for our system using the linear approximation method based on the van der Waals and electrostatic energy averages in the bound and free states during the MD simulation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27196711/s1, Tables S1–S11: Parameters and atom types for compounds 1–11; Figures S1–S11: RMSD graphs for compounds 1–11; Figure S12: Overview of the interaction pattern obtained with dynophore analysis for compounds 2–5; Figure S13: Overview of the interaction pattern obtained with dynophore analysis for compounds 7–9; Figure S14: Comparison between X-ray binding pose of harmine with the determined binding modes; Table S12: The average number of water molecules present in the MAO-A active site during the simulation.

Author Contributions

Conceptualization, A.P. (Alja Prah), A.P. (Andrej Perdih), and J.M.; methodology, A.P. (Andrej Perdih); investigation, T.G.; writing—original draft preparation, A.P. (Alja Prah); writing—review and editing, A.P. (Alja Prah), A.P. (Andrej Perdih), M.S.D., and J.M.; visualization, A.P. (Alja Prah); supervision, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research Agency, grant number P1-0012 and the young researchers grant awarded to Alja Prah.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank Barbara Herlah from the National Institute of Chemistry, Slovenia, for her help with the Apo Site Grid analysis.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Shih, J.C.; Chen, K.; Ridd, M.J. Monoamine oxidase: From genes to behavior. Annu. Rev. Neurosci. 1999, 22, 197–217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Miller, J.R.; Edmondson, D.E. Structure-activity relationships in the oxidation of para-substituted benzylamine analogues by recombinant human liver monoamine oxidase A. Biochemistry 1999, 38, 13670–13683. [Google Scholar] [CrossRef] [PubMed]
  3. Silverman, R.B. Radical ideas about monoamine-oxidase. Acc. Chem. Res. 1995, 28, 335–342. [Google Scholar] [CrossRef]
  4. Abad, E.; Zenn, R.K.; Kästner, J. Reaction Mechanism of Monoamine Oxidase from QM/MM Calculations. J. Phys. Chem. B 2013, 117, 14238–14246. [Google Scholar] [CrossRef] [PubMed]
  5. Zapata-Torres, G.; Fierro, A.; Barriga-Gonzalez, G.; Salgado, J.C.; Celis-Barros, C. Revealing monoamine oxidase B catalytic mechanisms by means of the quantum chemical cluster approach. J. Chem. Inf. Model. 2015, 55, 1349–1360. [Google Scholar] [CrossRef] [Green Version]
  6. Ralph, E.C.; Hirschi, J.S.; Anderson, M.A.; Cleland, W.W.; Singleton, D.A.; Fitzpatrick, P.F. Insights into the mechanism of flavoprotein-catalyzed amine oxidation from nitrogen isotope effects on the reaction of N-methyltryptophan oxidase. Biochemistry 2007, 46, 7655–7664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Erdem, S.S.; Karahan, O.; Yildiz, I.; Yelekci, K. A computational study on the amine-oxidation mechanism of monoamine oxidase: Insight into the polar nucleophilic mechanism. Org. Biomol. Chem. 2006, 4, 646–658. [Google Scholar] [CrossRef]
  8. Cakir, K.; Erdem, S.S.; Atalay, V.E. ONIOM calculations on serotonin degradation by monoamine oxidase B: Insight into the oxidation mechanism and covalent reversible inhibition. Org. Biomol. Chem. 2016, 14, 9239–9252. [Google Scholar] [CrossRef]
  9. Akyuz, M.A.; Erdem, S.S. Computational modeling of the direct hydride transfer mechanism for the MAO catalyzed oxidation of phenethylamine and benzylamine: ONIOM (QM/QM) calculations. J. Neural Transm. 2013, 120, 937–945. [Google Scholar] [CrossRef]
  10. Atalay, V.E.; Erdem, S.S. A comparative computational investigation on the proton and hydride transfer mechanisms of monoamine oxidase using model molecules. Comput. Biol. Chem. 2013, 47, 181–191. [Google Scholar] [CrossRef] [PubMed]
  11. Kurtz, K.A.; Rishavy, M.A.; Cleland, W.W.; Fitzpatrick, P.F. Nitrogen isotope effects as probes of the mechanism of D-amino acid oxidase. J. Am. Chem. Soc. 2000, 122, 12896–12897. [Google Scholar] [CrossRef]
  12. Fitzpatrick, P.F. Oxidation of amines by flavoproteins. Arch. Biochem. Biophys. 2010, 493, 13–25. [Google Scholar] [CrossRef] [Green Version]
  13. Vianello, R.; Repič, M.; Mavri, J. How are biogenic amines metabolized by monoamine oxidases? Eur. J. Org. Chem. 2012, 36, 7057–7065. [Google Scholar] [CrossRef]
  14. Fowler, J.S.; Volkow, N.D.; Wang, G.J.; Pappas, N.; Logan, J.; Shea, C.; Alexoff, D.; MacGregor, R.R.; Schlyer, D.J.; Zezulkova, I.; et al. Brain monoamine oxidase A inhibition in cigarette smokers. Proc. Natl. Acad. Sci. USA 1996, 93, 14065–14069. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Fowler, J.S.; Volkow, N.D.; Wang, G.J.; Pappas, N.; Logan, J.; MacGregor, R.; Alexoff, D.; Shea, C.; Schlyer, D.; Wolf, A.P.; et al. Inhibition of monoamine oxidase B in the brains of smokers. Nature 1996, 379, 733–736. [Google Scholar] [CrossRef]
  16. Castagnoli, K.; Murugesan, T. Tobacco Leaf, Smoke and Smoking, MAO Inhibitors, Parkinson’s Disease and Neuroprotection; Are There Links? Neurotoxicology 2004, 25, 279–291. [Google Scholar] [CrossRef]
  17. Herraiz, T. Relative exposure to beta-carbolines norharman and harman from foods and tobacco smoke. Food Addit. Contam. 2004, 21, 1041–1050. [Google Scholar] [CrossRef] [PubMed]
  18. Herraiz, T.; Chaparro, C. Human monoamine oxidase is inhibited by tobacco smoke: Beta-carboline alkaloids act as potent and reversible inhibitors. Biochem. Biophys. Res. Commun. 2005, 326, 378–386. [Google Scholar] [CrossRef] [PubMed]
  19. Piechowska, P.; Zawirska-Wojtasiak, R.; Mildner-Szkudlarz, S. Bioactive beta-Carbolines in Food: A Review. Nutrients 2019, 11, 814. [Google Scholar] [CrossRef] [Green Version]
  20. Moura, D.J.; Richter, M.F.; Boeira, J.M.; Pegas Henriques, J.A.; Saffi, J. Antioxidant properties of beta-carboline alkaloids are related to their antimutagenic and antigenotoxic activities. Mutagenesis 2007, 22, 293–302. [Google Scholar] [CrossRef] [PubMed]
  21. Freedland, C.S.; Mansbach, R.S. Behavioral profile of constituents in ayahuasca, an Amazonian psychoactive plant mixture. Drug Alcohol Depend. 1999, 54, 183–194. [Google Scholar] [CrossRef]
  22. Boeira, J.M.; Viana, A.F.; Picada, J.N.; Henriques, J.A.P. Genotoxic and recombinogenic activities of the two β-carboline alkaloids harman and harmine in Saccharomyces cerevisiae. Mutat. Res. Fundam. Mol. Mech. Mutagen. 2002, 500, 39–48. [Google Scholar] [CrossRef]
  23. Husbands, S.M.; Glennon, R.A.; Gorgerat, S.; Gough, R.; Tyacke, R.; Crosby, J.; Nutt, D.J.; Lewis, J.W.; Hudson, A.L. β-carboline binding to imidazoline receptors. Drug Alcohol Depend. 2001, 64, 203–208. [Google Scholar] [CrossRef]
  24. Squires, P.E.; Hills, C.E.; Rogers, G.J.; Garland, P.; Farley, S.R.; Morgan, N.G. The putative imidazoline receptor agonist, harmane, promotes intracellular calcium mobilisation in pancreatic beta-cells. Eur. J. Pharmacol. 2004, 501, 31–39. [Google Scholar] [CrossRef] [PubMed]
  25. Grella, B.; Dukat, M.; Young, R.; Teitler, M.; Herrick-Davis, K.; Gauthier, C.B.; Glennon, R.A. Investigation of hallucinogenic and related β-carbolines. Drug Alcohol Depend. 1998, 50, 99–107. [Google Scholar] [CrossRef]
  26. Loew, G.H.; Nienow, J.; Lawson, J.A.; Toll, L.; Uyeno, E.T. Theoretical structure-activity studies of beta-carboline analogs. Requirements for benzodiazepine receptor affinity and antagonist activity. Mol. Pharmacol. 1985, 28, 17–31. [Google Scholar]
  27. Glennon, R.A.; Dukat, M.; Grella, B.; Hong, S.-S.; Costantino, L.; Teitler, M.; Smith, C.; Egan, C.; Davis, K.; Mattson, M.V. Binding of β-carbolines and related agents at serotonin (5-HT2 and 5-HT1A), dopamine (D2) and benzodiazepine receptors. Drug Alcohol Depend. 2000, 60, 121–132. [Google Scholar] [CrossRef]
  28. Pimpinella, G.; Palmery, M. Interaction of β-carbolines with central dopaminergic transmission in mice: Structure-activity relationships. Neurosci. Lett. 1995, 189, 121–124. [Google Scholar] [CrossRef]
  29. Matsubara, K.; Gonda, T.; Sawada, H.; Uezono, T.; Kobayashi, Y.; Kawamura, T.; Ohtaki, K.; Kimura, K.; Akaike, A. Endogenously occurring beta-carboline induces parkinsonism in nonprimate animals: A possible causative protoxin in idiopathic Parkinson’s disease. J. Neurochem. 1998, 70, 727–735. [Google Scholar] [CrossRef]
  30. Hayashi, K.; Nagao, M.; Sugimura, T. Interactions of norharman and harman with DNA. Nucleic Acids Res. 1977, 4, 3679–3685. [Google Scholar] [CrossRef] [Green Version]
  31. Funayama, Y.; Nishio, K.; Wakabayashi, K.; Nagao, M.; Shimoi, K.; Ohira, T.; Hasegawa, S.; Saijo, N. Effects of β- and γ-carboline derivatives on DNA topoisomerase activities. Mutat. Res. Fundam. Mol. Mech. Mutagen. 1996, 349, 183–191. [Google Scholar] [CrossRef]
  32. Balón, M.; Muñoz, M.A.; Carmona, C.; Guardado, P.; Galán, M. A fluorescence study of the molecular interactions of harmane with the nucleobases, their nucleosides and mononucleotides. Biophys. Chem. 1999, 80, 41–52. [Google Scholar] [CrossRef]
  33. Pari, K.; Sundari, C.S.; Chandani, S.; Balasubramanian, D. Beta-carbolines that accumulate in human tissues may serve a protective role against oxidative stress. J. Biol. Chem. 2000, 275, 2455–2462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Tsuchiya, H.; Sato, M.; Watanabe, I. Antiplatelet activity of soy sauce as functional seasoning. J. Agric. Food. Chem. 1999, 47, 4167–4174. [Google Scholar] [CrossRef] [PubMed]
  35. Ancolio, C.; Azas, N.; Mahiou, V.; Ollivier, E.; Di Giorgio, C.; Keita, A.; Timon-David, P.; Balansard, G. Antimalarial activity of extracts and alkaloids isolated from six plants used in traditional medicine in Mali and Sao Tome. Phytother. Res. 2002, 16, 646–649. [Google Scholar] [CrossRef]
  36. Ferraz, C.A.A.; de Oliveira Junior, R.G.; Picot, L.; da Silva Almeida, J.R.G.; Nunes, X.P. Pre-clinical investigations of beta-carboline alkaloids as antidepressant agents: A systematic review. Fitoterapia 2019, 137, 104196. [Google Scholar] [CrossRef]
  37. Aricioglu, F.; Altunbas, H. Harmane induces anxiolysis and antidepressant-like effects in rats. Ann. N. Y. Acad. Sci. 2003, 1009, 196–201. [Google Scholar] [CrossRef] [PubMed]
  38. Smith, K.L.; Ford, G.K.; Jessop, D.S.; Finn, D.P. Behavioural, neurochemical and neuroendocrine effects of the endogenous beta-carboline harmane in fear-conditioned rats. J. Psychopharm. 2013, 27, 162–170. [Google Scholar] [CrossRef] [PubMed]
  39. Farzin, D.; Mansouri, N. Antidepressant-like effect of harmane and other beta-carbolines in the mouse forced swim test. Eur. Neuropsychopharmacol. 2006, 16, 324–328. [Google Scholar] [CrossRef] [PubMed]
  40. Adell, A.; Biggs, T.A.; Myers, R.D. Action of harman (1-methyl-β-carboline) on the brain: Body temperature and in vivo efflux of 5-HT from hippocampus of the rat. Neuropharmacology 1996, 35, 1101–1107. [Google Scholar] [CrossRef]
  41. Fortunato, J.J.; Reus, G.Z.; Kirsch, T.R.; Stringari, R.B.; Stertz, L.; Kapczinski, F.; Pinto, J.P.; Hallak, J.E.; Zuardi, A.W.; Crippa, J.A.; et al. Acute harmine administration induces antidepressive-like effects and increases BDNF levels in the rat hippocampus. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 2009, 33, 1425–1430. [Google Scholar] [CrossRef] [PubMed]
  42. Fortunato, J.J.; Reus, G.Z.; Kirsch, T.R.; Stringari, R.B.; Fries, G.R.; Kapczinski, F.; Hallak, J.E.; Zuardi, A.W.; Crippa, J.A.; Quevedo, J. Effects of beta-carboline harmine on behavioral and physiological parameters observed in the chronic mild stress model: Further evidence of antidepressant properties. Brain Res. Bull. 2010, 81, 491–496. [Google Scholar] [CrossRef]
  43. WHO. Depression. Available online: https://www.who.int/news-room/fact-sheets/detail/depression (accessed on 13 September 2022).
  44. Cao, R.; Peng, W.; Wang, Z.; Xu, A. Beta-Carboline alkaloids: Biochemical and pharmacological functions. Curr. Med. Chem. 2007, 14, 479–500. [Google Scholar] [CrossRef] [PubMed]
  45. Varela, A.P.; Burrows, H.D.; Douglas, P.; da Graça Miguel, M. Triplet state studies of β-carbolines. J. Photochem. Photobiol. A Chem. 2001, 146, 29–36. [Google Scholar] [CrossRef]
  46. May, T.; Rommelspacher, H.; Pawlik, M. [3H]harman binding experiments. I: A reversible and selective radioligand for monoamine oxidase subtype A in the CNS of the rat. J. Neurochem. 1991, 56, 490–499. [Google Scholar] [CrossRef] [PubMed]
  47. Son, S.Y.; Ma, A.; Kondou, Y.; Yoshimura, M.; Yamashita, E.; Tsukihara, T. Structure of human monoamine oxidase A at 2.2-angstrom resolution: The control of opening the entry for substrates/inhibitors. Proc. Natl. Acad. Sci. USA 2008, 105, 5739–5744. [Google Scholar] [CrossRef] [Green Version]
  48. Kim, H.; Sablin, S.O.; Ramsay, R.R. Inhibition of monoamine oxidase A by beta-carboline derivatives. Arch. Biochem. Biophys. 1997, 337, 137–142. [Google Scholar] [CrossRef] [PubMed]
  49. Jones, G.; Willett, P.; Glen, R.C.; Leach, A.R.; Taylor, R. Development and validation of a genetic algorithm for flexible docking. J. Mol. Biol. 1997, 267, 727–748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Wolber, G.; Langer, T. LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters. J. Chem. Inf. Model. 2005, 45, 160–169. [Google Scholar] [CrossRef]
  51. Wang, R.; Lu, Y.; Wang, S. Comparative evaluation of 11 scoring functions for molecular docking. J. Med. Chem. 2003, 46, 2287–2303. [Google Scholar] [CrossRef] [PubMed]
  52. King, E.; Aitchison, E.; Li, H.; Luo, R. Recent Developments in Free Energy Calculations for Drug Discovery. Front. Mol. Biosci. 2021, 8, 712085. [Google Scholar] [CrossRef] [PubMed]
  53. van Lipzig, M.M.; ter Laak, A.M.; Jongejan, A.; Vermeulen, N.P.; Wamelink, M.; Geerke, D.; Meerman, J.H. Prediction of ligand binding affinity and orientation of xenoestrogens to the estrogen receptor by molecular dynamics simulations and the linear interaction energy method. J. Med. Chem. 2004, 47, 1018–1030. [Google Scholar] [CrossRef] [PubMed]
  54. Stjernschantz, E.; Oostenbrink, C. Improved ligand-protein binding affinity predictions using multiple binding modes. Biophys. J. 2010, 98, 2682–2691. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Wang, W.; Wang, J.; Kollman, P.A. What determines the van der Waals coefficient in the LIE (linear interaction energy) method to estimate binding free energies using molecular dynamics simulations? Proteins Struct. Funct. Genet. 1999, 34, 395–402. [Google Scholar] [CrossRef]
  56. Bock, A.; Bermudez, M.; Krebs, F.; Matera, C.; Chirinda, B.; Sydow, D.; Dallanoce, C.; Holzgrabe, U.; De Amici, M.; Lohse, M.J.; et al. Ligand Binding Ensembles Determine Graded Agonist Efficacies at a G Protein-coupled Receptor. J. Biol. Chem. 2016, 291, 16375–16389. [Google Scholar] [CrossRef] [Green Version]
  57. Bermudez, M.; Rakers, C.; Wolber, G. Structural Characteristics of the Allosteric Binding Site Represent a Key to Subtype Selective Modulators of Muscarinic Acetylcholine Receptors. Mol. Inform. 2015, 34, 526–530. [Google Scholar] [CrossRef]
  58. Sydow, D. Dynophores: Novel Dynamic Pharmacophores; Lebenswissenschaftliche Fakultät: Berlin, Germany, 2015. [Google Scholar]
  59. Akyüz, M.A.; Erdem, S.S.; Edmondson, D.E. The aromatic cage in the active site of monoamine oxidase B: Effect on the structural and electronic properties of bound benzylamine and p-nitrobenzylamine. J. Neural Transm. 2007, 114, 693–698. [Google Scholar] [CrossRef]
  60. Harvey, A.L. Natural products in drug discovery. Drug Discov. Today 2008, 13, 894–901. [Google Scholar] [CrossRef]
  61. Atanasov, A.G.; Zotchev, S.B.; Dirsch, V.M.; Supuran, C.T.; The International Natural Product Sciences Taskforce. Natural products in drug discovery: Advances and opportunities. Nat. Rev. Drug Discov. 2021, 20, 200–216. [Google Scholar] [CrossRef]
  62. Reniers, J.; Robert, S.; Frederick, R.; Masereel, B.; Vincent, S.; Wouters, J. Synthesis and evaluation of beta-carboline derivatives as potential monoamine oxidase inhibitors. Bioorg. Med. Chem. 2011, 19, 134–144. [Google Scholar] [CrossRef]
  63. Aqvist, J.; Medina, C.; Samuelsson, J.E. A new method for predicting binding affinity in computer-aided drug design. Protein Eng. 1994, 7, 385–391. [Google Scholar] [CrossRef] [PubMed]
  64. Jorgensen, W.L.; Maxwell, D.S.; TiradoRives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  65. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  66. Case, D.A.; Ben-Shalom, I.Y.; Brozell, S.R.; Cerutti, D.S.; Cheatham, T.; Cruzeiro, V.W.D.; Darden, T.A.; Duke, R.E.; Ghoreishi, D.; Gilson, M.K.; et al. AMBER 2018; University of California: San Francisco, CA, USA, 2018. [Google Scholar]
  67. Perrin, D.D.; Dempsey, B.; Serjeant, E.P. pKa Prediction for Organic Acids and Bases; Chapman and Hall: London, UK; New York, NY, USA, 1981; 146p. [Google Scholar]
  68. Poberžnik, M.; Purg, M.; Repič, M.; Mavri, J.; Vianello, R. Empirical valence bond simulations of the hydride-transfer step in the monoamine oxidase A catalyzed metabolism of noradrenaline. J. Phys. Chem. B 2016, 120, 11419–11427. [Google Scholar] [CrossRef]
  69. Marelius, J.; Kolmodin, K.; Feierberg, I.; Aqvist, J. Q: A molecular dynamics program for free energy calculations and empirical valence bond simulations in biomolecular systems. J. Mol. Graph. Model. 1998, 16, 213–225, 261. [Google Scholar] [CrossRef]
  70. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. Modell. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  71. Aqvist, J.; Luzhkov, V.B.; Brandsdal, B.O. Ligand binding affinities from MD simulations. Acc. Chem. Res. 2002, 35, 358–365. [Google Scholar] [CrossRef]
  72. Perdih, A.; Bren, U.; Solmajer, T. Binding free energy calculations of N-sulphonyl-glutamic acid inhibitors of MurD ligase. J. Mol. Model. 2009, 15, 983–996. [Google Scholar] [CrossRef] [PubMed]
  73. Aqvist, J.; Marelius, J. The linear interaction energy method for predicting ligand binding free energies. Comb. Chem. High Throughput Screen. 2001, 4, 613–626. [Google Scholar] [CrossRef]
Figure 1. Left: monoamine oxidase A enzyme. Top right: Active site of monoamine oxidase A with the flavin adenine nucleotide (FAD) cofactor and the bound β-carboline harmine (HRM). Bottom right: basic tricyclic structural element and the numbering of the β-carbolines.
Figure 1. Left: monoamine oxidase A enzyme. Top right: Active site of monoamine oxidase A with the flavin adenine nucleotide (FAD) cofactor and the bound β-carboline harmine (HRM). Bottom right: basic tricyclic structural element and the numbering of the β-carbolines.
Molecules 27 06711 g001
Figure 2. (A) 3D (top) and 2D (bottom) structure-based pharmacophore model for harmine (compound 6) docked in the MAO-A active site. (B)) 3D (top) and 2D (bottom) structure-based pharmacophore model for 6-methoxy harmalan (compound 10) docked in the MAO-A active site. Yellow spheres represent favorable areas of the ligand for hydrophobic interactions while green and red arrows denote the identified H-bond donors and acceptor interactions, respectively. (C) GoldScore, ChemScore and ChemPLP scoring functions compared to experimentally observed binding free energy of compounds 1–11.
Figure 2. (A) 3D (top) and 2D (bottom) structure-based pharmacophore model for harmine (compound 6) docked in the MAO-A active site. (B)) 3D (top) and 2D (bottom) structure-based pharmacophore model for 6-methoxy harmalan (compound 10) docked in the MAO-A active site. Yellow spheres represent favorable areas of the ligand for hydrophobic interactions while green and red arrows denote the identified H-bond donors and acceptor interactions, respectively. (C) GoldScore, ChemScore and ChemPLP scoring functions compared to experimentally observed binding free energy of compounds 1–11.
Molecules 27 06711 g002
Figure 3. (A) Outline of the defined bound and free state for compound 2. (B) Interaction energy decomposition for compound 2: van der Waals (top, vdW) and electrostatic (bottom, el. stat.) energy component between compound 2 and its surroundings—water in the free state (blue) and solvated protein in the bound state (orange) during the whole molecular dynamics simulation. (C) Experimental- and LIE-calculated binding free energies together with error bars for compounds 1–11.
Figure 3. (A) Outline of the defined bound and free state for compound 2. (B) Interaction energy decomposition for compound 2: van der Waals (top, vdW) and electrostatic (bottom, el. stat.) energy component between compound 2 and its surroundings—water in the free state (blue) and solvated protein in the bound state (orange) during the whole molecular dynamics simulation. (C) Experimental- and LIE-calculated binding free energies together with error bars for compounds 1–11.
Molecules 27 06711 g003
Figure 4. Overview of the interaction pattern obtained through dynophore analysis for β-carboline compound 1. (A) Calculated percentage of the occurrence of a dynamic pharmacophore element based on all frames divided into element-interacting amino acid pairs (for compounds 2–5, see Supplementary Materials, Figure S12). (B) The presence of hydrophobic interactions of the pyridine ring with various active site amino acid residues during the molecular dynamics simulation. (C) RMSD of compound 1 during the molecular dynamics simulation (for compounds 2–5, see Supplementary Materials, Figures S2–S5). (D) Distance between the pyrrole nitrogen and residue Asn181 during the molecular dynamics simulation for compounds 1–5. The distance is shown as a block average over 1 ns (1000 total blocks of 0.1 ns cover the entire 100 ns MD simulation).
Figure 4. Overview of the interaction pattern obtained through dynophore analysis for β-carboline compound 1. (A) Calculated percentage of the occurrence of a dynamic pharmacophore element based on all frames divided into element-interacting amino acid pairs (for compounds 2–5, see Supplementary Materials, Figure S12). (B) The presence of hydrophobic interactions of the pyridine ring with various active site amino acid residues during the molecular dynamics simulation. (C) RMSD of compound 1 during the molecular dynamics simulation (for compounds 2–5, see Supplementary Materials, Figures S2–S5). (D) Distance between the pyrrole nitrogen and residue Asn181 during the molecular dynamics simulation for compounds 1–5. The distance is shown as a block average over 1 ns (1000 total blocks of 0.1 ns cover the entire 100 ns MD simulation).
Molecules 27 06711 g004
Figure 5. Overview of the interaction pattern obtained through dynophore analysis for β-carboline harmine 6. (A) Calculated percentage of the occurrence of a dynamic pharmacophore element based on all frames divided into element-interacting amino acids pairs (for compounds 7–9, see Supplementary Materials, Figure S13). (B) The presence of hydrophobic interactions of the pyridine ring with various active site amino acid residues during the molecular dynamics simulation. (C) RMSD of compound 6 during the molecular dynamics simulation (for all compounds, see Supplementary Materials, Figures S1–S11). (D) Distance between the pyrrole nitrogen and Asn181 during the molecular dynamics simulation.
Figure 5. Overview of the interaction pattern obtained through dynophore analysis for β-carboline harmine 6. (A) Calculated percentage of the occurrence of a dynamic pharmacophore element based on all frames divided into element-interacting amino acids pairs (for compounds 7–9, see Supplementary Materials, Figure S13). (B) The presence of hydrophobic interactions of the pyridine ring with various active site amino acid residues during the molecular dynamics simulation. (C) RMSD of compound 6 during the molecular dynamics simulation (for all compounds, see Supplementary Materials, Figures S1–S11). (D) Distance between the pyrrole nitrogen and Asn181 during the molecular dynamics simulation.
Molecules 27 06711 g005
Figure 6. Overview of the interaction pattern obtained by dynophore analysis for compounds 10 (top) and 11 (bottom). Calculated percentage of the occurrence of a dynamic pharmacophore element based on all frames divided into element-interacting amino acids pairs for compound 10 (A) and compound 11 (C). Calculated RMSD value for the duration of the molecular dynamics simulation for compound 10 (B) and compound 11 (D).
Figure 6. Overview of the interaction pattern obtained by dynophore analysis for compounds 10 (top) and 11 (bottom). Calculated percentage of the occurrence of a dynamic pharmacophore element based on all frames divided into element-interacting amino acids pairs for compound 10 (A) and compound 11 (C). Calculated RMSD value for the duration of the molecular dynamics simulation for compound 10 (B) and compound 11 (D).
Molecules 27 06711 g006
Figure 7. (A) Calculated molecular interaction fields in the MAO-A active site using different molecular probes. Regions of the active site with amino acids capable of acting as hydrogen bond donors (green), acceptors (red) or forming hydrophobic interactions (yellow). (B) General β-carboline structure with proposed possibilities for rational structure-based optimization of MAO-A inhibitory activity.
Figure 7. (A) Calculated molecular interaction fields in the MAO-A active site using different molecular probes. Regions of the active site with amino acids capable of acting as hydrogen bond donors (green), acceptors (red) or forming hydrophobic interactions (yellow). (B) General β-carboline structure with proposed possibilities for rational structure-based optimization of MAO-A inhibitory activity.
Molecules 27 06711 g007
Figure 8. Thermodynamic cycle used in LIE binding free energy calculations. P stands for protein, L for ligand and W for water. The dashed line around the ligand represents decoupling (turning the non-bonding electrostatic and van der Waals interactions off).
Figure 8. Thermodynamic cycle used in LIE binding free energy calculations. P stands for protein, L for ligand and W for water. The dashed line around the ligand represents decoupling (turning the non-bonding electrostatic and van der Waals interactions off).
Molecules 27 06711 g008
Table 1. Average van der Waals and electrostatic interaction energy between β-carbolines 1–11 (L) and their surroundings (solvated protein P or aqueous solution W) obtained from MD simulations. Standard deviations are also given and were calculated from ten subsequent molecular dynamics simulation steps.
Table 1. Average van der Waals and electrostatic interaction energy between β-carbolines 1–11 (L) and their surroundings (solvated protein P or aqueous solution W) obtained from MD simulations. Standard deviations are also given and were calculated from ten subsequent molecular dynamics simulation steps.
Compound V L P v d W [ kcal   mol −1 ] V L W v d W [ kcal   mol −1 ] V L P v d W V L W v d W [ kcal   mol −1 ] V L P e l [ kcal   mol −1 ] V L W e l [ kcal   mol −1 ] V L P e l V L W e l [ kcal   mol −1 ]
1−28.4 ± 0.5−15.4 ± 0.1−13.1 ± 0.5−15.1 ± 1.6−22.2 ± 0.1+7.1 ± 1.6
2 *−36.7 ± 0.5−17.9 ± 0.1−18.9 ± 0.5−84.1 ± 2.2−98.0 ± 0.2+13.9 ± 2.2
3−32.5 ± 0.3−16.8 ± 0.1−15.7 ± 0.3−15.7 ± 0.9−23.6 ± 0.1+7.9 ± 0.9
4 *−38.2 ± 0.8−19.4 ± 0.1−17.7 ± 0.8−82.9 ± 1.6−96.9 ± 0.1+14.0 ± 1.6
5 *−39.8 ± 1.0−22.7 ± 0.1−17.1 ± 1.0−77.6 ± 0.7−86.5 ± 0.2+8.9 ± 0.7
6−36.9 ± 0.4−19.6 ± 0.1−17.2 ± 0.4−17.0 ± 1.2−24.5 ± 0.1+7.2 ± 1.2
7 *−42.3 ± 1.0−22.0 ± 0.1−20.3 ± 1.0−78.9 ± 4.1−96.7 ± 0.2+17.8 ± 4.1
8 *−43.9 ± 0.7−25.3 ± 0.1−18.6 ± 0.7−74.5 ± 1.6−87.5 ± 0.2+12.9 ± 1.6
9−41.9 ± 1.0−20.6 ± 0.1−21.2 ± 1.0−9.7 ± 0.6−23.4 ± 0.1+13.7 ± 0.6
10−39.5 ± 0.6−20.7 ± 0.1−18.8 ± 0.6−12.9 ± 0.7−23.3 ± 0.1+10.4 ± 0.7
11 *−41.0 ± 0.2−22.6 ± 0.1−18.4 ± 0.2−81.1 ± 1.1−96.8 ± 0.2+15.7 ± 1.1
* Inhibitors that are charged.
Table 2. Calculated binding free energies of β-carbolines 1–11 compared with the available experimental results (obtained from ref. [48]) and average RMSD values for each ligand (for RMSD graphs of all compounds along the duration of molecular dynamics simulations see Supplementary Materials, Figures S1–S11).
Table 2. Calculated binding free energies of β-carbolines 1–11 compared with the available experimental results (obtained from ref. [48]) and average RMSD values for each ligand (for RMSD graphs of all compounds along the duration of molecular dynamics simulations see Supplementary Materials, Figures S1–S11).
CompoundStructure K i [ µ M ] G b i n d i n g e x p [ kcal   mol −1 ] G b i n d i n g c a l c [ kcal   mol −1 ] RMSD
[Å]
1Molecules 27 06711 i0013.34 ± 0.10−7.6 ± 0.02−7.9 ± 0.53.7
2Molecules 27 06711 i0021.43 ± 0.09−8.1 ± 0.04−9.3 ± 0.61.1
3Molecules 27 06711 i0030.26 ± 0.024−9.1 ± 0.06−9.1 ± 0.31.9
4Molecules 27 06711 i0040.68 ± 0.026−8.6 ± 0.02−8.6 ± 0.61.5
5Molecules 27 06711 i0050.16 ± 0.036−9.4 ± 0.14−9.6 ± 0.61.6
6Molecules 27 06711 i0060.005 ± 0.0002−11.5 ± 0.02−10.0 ± 0.42.3
7Molecules 27 06711 i0070.069 ± 0.008−9.9 ± 0.07−9.0 ± 1.21.7
8Molecules 27 06711 i0080.015 ± 0.0008−10.9 ± 0.03−9.4 ± 0.61.4
9Molecules 27 06711 i0090.048 ± 0.007−10.2 ± 0.09−10.6 ± 0.61.9
10Molecules 27 06711 i0100.39 ± 0.052−8.9 ± 0.08−10.1 ± 0.43.7
11Molecules 27 06711 i0111.23 ± 0.27−8.2 ± 0.01−8.5 ± 0.33.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Prah, A.; Gavranić, T.; Perdih, A.; Sollner Dolenc, M.; Mavri, J. Computational Insights into β-Carboline Inhibition of Monoamine Oxidase A. Molecules 2022, 27, 6711. https://doi.org/10.3390/molecules27196711

AMA Style

Prah A, Gavranić T, Perdih A, Sollner Dolenc M, Mavri J. Computational Insights into β-Carboline Inhibition of Monoamine Oxidase A. Molecules. 2022; 27(19):6711. https://doi.org/10.3390/molecules27196711

Chicago/Turabian Style

Prah, Alja, Tanja Gavranić, Andrej Perdih, Marija Sollner Dolenc, and Janez Mavri. 2022. "Computational Insights into β-Carboline Inhibition of Monoamine Oxidase A" Molecules 27, no. 19: 6711. https://doi.org/10.3390/molecules27196711

APA Style

Prah, A., Gavranić, T., Perdih, A., Sollner Dolenc, M., & Mavri, J. (2022). Computational Insights into β-Carboline Inhibition of Monoamine Oxidase A. Molecules, 27(19), 6711. https://doi.org/10.3390/molecules27196711

Article Metrics

Back to TopTop