Exploring how changes to the steroidal core alter oleogelation capability in sterol: γ-oryzanol blends.

Oleogels based on sterols such as ß-sitosterol blended with the sterol ester γ-oryzanol are a very interesting class of systems, but there are aspects of their formation and structure that remain elusive. It has previously been shown that a methyl group on the C30 position of the sterol-ester plays an important role in gelation. This work explored the effect that having C30 methyl groups on both the sterol and the sterol-ester had on the gelation process and subsequent gel structure. Lanosterol and saponified γ-oryzanol (which was synthesized as part of this study) were identified as materials of interest, as both feature a methyl group on the C30 position of their steroidal cores. It was observed that both sterols formed gels when blended with γ-oryzanol, and also that lanosterol gelled sunflower oil without the addition of γ-oryzanol. All of these gels were significantly weaker than that formed by ß-sitosterol blended with γ-oryzanol. To explore why, molecular docking simulations along with AFM and SAXS were used to examine these gels on a broad range of length scales. The results suggest that saponified γ-oryzanol-γ-oryzanol gels have a very similar structure to that of ß-sitosterol-γ-oryzanol gels. Lanosterol-γ-oryzanol gels and pure lanosterol gel, however, form with a totally different structure facilitated by the head-to-tail stacking motif exhibited by lanosterol. These results give further evidence that relatively slight changes to the molecular structure of gelators can result in significant differences in subsequent gel properties.

Replacement of heme by soluble guanylate cyclase (sGC) activators abolishes heme-nitric oxide/oxygen (H-NOX) domain structural plasticity.

The gasotransmitter nitric oxide (NO) is a critical endogenous regulator of homeostasis, in major part via the generation of cGMP (cyclic guanosine monophosphate) from GTP (guanosine triphosphate) by NO's main physiological receptor, the soluble guanylate cyclase (sGC). sGC is a heterodimer, composed of an α1 and a ß1 subunit, of which the latter contains the heme-nitric oxide/oxygen (H-NOX) domain, responsible for NO recognition, binding and signal initiation. The NO/sGC/cGMP axis is dysfunctional in a variety of diseases, including hypertension and heart failure, especially since oxidative stress results in heme oxidation, sGC unresponsiveness to NO and subsequent degradation. As a central player in this axis, sGC is the focus of intense research efforts aiming to develop therapeutic molecules that enhance its activity. A class of drugs named sGC "activators" aim to replace the oxidized heme of the H-NOX domain, thus stabilizing the enzyme and restoring its activity. Although numerous studies outline the pharmacology and binding behavior of these compounds, the static 3D models available so far do not allow a satisfactory understanding of the structural basis of sGC's activation mechanism by these drugs. Herein, application NMR describes different conformational states during the replacement of the heme by a sGC activators. We show that the two sGC activators (BAY 58-2667 and BAY 60-2770) significantly decrease the conformational plasticity of the recombinant H-NOX protein domain of Nostoc sp. cyanobacterium, rendering it a lot more rigid compared to the heme-occupied H-NOX. NMR methodology also reveals, for the first time, a surprising bi-directional competition between reduced heme and these compounds, pointing to a highly dynamic regulation of the H-NOX domain. This competitive, bi-directional mode of interaction is also confirmed by monitoring cGMP generation in A7r5 vascular smooth muscle cells by these activators. We show that, surprisingly, heme's redox state impacts differently the bioactivity of these two structurally similar compounds. In all, by NMR-based and functional approaches we contribute unique experimental insight into the dynamic interaction of sGC activators with the H-NOX domain and its dependence on the heme redox status, with the ultimate goal to permit a better design of such therapeutically important molecules.

Backbone and side chain NMR assignments of the H-NOX domain from Nostoc sp. in complex with BAY58-2667 (cinaciguat).

