To obtain insights into the dynamics of nutrient exchange in arbuscular mycorrhizal (AM) symbiosis, we modelled mathematically the two-membrane system at the plant–fungus interface and simulated its ...dynamics.
In computational cell biology experiments, the full range of nutrient transport pathways was tested for their ability to exchange phosphorus (P)/carbon (C)/nitrogen (N) sources.
As a result, we obtained a thermodynamically justified, independent and comprehensive model of the dynamics of the nutrient exchange at the plant–fungus contact zone. The predicted optimal transporter network coincides with the transporter set independently confirmed in wet-laboratory experiments previously, indicating that all essential transporter types have been discovered.
The thermodynamic analyses suggest that phosphate is released from the fungus via proton-coupled phosphate transporters rather than anion channels. Optimal transport pathways, such as cation channels or proton-coupled symporters, shuttle nutrients together with a positive charge across the membranes. Only in exceptional cases does electroneutral transport via diffusion facilitators appear to be plausible. The thermodynamic models presented here can be generalized and adapted to other forms of mycorrhiza and open the door for future studies combining wet-laboratory experiments with computational simulations to obtain a deeper understanding of the investigated phenomena.
In the last years, the interactions of flavonoids with protein kinases (PKs) have been described by using crystallographic experiments. Interestingly, different orientations have been found for one ...flavonoid inside different PKs and different chemical substitutions lead to different orientations of the flavonoid scaffold inside one PK. Accordingly, orientation predictions of novel analogues could help to the design of flavonoids with high PK inhibitory activities. With this in mind, we studied the binding modes of 37 flavonoids (flavones and chalcones) inside the cyclin-dependent PK CDK1 using docking experiments. We found that the compounds under study adopted two different orientations into the active site of CDK1 (orientations I and II in the manuscript). In addition, quantitative structure-activity relationship (QSAR) models using CoMFA and CoMSIA methodologies were constructed to explain the trend of the CDK1 inhibitory activities for the studied flavonoids. Template-based and docking-based alignments were used. Models developed starting from docking-based alignment were applied for describing the whole dataset and compounds with orientation I. Adequate R2 and Q2 values were obtained by each method; interestingly, only hydrophobic and hydrogen bond donor fields describe the differential potency of the flavonoids as CDK1 inhibitors for both defined alignments and subsets. Our current application of docking and QSAR together reveals important elements to be drawn for the design of novel flavonoids with increased PK inhibitory activities.
Two Pore Channels (TPCs) are cation-selective voltage- and ligand-gated ion channels in membranes of intracellular organelles of eukaryotic cells. In plants, the TPC1 subtype forms the slowly ...activating vacuolar (SV) channel, the most dominant ion channel in the vacuolar membrane. Controversial reports about the permeability properties of plant SV channels fueled speculations about the physiological roles of this channel type. TPC1 is thought to have high Ca2+ permeability, a conclusion derived from relative permeability analyses using the Goldman–Hodgkin–Katz (GHK) equation. Here, we investigated in computational analyses the properties of the permeation pathway of TPC1 from Arabidopsis thaliana. Using the crystal structure of AtTPC1, protein modeling, molecular dynamics (MD) simulations, and free energy calculations, we identified a free energy minimum for Ca2+, but not for K+, at the luminal side next to the selectivity filter. Residues D269 and E637 coordinate in particular Ca2+ as demonstrated in in silico mutagenesis experiments. Such a Ca2+-specific coordination site in the pore explains contradicting data for the relative Ca2+/K+ permeability and strongly suggests that the Ca2+ permeability of SV channels is largely overestimated from relative permeability analyses. This conclusion was further supported by in silico electrophysiological studies showing a remarkable permeation of K+ but not Ca2+ through the open channel.
The high rates of morbidity and mortality due to fungal infections are associated with a limited antifungal arsenal and the high toxicity of drugs. Therefore, the identification of novel drug targets ...is challenging due to the several resemblances between fungal and human cells. Here, we report the in vitro antifungal evaluation of two acylphenols series, namely 2-acyl-1,4-benzo- and 2-acyl-1,4-naphthohydroquinones. The antifungal properties were assessed on diverse
and filamentous fungi strains through the halo of inhibition (HOI) and minimal inhibitory concentration (MIC). The antifungal activities of 2-acyl-1,4-benzohydroquinone derivatives were higher than those of the 2-acyl-1,4-naphthohydroquinone analogues. The evaluation indicates that 2-octanoylbenzohydroquinone
is the most active member of the 2-acylbenzohydroquinone series, with MIC values ranging from 2 to 16 μg/mL. In some fungal strains (i.e.,
and
), such MIC values of compound
(2 and 4 μg/mL) were comparable to that obtained by amphotericin B (1 μg/mL). The compound
was evaluated for its antioxidant activity by means of FRAP, ABTS and DPPH assays, showing moderate activity as compared to standard antioxidants. Molecular docking studies of compound
and ADMET predictions make this compound a potential candidate for topical pharmacological use. The results obtained using the most active acylbenzohydroquinones are promising because some evaluated
strains are known to have decreased sensitivity to standard antifungal treatments.
Aroma and flavor are important factors of fruit quality and consumer preference. The specific pattern of aroma is generated during ripening by the accumulation of volatiles compounds, which are ...mainly esters. Alcohol acyltransferase (AAT) (EC 2.3.1.84) catalyzes the esterification reaction of aliphatic and aromatic alcohols and acyl-CoA into esters in fruits and flowers. In Fragaria x ananassa, there are different volatiles compounds that are obtained from different alcohol precursors, where octanol and hexanol are the most abundant during fruit ripening. At present, there is not structural evidence about the mechanism used by the AAT to synthesize esters. Experimental data attribute the kinetic role of this enzyme to 2 amino acidic residues in a highly conserved motif (HXXXD) that is located in the middle of the protein. With the aim to understand the molecular and energetic aspects of volatiles compound production from F. x ananassa, we first studied the binding modes of a series of alcohols, and also different acyl-CoA substrates, in a molecular model of alcohol acyltransferase from Fragaria x ananassa (SAAT) using molecular docking. Afterwards, the dynamical behavior of both substrates, docked within the SAAT binding site, was studied using routine molecular dynamics (MD) simulations. In addition, in order to correlate the experimental and theoretical data obtained in our laboratories, binding free energy calculations were performed; which previous results suggested that octanol, followed by hexanol, presented the best affinity for SAAT. Finally, and concerning the SAAT molecular reaction mechanism, it is suggested from molecular dynamics simulations that the reaction mechanism may proceed through the formation of a ternary complex, in where the Histidine residue at the HXXXD motif deprotonates the alcohol substrates. Then, a nucleophilic attack occurs from alcohol charged oxygen atom to the carbon atom at carbonyl group of the acyl CoA. This mechanism is in agreement with previous results, obtained in our group, in alcohol acyltransferase from Vasconcellea pubescens (VpAAT1).
Human arginase I (hARGI) is an important enzyme involved in the urea cycle; its overexpression has been associated to cardiovascular and cerebrovascular diseases. In the last years, several ...congeneric sets of hARGI inhibitors have been reported with possible beneficial roles for the cardiovascular system. At the same time, crystallographic data have been reported including hARGI⁻inhibitor complexes, which can be considered for the design of novel inhibitors. In this work, the structure⁻activity relationship (SAR) of Cα substituted 2(S)-amino-6-boronohexanoic acid (ABH) derivatives as hARGI inhibitors was studied by using a three-dimensional quantitative structure⁻activity relationships (3D-QSAR) method. The predictivity of the obtained 3D-QSAR model was demonstrated by using internal and external validation experiments. The best model revealed that the differential hARGI inhibitory activities of the
derivatives can be described by using steric and electrostatic fields; the local effects of these fields in the activity are presented. In addition, binding modes of the above-mentioned compounds inside the hARGI binding site were obtained by using molecular docking. It was found that
derivatives adopted the same orientation reported for
within the hARGI active site, with the substituents at Cα exposed to the solvent with interactions with residues at the entrance of the binding site. The hARGI residues involved in chemical interactions with inhibitors were identified by using an interaction fingerprints (IFPs) analysis.
To fire action-potential-like electrical signals, the vacuole membrane requires the two-pore channel TPC1, formerly called SV channel. The TPC1/SV channel functions as a depolarization-stimulated, ...non-selective cation channel that is inhibited by luminal Ca2+. In our search for species-dependent functional TPC1 channel variants with different luminal Ca2+ sensitivity, we found in total three acidic residues present in Ca2+ sensor sites 2 and 3 of the Ca2+-sensitive AtTPC1 channel from Arabidopsis thaliana that were neutral in its Vicia faba ortholog and also in those of many other Fabaceae. When expressed in the Arabidopsis AtTPC1-loss-of-function background, wild-type VfTPC1 was hypersensitive to vacuole depolarization and only weakly sensitive to blocking luminal Ca2+. When AtTPC1 was mutated for these VfTPC1-homologous polymorphic residues, two neutral substitutions in Ca2+ sensor site 3 alone were already sufficient for the Arabidopsis At-VfTPC1 channel mutant to gain VfTPC1-like voltage and luminal Ca2+ sensitivity that together rendered vacuoles hyperexcitable. Thus, natural TPC1 channel variants exist in plant families which may fine-tune vacuole excitability and adapt it to environmental settings of the particular ecological niche.
It is evident that 99mTc causes radical-mediated DNA damage due to Auger electrons, which were emitted simultaneously with the known γ-emission of 99mTc. We have synthesized a series of new ...99mTc-labeled pyrene derivatives with varied distances between the pyrene moiety and the radionuclide. The pyrene motif is a common DNA intercalator and allowed us to test the influence of the radionuclide distance on damages of the DNA helix. In general, pUC 19 plasmid DNA enables the investigation of the unprotected interactions between the radiotracers and DNA that results in single-strand breaks (SSB) or double-strand breaks (DSB). The resulting DNA fragments were separated by gel electrophoresis and quantified by fluorescent staining. Direct DNA damage and radical-induced indirect DNA damage by radiolysis products of water were evaluated in the presence or absence of the radical scavenger DMSO. We demonstrated that Auger electrons directly induced both SSB and DSB in high efficiency when 99mTc was tightly bound to the plasmid DNA and this damage could not be completely prevented by DMSO, a free radical scavenger. For the first time, we were able to minimize this effect by increasing the carbon chain lengths between the pyrene moiety and the 99mTc nuclide. However, a critical distance between the 99mTc atom and the DNA helix could not be determined due to the significantly lowered DSB generation resulting from the interaction which is dependent on the type of the 99mTc binding motif. The effect of variable DNA damage caused by the different chain length between the pyrene residue and the Tc-core as well as the possible conformations of the applied Tc-complexes was supplemented with molecular dynamics (MD) calculations. The effectiveness of the DNA-binding 99mTc-labeled pyrene derivatives was demonstrated by comparison to non-DNA-binding 99mTcO4-, since nearly all DNA damage caused by 99mTcO4- was prevented by incubating with DMSO.
The human Sigma1 receptor (S1R), which has been identified as a target with an important role in neuropsychological disorders, was first crystallized 3 years ago. Since S1R structure has no relation ...with another previous crystallized structures, the presence of the new crystal is an important hallmark for the design of agonists and antagonists against this important target. Some years ago, our group identified RC-33, a potent and selective S1R agonist, endowed with neuroprotective properties. In this work, drawing on new structural information, we studied the interactions of RC-33 and its analogs with the S1R binding site by using computational methods such as docking, interaction fingerprints, and receptor-guided alignment three dimensional quantitative structure-activity relationship (3D-QSAR). We found that RC-33 and its analogs adopted similar orientations within S1R binding site, with high similitude with orientations of the crystallized ligands; such information was used for identifying the residues involved in chemical interactions with ligands. Furthermore, the structure-activity relationship of the studied ligands was adequately described considering classical QSAR tests. All relevant aspects of the interactions between the studied compounds and S1R were covered here, through descriptions of orientations, binding interactions, and features that influence differential affinities. In this sense, the present results could be useful in the future design of novel S1R modulators.
Two-pore domain (K2P) channels are twofold symmetric K
channels which control cell excitability by enabling the leak of potassium ions from cells in response to physicochemical stimuli. ...Crystallization of K2P channels revealed the presence of several structural features, which include an external cap. In the available crystallographic structures, the cap is present as non-domain-swapped (NDS) and domain-swapped (DS) chain conformations, where DS chain conformation exchanges two opposing outer helices 180° around the channel. In this work, energy differences between the residues located at the highest point of the cap in NDS and DS conformations were evaluated for TRAAK, a K2P channel that was crystallized in both conformations. Results indicated a preference for DS conformation, but this result is not extensible to TASK K2P channels.