Soluble guanylate cyclase (sGC) enzyme is activated by the gaseous signaling agent nitric oxide (NO) and triggers the conversion of GTP (guanosine 5'-triphosphate) to cGMP (cyclic guanylyl monophosphate). It contains the heme binding H-NOX (heme-nitric oxide/oxygen binding) domain which serves as the sensor of NO and it is highly conserved across eukaryotes and bacteria as well. Many research studies focus on the synthesis of chemical compounds bearing possible therapeutic action, which mimic the heme moiety and activate the sGC enzyme. In this study, we report a preliminary solution NMR (Nuclear Magnetic Resonance) study of the H-NOX domain from Nostoc sp. cyanobacterium in complex with the chemical sGC activator cinaciguat (BAY58-2667). An almost complete sequence-specific assignment of its 1H, 15N and 13C resonances was obtained and its secondary structure predicted by TALOS+.

Molecular Interactions behind the Self-Assembly and Microstructure of Mixed Sterol Organogels.

In this work, we have employed docking and atomistic molecular dynamics (MD) simulations supported by complementary experiments using atomic force microscopy, rheology, and spectroscopy to investigate the self-assembled structure of ß-sitosterol and γ-oryzanol molecules into cylindrical tubules in a nonaqueous solvent. Docking models of several phytosterols, including sitosterol, with oryzanol and other sterol esters demonstrate that for systems to form tubules, the phytosterol sterane group must be stacked in a wedge shape with the ester sterane group and a hydrogen bond must form between the hydroxyl group of the phytosterol and the carbonyl group of the ester. MD of the self-assembled structure were initiated with the molecules in a roughly cylindrical configuration, as suggested from previous experimental studies, and the configurations were found to be stable during 50 ns simulations. We performed MD simulations of two tubules in proximity to better understand the aggregation of these fibrils and how the fibrils interact in order to stick together. We found that an interfibril network of noncovalent bonds, in particular van der Waals and π-π contacts, which is formed between the ferulic acid groups of oryzanol through the hydroxyl, methoxy, and aromatic groups, is responsible for the surface-to-surface interactions between fibrils; an observation supported by molecular spectroscopy. We believe that these interactions are of primary importance in creating a strong organogel network.

Stable emulsions of droplets in a solid edible organogel matrix.

Sitosterol and oryzanol self-assemble to form very firm gels in a range of organic solvents. However, due to the formation of sitosterol hydrate crystals, these gels are unstable in the presence of water, prohibiting the dispersal of water droplets throughout the gel matrix. We demonstrate that by using glycerol as the polar phase rather than water, droplets may be dispersed throughout the oil phase without disrupting the self-assembly of the gel. As increasing volumes of water are added to the glycerol, the G' values decrease. This can be correlated to both a drop in water activity, and also the stability of the fibrils in the presence of glycerol compared to water, as elucidated by molecular dynamics simulations. We explore how changing the total volume of polar droplets, and changing the water content of these droplets alters the strength of 15% w/w sterol gels. We find that gels exhibit G' values of ∼1 × 107 Pa even with ∼30% w/w glycerol dispersed throughout the matrix. At higher glycerol loadings, complex multiple emulsion morphologies can form.

The development of phytosterol-lecithin mixed micelles and organogels.

We demonstrate that by mixing the phytosterol-ester oryzanol with lecithin in an organic solvent, both components may be dispersed at much higher concentrations than they may be individually. Dynamic light scattering and molecular dynamics simulations show that the mechanism for this is the formation of r ∼ 4 nm mixed micelles. Infrared spectroscopy and simulations suggest that these micelles are formed due in part to hydrogen bonding of the phosphate of the lecithin head-group, and the phenol group of the oryzanol. Rheology shows that by mixing these materials at an equimolar ratio, highly viscous suspensions are created. Furthermore, by adding water to these samples, a solid-like gel may be formed which offers mechanical properties close to those desired for a margarine type spread, whilst still solubilizing the oryzanol.

Microstructure of ß-Sitosterol:γ-Oryzanol Edible Organogels.

Rheology and atomic force microscopy (AFM) were employed to examine the microstructure of ß-sitosterol:γ-oryzanol organogels in sunflower oil. Using time-resolved rheology, we followed gel formation, paying specific attention to the fibril aggregation process, which had not been studied in detail previously for this system. Using AFM, we observed gel structures directly and obtained detailed information on the gel structure, far exceeding previous studies. Our analysis suggests that though gels are formed by the self-assembly and aggregation of one-dimensional fibrils, the manner in which these fibrils aggregate into ribbons results in complex structures of higher dimensionality. We emphasize that it is a surprise to find ribbons and not twisted strands. Comparing AFM images of 10% w/w and 20% w/w gelator systems, we observed differences in the degree of branching which are consistent with the rheology. We also observed the individual self-assembled fibrils which make up these gels with much greater clarity than in previous microscopy studies, and the fibril diameters of ∼9.8 nm we measured agree excellently with those obtained from existing small-angle neutron scattering data. These results provide new insight into the structure and formation kinetics of this important organogel system.

SEPIa, a knowledge-driven algorithm for predicting conformational B-cell epitopes from the amino acid sequence.

BACKGROUND: The identification of immunogenic regions on the surface of antigens, which are able to be recognized by antibodies and to trigger an immune response, is a major challenge for the design of new and effective vaccines. The prediction of such regions through computational immunology techniques is a challenging goal, which will ultimately lead to a drastic limitation of the experimental tests required to validate their efficiency. However, current methods are far from being sufficiently reliable and/or applicable on a large scale. RESULTS: We developed SEPIa, a B-cell epitope predictor from the protein sequence, which is sufficiently fast to be applicable on a large scale. The originality of SEPIa lies in the combination of two classifiers, a naïve Bayesian and a random forest classifier, through a voting algorithm that exploits the advantages of both. It is based on 13 sequence-based features, whose values in a 9-residue sequence window are compiled to predict the epitope/non-epitope state of the central residue. The features are related to the type of amino acid, its conservation in homologous proteins, and its tendency of being exposed to the solvent, soluble, flexible, and disordered. The highest signal is obtained from statistical amino acid preferences, but all 13 features contribute non-negligibly in the predictor. SEPIa's average prediction accuracy is limited, with an AUC score (area under the receiver operating characteristic curve) that reaches 0.65 both in 10-fold cross-validation and on an independent test set. It is nevertheless slightly higher than that of other methods evaluated on the same test set. CONCLUSIONS: SEPIa was applied to a test protein whose epitopes are known, human ß2 adrenergic G-protein-coupled receptor, with promising results. Although the actual AUC score is rather low, many of the predicted epitopes cluster together and overlap the experimental epitope region. The reasons underlying the limitations of SEPIa and of all other B-cell epitope predictors are discussed.

Cation-π, amino-π, π-π, and H-bond interactions stabilize antigen-antibody interfaces.

The identification of immunogenic regions on the surface of antigens, which are able to stimulate an immune response, is a major challenge for the design of new vaccines. Computational immunology aims at predicting such regions--in particular B-cell epitopes--but is far from being reliably applicable on a large scale. To gain understanding into the factors that contribute to the antigen-antibody affinity and specificity, we perform a detailed analysis of the amino acid composition and secondary structure of antigen and antibody surfaces, and of the interactions that stabilize the complexes, in comparison with the composition and interactions observed in other heterodimeric protein interfaces. We make a distinction between linear and conformational B-cell epitopes, according to whether they consist of successive residues along the polypeptide chain or not. The antigen-antibody interfaces were shown to differ from other protein-protein interfaces by their smaller size, their secondary structure with less helices and more loops, and the interactions that stabilize them: more H-bond, cation-π, amino-π, and π-π interactions, and less hydrophobic packing; linear and conformational epitopes can clearly be distinguished. Often, chains of successive interactions, called cation/amino-π and π-π chains, are formed. The amino acid composition differs significantly between the interfaces: antigen-antibody interfaces are less aliphatic and more charged, polar and aromatic than other heterodimeric protein interfaces. Moreover, paratopes and epitopes-albeit to a lesser extent-have amino acid compositions that are distinct from general protein surfaces. This specificity holds promise for improving B-cell epitope prediction.

Zn(II)/pyridyloxime complexes as potential reactivators of OP-inhibited acetylcholinesterase: in vitro and docking simulation studies.

In order to investigate the ability of metal complexes to act as reactivators of organophosphorus compounds (OP)-inhibited acetylcholinesterase (AChE), we have synthesized and crystallographically characterized three novel mononuclear Zn(II) complexes formulated as [ZnCl2{(4-py)CHNOH}2] (1), [ZnBr2{(4-py)CHNOH}2] (2) and [Zn(O2CMe)2{(4-py)CHNOH}2]∙2MeCN (3∙2MeCN), where (4-py)CHNOH is 4-pyridinealdoxime. Their reactivation potency was tested in vitro with a slight modification of the Ellman's method using Electric eel acetylcholinesterase and the insecticide paraoxon (diethyl 4-nitrophenyl phosphate) as inhibitor. The activity of the already reported complex [Zn2(O2CPh)2{(4-py)CHNOH}2]·2MeCN (4·2MeCN) and of the clinically used drug obidoxime 1,1'-[oxybis(methylene)]bis{4-[(E)- (hydroxyimino)methyl]pyridinium} was also examined. The results of the in vitro experiments demonstrate moderate reactivation of the metal complexes compared to the drug obidoxime. On the other hand, it is clearly shown that the metal complex is the responsible molecular entity for the observed activity, as the reactivation efficacy of the organic ligand (4-pyridinealdoxime) is found to be inconsequential. Docking simulation studies were performed in the light of predicted complex-enzyme interactions using the paraoxon-inhibited enzyme along with the four Zn(II) complexes and obidoxime as a reference reactivator. The results showed that the three mononuclear metal complexes possess the required characteristics to be accommodated into the active site of AChE, while the entrance of the dinuclear Zn(II) compound is unsuccessful. An interesting outcome of docking simulations is the fact that the mononuclear compounds accommodate into the active site of AChE in a similar mode as obidoxime.

State-of-the-art technology in modern computer-aided drug design.

The quest for small drug-like compounds that selectively inhibit the function of biological targets has always been a major focus in the pharmaceutical industry and in academia as well. High-throughput screening of compound libraries requires time, cost and resources. Therefore, the use of alternative methods is necessary for facilitating lead discovery. Computational techniques that dock small molecules into macromolecular targets and predict the affinity and activity of the small molecule are widely used in drug design and discovery, and have become an integral part of the industrial and academic research. In this review, we present an overview of some state-of-the-art technologies in modern drug design that have been developed for expediting the search for novel drug candidates.

A computational investigation on the role of glycosylation in the binding of alpha1 nicotinic acetylcholine receptor with two alpha-neurotoxins.

Based on the crystal structure of the extracellular domain (ECD) of the mouse nicotinic acetylcholine receptor (nAChR) alpha1 subunit bound to α-bungarotoxin (α-Btx) we have generated in silico models of the human nAChR α1 bound to α-Btx and α-cobratoxin (α-Cbtx), both in the presence and in the absence of the N-linked carbohydrate chain. To gain further insight into the structural role of glycosylation molecular dynamics (MD) simulations were carried out in explicit solvent so as to compare the conformational dynamics of the binding interface between nAChR α1 and the two toxins. An interesting observation during the course of the MD simulations is the strengthening of the receptor-toxin interaction in the presence of the carbohydrate chain, mediated through a shift in the position of the sugars towards the bound toxin. Critical protein-sugar interactions implicate residues Ser187 and Trp184 of nAChR and Thr6, Ser9, and Thr15 of α-Btx, as well as Thr6 and Pro7 of α-Cbtx. Analysis of the predicted residue-specific intermolecular interactions is intended to inspire biophysical studies on the functional role of glycosylation in the gating mechanism.

Conformational dynamics of the anthrax lethal factor catalytic center.

Anthrax lethal factor (LF) is a zinc-metalloprotease that together with the protective antigen constitutes anthrax lethal toxin, which is the most prominent virulence factor of the anthrax disease. The solution nuclear magnetic resonance and in silico conformational dynamics of the 105 C-terminal residues of the LF catalytic core domain in its apo form are described here. The polypeptide adopts a compact structure even in the absence of the Zn(2+) cofactor, while the 40 N-terminal residues comprising the metal ligands and residues that participate in substrate and inhibitor recognition exhibit more flexibility than the C-terminal region.

Purification and biophysical characterization of the core protease domain of anthrax lethal factor.

Anthrax lethal toxin (LeTx) stands for the major virulence factor of the anthrax disease. It comprises a 90kDa highly specific metalloprotease, the anthrax lethal factor (LF). LF possesses a catalytic Zn(2+) binding site and is highly specific against MAPK kinases, thus representing the most potent native biomolecule to alter and inactivate MKK [MAPK (mitogen-activated protein kinase) kinases] signalling pathways. Given the importance of the interaction between LF and substrate for the development of anti-anthrax agents as well as the potential treatment of nascent tumours, the analysis of the structure and dynamic properties of the LF catalytic site are essential to elucidate its enzymatic properties. Here we report the recombinant expression and purification of a C-terminal part of LF (LF(672-776)) that harbours the enzyme's core protease domain. The biophysical characterization and backbone assignments ((1)H, (13)C, (15)N) of the polypeptide revealed a stable, well folded structure even in the absence of Zn(2+), suitable for high resolution structural analysis by NMR.

A computational approach to the study of the binding mode of dual ACE/NEP inhibitors.

Combined blockade of the renin-angiotensin-aldosterone system (RAAS) is an attractive therapeutic strategy for the treatment of cardiovascular diseases. Vasopeptidase inhibitors are a group of compounds capable of inhibiting more than one enzyme, which leads to potentiation of natriuretic peptide actions and suppression of the RAAS. In this study, molecular modeling has been used to elucidate key structural features that govern the binding and/or selectivity of a single compound toward the zinc catalytic sites of the N- and C-domains of the angiotensin-converting enzyme (ACE) and the neutral endopeptidase (NEP). Eleven dual inhibitors were categorized in three classes, according to their zinc binding groups. Analysis of their docked conformers revealed the molecular environment of the catalytic sites and the specific interactions between the inhibitors and amino acid residues that are important for selectivity and cooperativity. In addition, inhibitors were predicted to bind to the C-domain of the ACE with greater affinity than the N-domain, with an average difference in the free energy of binding approximately 2-3 kcal mol(-1). Residues that were identified to actively participate in the binding and stabilizing of the enzyme-inhibitor complexes were analyzed in a consensus way for both the ACE and the NEP. These atomic-level insights into enzyme-ligand binding can be used to drive new structure-based drug design processes in the quest for more selective and effective vasopeptidase inhibitors.

Study of a lipophilic captopril analogue binding to angiotensin I converting enzyme.

Human ACE is a central component of the renin-angiotensin system and a major therapeutic target for cardiovascular diseases. The somatic form of the enzyme (sACE) comprises two homologous metallopeptidase domains (N and C), each bearing a zinc active site with similar but distinct substrate and inhibitor specificities. In this study, we present the biological activity of silacaptopril, a silylated analogue of captopril, and its binding affinity towards ACE. Based on the recently determined crystal structures of both the ACE domains, a series of docking calculations were carried out in order to study the structural characteristics and the binding properties of silacaptopril and its analogues with ACE.

Insights into the anthrax lethal factor-substrate interaction and selectivity using docking and molecular dynamics simulations.

The anthrax toxin of the bacterium Bacillus anthracis consists of three distinct proteins, one of which is the anthrax lethal factor (LF). LF is a gluzincin Zn-dependent, highly specific metalloprotease with a molecular mass of approximately 90 kDa that cleaves most isoforms of the family of mitogen-activated protein kinase kinases (MEKs/MKKs) close to their amino termini, resulting in the inhibition of one or more signaling pathways. Previous studies on the crystal structures of uncomplexed LF and LF complexed with the substrate MEK2 or a MKK-based synthetic peptide provided structure-activity correlations and the basis for the rational design of efficient inhibitors. However, in the crystallographic structures, the substrate peptide was not properly oriented in the active site because of the absence of the catalytic zinc atom. In the current study, docking and molecular dynamics calculations were employed to examine the LF-MEK/MKK interaction along the catalytic channel up to a distance of 20 A from the zinc atom. This residue-specific view of the enzyme-substrate interaction provides valuable information about: (i) the substrate selectivity of LF and its inactivation of MEKs/MKKs (an issue highly important not only to anthrax infection but also to the pathogenesis of cancer), and (ii) the discovery of new, previously unexploited, hot-spots of the LF catalytic channel that are important in the enzyme/substrate binding and interaction.

Soluble guanylyl cyclase activation by HMR-1766 (ataciguat) in cells exposed to oxidative stress.

Many vascular diseases are characterized by increased levels of ROS that destroy the biological activity of nitric oxide and limit cGMP formation. In the present study, we investigated the cGMP-forming ability of HMR-1766 in cells exposed to oxidative stress. Pretreatment of smooth muscle cells with H(2)O(2) reduced cGMP production stimulated by sodium nitroprusside (SNP) or BAY 41-2272. However, pretreatment with H(2)O(2) significantly increased HMR-1766 responses. Similar results were obtained with SIN-1, menadione, and rotenone. In addition, HMR-1766 was more effective in stimulating heme-free sGC compared with the wild-type enzyme. Interestingly, in cells expressing heme-free sGC, H(2)O(2) inhibited instead of potentiated HMR-1766 responses, suggesting that the ROS-induced enhancement of cGMP formation was heme dependent. Moreover, using truncated forms of sGC, we observed that the NH(2)-terminus of the beta(1)-subunit is required for the action of HMR-1766. Finally, to study tolerance development to HMR-1766, cells were pretreated with this sGC activator and reexposed to HMR-1766 or SNP. Results from these experiments demonstrated lack of tolerance development to HMR-1766 as well as lack of cross-tolerance with SNP. We conclude that HMR-1766 is an improved sGC activator as it has the ability to activate oxidized/heme-free sGC and is resistant to the development of tolerance; these observations make HMR-1766 a promising agent for treating diseases associated with increased vascular tone combined with enhanced ROS production.

Low molecular weight inhibitors of the protease anthrax lethal factor.

Anthrax Lethal Factor (LF) is a zinc-dependent metalloprotease that together with the protective antigen constitute the anthrax lethal toxin, the most prominent virulence factor of the disease anthrax. This review summarizes the current knowledge on anthrax toxicity and defense in relation to LF. Particular emphasis is placed on the structural aspects of LF, the properties of its substrates and the achievements in the design of low molecular weight inhibitors of the catalytic activity of the metalloenzyme.