In silico study of subtilisin-like protease 1 (SUB1) from different Plasmodium species in complex with peptidyl-difluorostatones and characterization of potent pan-SUB1 inhibitors

Plasmodium falciparum subtilisin-like protease 1 (SUB1) is a novel target for the development of innovative antimalarials. We recently described the first potent difluorostatone-based inhibitors of the enzyme ((4S)-(N-((N-acetyl-l-lysyl)-l-isoleucyl-l-threonyl-l-alanyl)-2,2-difluoro-3-oxo-4-aminopentanoyl)glycine (1) and (4S)-(N-((N-acetyl-l-isoleucyl)-l-threonyl-l-alanylamino)-2,2-difluoro-3-oxo-4-aminopentanoyl)glycine (2)). As a continuation of our efforts towards the definition of the molecular determinants of enzyme-inhibitor interaction, we herein propose the first comprehensive computational investigation of the SUB1 catalytic core from six different Plasmodium species, using homology modeling and molecular docking approaches. Investigation of the differences in the binding sites as well as the interactions of our inhibitors 1,2 with all SUB1 orthologues, allowed us to highlight the structurally relevant regions of the enzyme that could be targeted for developing pan-SUB1 inhibitors. According to our in silico predictions, compounds 1,2 have been demonstrated to be potent inhibitors of SUB1 from all three major clinically relevant Plasmodium species (P. falciparum, P. vivax, and P. knowlesi). We next derived multiple structure-based pharmacophore models that were combined in an inclusive pan-SUB1 pharmacophore (SUB1-PHA). This latter was validated by applying in silico methods, showing that it may be useful for the future development of potent antimalarial agents.     

In silico and in vitro screening to identify structurally diverse non-azole CYP51 inhibitors as potent antifungal agent

The problem of resistance to azole class of antifungals is a serious cause of concern to the medical fraternity and thus there is an urgent need to identify non-azole scaffolds with high affinity for lanosterol 14α-demethylase (CYP51). In view of this we have attempted to identify novel non-azole CYP51 inhibitors through the application of pharmacophore based virtual screening and in vitro evaluation. A rigorously validated pharmacophore model comprising of 2 hydrogen bond acceptor and 2 hydrophobic features has been developed and used to mine NCI database. Out of 265 retrieved hits, NSC 1215 and 1520 have been chosen on the basis of Lipinski’s rule of five, fit and estimated values. Both the hits were docked into the active site of CYP51. In view of high fit value and CDocker score, NSC 1215 and 1520 have been subjected to in vitro microbiological assay. The result reveals that NSC 1215 and 1520 are active against Candida albicans, Candida parapsilosis, Candida tropicalis, and Aspergillus niger. In addition to this the absorption characteristics of both the hits have also been determined using the rat sac technique and permeation in order of NSC 1520 > NSC 1215 has been observed.     

Comparison and analysis of the structures and binding modes of antifungal SE and CYP51 inhibitors

With the abuse of clinical broad-spectrum antimicrobial agents, immunosuppressive agents, chemotherapy drugs, the emergence of pathogenic fungi resistance is more and more frequent. However, there is still no effective treatment for the fungal resistance. Squalenee epoxidase (SE) and 14 α-demethylase (CYP51) are important antifungal drug targets. In order to achieve a deeper insight into the structural characteristics and the action modes of SE and CYP51inhibitors, the homology model of SE (Candida albicans) was constructed using monooxygenase of Pseudomonas aeruginosa as template, and the reliability of model was confirmed by Ramachandran plots and Verify 3D. Subsequently, the molecular superimposition and molecular docking were performed, and the pharmacophore model based on the CYP51 receptor structure was constructed. The results indicate that SE and CYP51 inhibitors have common structural feature with two parts of essential fragments, which are mainly composed of aromatic groups. In addition, the fragment structures of inhibitors are combined in the similar hydrophobic pockets through the hydrophobic forces. The present study provides a deeper perspective to understand the characteristics and docking modes of SE and CYP51 inhibitors. It can be used to guide the optimization and design of SE and CYP51 inhibitors. In addition, it also provides the oretical support for the development of dual target antifungal inhibitors (SE and CYP51), which can help us solve the problem of fungi resistance.     

Structural insights into the camel milk lactoperoxidase: Homology modeling and molecular dynamics simulation studies

Lactoperoxodase (LPO) is a heme peroxidase enzyme present in mammalian milk. It is an antimicrobial protein with wide range of industrial applications. Although the three dimensional structure of LPO from various mammalian species has been reported, but its structure from camel source is still unknown. So far, the crystallization attempts have not been successful in determining camel LPO (cLPO) structure. Herein, we developed the three dimensional structure of cLPO by homology modeling approach using prime module available in Schrodinger suite. The developed model in complex with ligand hypothiocyanate (OSCN⁻) was further validated by Ramachandran plot followed by molecular dynamics (MD) simulation studies using Desmond module of Schrodinger. cLPO model exhibited overall structural similarity with template crystal structure, however, it displayed different interaction pattern of amino acid residues with ligand OSCN⁻ in comparison to template crystal structure. Moreover, the ligand binding site environment in cLPO is more polar, less hydrophobic, and harbours more number of charged residues than template crystal structure. The substrate binding pocket environment of cLPO shows a considerable difference from template crystal structure. This subsequently resulted in dissimilar behaviour of ligand during the course of MD simulation studies.     

Prediction of dual agents as an activator of mutant p53 and inhibitor of Hsp90 by docking,molecular dynamic simulation and virtual screening

Heat shock protein90s (Hsp90s) play a crucial role in the development of cancer, and their inhibitors are a main target for tumor suppression. P53 also is a tumor suppressor, but in cancer cells, mutations in the p53 gene lead to the inactivation and accumulation of protein. For instance, the ninth p53 cancer mutation, Y220C, destabilizes the p53 core domain. Small molecules have been assumed to bind to Y220C DNA-binding domain and reactivate cellular mutant p53 functions. In this study, one of the mutant p53 activators is suggested as an Hsp90 inhibitor according to a pyrazole scaffold. To confirm a new ligand as a dual agent, molecular docking and molecular dynamic simulations were performed on both proteins (p53 and Hsp90). Molecular dynamic simulations were also conducted to evaluate the obtained results on the other two pyrazole structures, one known as Hsp90 inhibitor and the other as the reported mutant p53 activator. The findings indicate that the new ligand was stable in the active site of both proteins. Finally, a virtual screening was performed on ZINC database, and a set of new dual agents was proposed according to the new ligand scaffold.     

Identification and characterization of a glycosaminoglycan binding site on interleukin-10 via molecular simulation methods

The biological function of the pleiotropic cytokine interleukin-10 (IL-10), which has an essential role in inflammatory processes, is known to be affected by glycosaminoglycans (GAGs). GAGs are essential constituents of the extracellular matrix with an important role in modulating the biological function of many proteins. The molecular mechanisms governing the IL-10–GAG interaction, though, are unclear so far. In particular, detailed knowledge about GAG binding sites and recognition mode on IL-10 is lacking, despite of its imminent importance for understanding the functional consequences of IL-10–GAG interaction. In the present work, we report a GAG binding site on IL-10 identified by applying computational methods based on Coulomb potential calculations and specialized molecular dynamics simulations. The identified GAG binding site is constituted of several positively charged residues, which are conserved among species. Exhaustive conformational space sampling of a series of GAG ligands binding to IL-10 led to the observation of two GAG binding modes in the predicted binding site, and to the identification of IL-10 residues R104, R106, R107, and K119 as being most important for molecular GAG recognition. In silico mutation as well as single-residue energy decomposition and detailed analysis of hydrogen-bonding behavior led to the conclusion that R107 is most essential and assumes a unique role in IL-10–GAG interaction. This structural and dynamic characterization of GAG-binding to IL-10 represents an important step for further understanding the modulation of the biological function of IL-10.     

Computational docking,molecular dynamics simulation and subsite structure analysis of a maltogenic amylase from Bacillus lehensis G1 provide insights into substrate and product specificity

Maltogenic amylase (MAG1) from Bacillus lehensis G1 displayed the highest hydrolysis activity on β-cyclodextrin (β-CD) to produce maltose as a main product and exhibited high transglycosylation activity on malto-oligosaccharides with polymerization degree of three and above. These substrate and product specificities of MAG1 were elucidated from structural point of view in this study. A three-dimensional structure of MAG1 was constructed using homology modeling. Docking of β-CD and malto-oligosaccharides was then performed in the MAG1 active site. An aromatic platform in the active site was identified which is responsible in substrate recognition especially in determining the enzyme’s preference toward β-CD. Molecular dynamics (MD) simulation showed MAG1 structure is most stable when docked with β-CD and least stable when docked with maltose. The docking analysis and MD simulation showed that the main subsites for substrate stabilization in the active site are −2, −1, +1 and +2. A bulky residue, Trp359 at the +2 subsite was identified to cause steric interference to the bound linear malto-oligosaccharides thus prevented it to occupy subsite +3, which can only be reached by a highly bent glucose molecule such as β-CD. The resulted modes of binding from docking simulation show a good correlation with the experimentally determined hydrolysis pattern. The subsite structure generated from this study led to a possible mode of action that revealed how maltose was mainly produced during hydrolysis. Furthermore, maltose only occupies subsite +1 and +2, therefore could not be hydrolyzed or transglycosylated by the enzyme. This important knowledge has paved the way for a novel structure-based molecular design for modulation of its catalytic activities.     

Combined in silico approaches for the identification of novel inhibitors of human islet amyloid polypeptide (hIAPP) fibrillation

Human islet amyloid polypeptide (hIAPP) is a natively unfolded polypeptide hormone of glucose metabolism, which is co-secreted with insulin by the β-cells of the pancreas. In patients with type 2 diabetes, IAPP forms amyloid fibrils because of diabetes-associated β-cells dysfunction and increasing fibrillation, in turn, lead to failure of secretory function of β-cells. This provides a target for the discovery of small organic molecules against protein aggregation diseases. However, the binding mechanism of these molecules with monomers, oligomers and fibrils to inhibit fibrillation is still an open question. In this work, ligand and structure-based in silico approaches were used to identify novel fibrillation inhibitors and/or fibril binding compounds. The best pharmacophore model was used as a 3D search query for virtual screening of a compound database to identify novel molecules having the potential to be therapeutic agents against protein aggregation diseases. Docking and molecular dynamics simulation studies were used to explore the interaction pattern and mechanism of the identified novel small molecules with predicted hIAPP structure, its aggregation prone conformation and fibril forming segments. We show that catechins with galloyl group and molecules having two to three planar apolar rings bind to hIAPP structures and fibril forming segments with greater affinity. The differences in binding affinities of different compounds against several fibril forming segments of the peptide suggest that a mixture of active compounds may be required for treatment of aggregation diseases.     

Comparative modeling and molecular docking analysis of white,brown and soft rot fungal laccases using lignin model compounds for understanding the structural and functional properties of laccases

Extrinsic catalytic properties of laccase enable it to oxidize a wide range of aromatic (phenolic and non-phenolic) compounds which makes it commercially an important enzyme. In this study, we have extensively compared and analyzed the physico-chemical, structural and functional properties of white, brown and soft rot fungal laccases using standard protein analysis software. We have computationally predicted the three-dimensional comparative models of these laccases and later performed the molecular docking studies using the lignin model compounds. We also report a customizable rapid and reliable protein modelling and docking pipeline for developing structurally and functionally stable protein structures. We have observed that soft rot fungal laccases exhibited comparatively higher structural variation (higher random coil) when compared to brown and white rot fungal laccases. White and brown rot fungal laccase sequences exhibited higher similarity for conserved domains of Trametes versicolor laccase, whereas soft rot fungal laccases shared higher similarity towards conserved domains of Melanocarpus albomyces laccase. Results obtained from molecular docking studies showed that aminoacids PRO, PHE, LEU, LYS and GLN were commonly found to interact with the ligands. We have also observed that white and brown rot fungal laccases showed similar docking patterns (topologically monomer, dimer and trimer bind at same pocket location and tetramer binds at another pocket location) when compared to soft rot fungal laccases. Finally, the binding efficiencies of white and brown rot fungal laccases with lignin model compounds were higher compared to the soft rot fungi. These findings can be further applied in developing genetically efficient laccases which can be applied in growing biofuel and bioremediation industries.     

The role of LasR active site amino acids in the interaction with the Acyl Homoserine Lactones (AHLs) analogues: A computational study

The present work combines molecular docking calculations, 3D-QSAR, molecular dynamics simulations and free binding energy calculations (MM/PBSA and MM/GBSA) in a set of 28 structural analogues of acyl homoserine lactones with Quorum Sensing antagonist activity. The aim of this work is to understand how ligand binds and is affected by the molecular microenvironment in the active site of the LasR receptor for pseudomonas aeruginosa. We also study the stability of the interaction to find key structural characteristics that explain the antagonist activities of this set of ligands. This information is relevant for the rational modification or design of molecules and their identification as powerful LasR modulators.The analysis of molecular docking simulations shows that the 28 analogues have a similar binding mode compared to the native ligand. The carbonyl groups belonging to the lactone ring and the amide group of the acyl chain are oriented towards the amino acids forming hydrogen bond like interactions. The difference in antagonist activity is due to location and orientation of the LasR side chains within the hydrophobic pocket in its binding site. Additionally, we carried out molecular dynamics simulations to understand the conformational changes in the ligand-receptor interaction and the stability of each complex. Results show a direct relationship among the interaction energies of the ligands and the activities as an antagonist of the LasR receptor.     

Effects of synergistic and non-synergistic anions on the iron binding site from serum transferrin: A molecular dynamic simulation analysis

The role of the anions in transferrin chemistry highlights the importance of the anion binding site in transferrin family. A synergistic anion as carbonate is an anion that is required for iron binding by transferrin while non-synergistic anions do not act as the synergistic anions to promote iron binding, but affect the iron binding and release. Some questions remain unclear about the difference between synergistic and non-synergistic anion functions. In the present work, molecular dynamic simulation techniques were employed in order to gain access into a molecular level understanding of the iron binding site of the human serum transferrin during the synergistic and non-synergistic anion binding. For this purpose, a comparative analysis was performed to illustrate the observed changes. In addition to the comparison between the synergistic and non-synergistic anions, structural differences between two synergistic anions, Carbonate and Oxalate were studied. Meanwhile,the simulation of the open (Apo), partially closed (Carbonate) and fully closed (Carbonate-Fe) forms of the transferrin structure allows a direct comparison between the iron binding site of these three states.On the basis of results, synergistic anions form high affinity binding site, while non-synergistic anions act like Apo state of the transferrin structure and change the proper conformation of the binding site. In order to act as a synergistic anion and form high affinity binding site, anion stereochemistry and interactions must be able to achieve a Carbonate-like configuration. Carbonate complex showed the highest binding affinity and electrostatic energy is the major favorable contributor to synergistic anion-transferrin interaction. Carbonate and Oxalatecomplexes as synergistic anions have many features in common, without a significant change in the transferrin structure. Only the residues in the vicinity of the binding site showed a little different conformation depending on whether the synergistic anion is Carbonate orOxalate.Finally, the results show thatASP63, GLY65 and HIS249 residues have the maximum displacement during the Carbonate and iron binding. ASP63 and HIS249 are the residues, which are coordinated to the iron and GLY65 is in the second shell residuesof the transferrin structure.     

Ole e 13 is the unique food allergen in olive: Structure-functional,substrates docking,and molecular allergenicity comparative analysis

Thaumatin-like proteins (TLPs) are enzymes with important functions in pathogens defense and in the response to biotic and abiotic stresses. Last identified olive allergen (Ole e 13) is a TLP, which may also importantly contribute to food allergy and cross-allergenicity to pollen allergen proteins. The goals of this study are the characterization of the structural-functionality of Ole e 13 with a focus in its catalytic mechanism, and its molecular allergenicity by extensive analysis using different molecular computer-aided approaches covering a) functional-regulatory motifs, b) comparative study of linear sequence, 2-D and 3D structural homology modeling, c) molecular docking with two different β-D-glucans, d) conservational and evolutionary analysis, e) catalytic mechanism modeling, and f) IgE-binding, B- and T-cell epitopes identification and comparison to other allergenic TLPs.Sequence comparison, structure-based features, and phylogenetic analysis identified Ole e 13 as a thaumatin-like protein. 3D structural characterization revealed a conserved overall folding among plants TLPs, with mayor differences in the acidic (catalytic) cleft. Molecular docking analysis using two β-(1,3)-glucans allowed to identify fundamental residues involved in the endo-1,3-β-glucanase activity, and defining E84 as one of the conserved residues of the TLPs responsible of the nucleophilic attack to initiate the enzymatic reaction and D107 as proton donor, thus proposing a catalytic mechanism for Ole e 13. Identification of IgE-binding, B- and T-cell epitopes may help designing strategies to improve diagnosis and immunotherapy to food allergy and cross-allergenic pollen TLPs.     

Computational analysis of Amsacrine resistance in human topoisomerase II alpha mutants (R487K and E571K) using homology modeling,docking and all-atom molecular dynamics simulation in explicit solvent

Amsacrine is an effective topoisomerase II enzyme inhibitor in acute lymphatic leukemia. Previous experimental studies have successfully identified two important mutations (R487K and E571K) conferring 100 and 25 fold resistance to Amsacrine respectively. Although the reduction of the cleavage ligand-DNA-protein ternary complex has been well thought as the major cause of drug resistance, the detailed energetic, structural and dynamic mechanisms remain to be elusive. In this study, we constructed human topoisomerase II alpha (hTop2α) homology model docked with Amsacrine based on crystal structure of human Top2β in complex with etoposide. This wild type complex was used to build the ternary complex with R487K and E571K mutants. Three 500 ns molecular dynamics simulations were performed on complex systems of wild type and two mutants. The detailed energetic, structural and dynamic analysis were performed on the simulation data. Our binding data indicated a significant impairment of Amsacrine binding energy in the two mutants compared with the wild type. The order of weakening (R487K > E571K) was in agreement with the order of experimental drug resistance fold (R489K > E571K). Our binding energy decomposition further indicated that weakening of the ligand-protein interaction rather than the ligand-DNA interaction was the major contributor of the binding energy difference between R487K and E571K. In addition, key residues contributing to the binding energy (ΔG) or the decrease of the binding energy (ΔΔG) were identified through the energy decomposition analysis. The change in ligand binding pose, dynamics of protein, DNA and ligand upon the mutations were thoroughly analyzed and discussed. Deciphering the molecular basis of drug resistance is crucial to overcome drug resistance using rational drug design.     

Comparative analysis of polyspecificity of the endogenous tRNA synthetase of different expression host towards photocrosslinking amino acids using an in silico approach

Photo-induced covalent crosslinking has emerged as the powerful strategy for analyzing and characterizing the protein–protein interaction and mapping protein 3D conformations. In the last decades, a number of photocrosslinking amino acids have been reported but only a few have been efficiently utilized for photocrosslinking purposes. Recently, incorporation of diazirine containing photoactivatable analogs such as photo-methionine, photo-leucine, photo-isoleucine and photo-lysine into target proteins were accomplished in live cells (Human A549cells, HEK 293) by depleting corresponding natural amino acid and supplementing these analogs in the medium. Likewise, incorporation of photo-methionine and photo-leucine is also reported in E. coli. Incorporation of these unnatural amino acids were demonstrated only in a limited number species, thereby conventional methods have been utilized for the protein–protein interaction study in other species. With this in mind, we studied in silico analysis of polyspecificity of four endogenous tRNA synthetases (LeuRS, IleRS, MetRS, and LysRS) from six different species such as Escherichia coli, Pseudomonas fluorescens, Corynebacterium glutamicum, Saccharomyces cerevisiae, Aspergillus oryzae and Homo sapiens towards its photocrosslinking amino acids. In addition, here we describe the active site similarity of different protein bio-factories. Based on the active site similarity and similar binding mode, we predicted that the endogenous tRNA synthetases of all the species are reactive to corresponding photoactivatable analogs. This is the first in silico study to demonstrate that the photocrosslinking unnatural amino acids are recognized by the endogenous tRNA synthetases of different protein expression biofactories.     

Modeling of flux,binding and substitution of urea molecules in the urea transporter dvUT

Urea transporters (UTs) are transmembrane proteins that transport urea molecules across cell membranes and play a crucial role in urea excretion and water balance. Modeling the functional characteristics of UTs helps us understand how their structures accomplish the functions at the atomic level, and facilitates future therapeutic design targeting the UTs. This study was based on the crystal structure of Desulfovibrio vulgaris urea transporter (dvUT). To model the binding behavior of urea molecules in dvUT, we constructed a cooperative binding model. To model the substitution of urea by the urea analogue N,N′-dimethylurea (DMU) in dvUT, we calculated the occupation probability of DMU along the urea pore and the ratio of the occupation probabilities of DMU at the external (S_ext) and internal (S_int) binding sites, and we established the mutual substitution rule for binding and substitution of urea and DMU. Based on these calculations and modelings, together with the use of the Monte Carlo (MC) method, we further modeled the urea flux in dvUT, equilibrium urea binding to dvUT, and the substitution of urea by DMU in the dvUT. Our modeling results are in good agreement with the existing experimental functional data. Furthermore, the modelings have discovered the microscopic process and mechanisms of those functional characteristics. The methods and the results would help our future understanding of the underlying mechanisms of the diseases associated with impaired UT functions and rational drug design for the treatment of these diseases.     

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy

With the present work we quantitatively studied the modellability of the inactive state of Class A G protein-coupled receptors (GPCRs). Specifically, we constructed models of one of the Class A GPCRs for which structures solved in the inactive state are available, namely the β₂ AR, using as templates each of the other class members for which structures solved in the inactive state are also available. Our results showed a detectable linear correlation between model accuracy and model/template sequence identity. This suggests that the likely accuracy of the homology models that can be built for a given receptor can be generally forecasted on the basis of the available templates. We also probed whether sequence alignments that allow for the presence of gaps within the transmembrane domains to account for structural irregularities afford better models than the classical alignment procedures that do not allow for the presence of gaps within such domains. As our results indicated, although the overall differences are very subtle, the inclusion of internal gaps within the transmembrane domains has a noticeable a beneficial effect on the local structural accuracy of the domain in question.     

Molecular modeling study on the effect of residues distant from the nucleotide-binding portion on RNA binding in Staphylococcus aureus Hfq

Hfq is an abundant RNA-binding bacterial protein that was first identified in E. coli as a required host factor for phage Qβ RNA replication. The pleiotrophic phenotype resulting from the deletion of Hfq predicates the importance of this protein. Two RNA-binding sites have been characterized: the proximal site which binds sRNA and mRNA and the distal site which binds poly(A) tails. Previous studies mainly focused on the key residues in the proximal site of the protein. A recent mutation study in E. coli Hfq showed that a distal residue Val43 is important for the protein function. Interestingly, when we analyzed the sequence and structure of Staphylococcus aureus Hfq using the CONSEQ server, the results elicited that more functional residues were located far from the nucleotide-binding portion (NBP). From the analysis seven individual residues Asp9, Leu12, Glu13, Lys16, Gln31, Gly34 and Asp40 were selected to investigate the conformational changes in Hfq–RNA complex due to point mutation effect of those residues using molecular dynamics simulations. Results showed a significant effect on Asn28 which is an already known highly conserved functionally important residue. Mutants D9A, E13A and K16A depicted effects on base stacking along with increase in RNA pore diameter, which is required for the threading of RNA through the pore for the post-translational modification. Further, the result of protein stability analysis by the CUPSAT server showed destabilizing effect in the most mutants. From this study we characterized a series of important residues located far from the NBP and provide some clues that those residues may affect sRNA binding in Hfq.     

Molecular simulations of lactose-bound and unbound forms of the FaeG adhesin reveal critical amino acids involved in sugar binding

F4 fimbriae are protein filaments found in enterotoxigenic Escherichia coli cells and are implicated in the process of bacterial infection due to their function as bacterial adhesins. These filaments are comprised from several proteins, but the bacterial adhesin FaeG, which is a lactose-binding protein, is the major subunit comprising F4 fimbriae. Crystal structures for three variants of the FaeG protein were recently solved, including the ad variant of FaeG that was crystallized in complex with lactose. However, the dynamics of the FaeG protein bound to lactose have not been explored previously using molecular dynamics simulations. Therefore, in order to study the dynamical interactions between the FaeG ad variant and lactose, we have carried out the first all-atom molecular dynamics simulations of this system. We have also probed the role of crystallographic water molecules on the stability of lactose in the FaeG binding site, and have simulated seven FaeG mutants to probe the influence of amino acid substitutions on the ability of FaeG to bind lactose effectively. Our simulations agree well with experimental results for the influence of mutations on lactose binding, provide dynamical insights into the interactions of FaeG with lactose, and also suggest the possibility of additional regions of the FaeG protein that may act as secondary lactose binding sites.     

Space and time basis function design for the method of moments in time‐domain analysis of wire and planar structures

In this article, an algorithm able to deal at the same time with wire frame and surface patch models for the method of moments in time domain is presented. After a unified theory combining both models, attention is focused on stability dependence issues on the time basis function chosen and on other algorithm parameters. An accurate analysis of time interpolation functions and of their time filtering properties is provided. © 2011 Wiley Periodicals, Inc. Int J RF and Microwave CAE, 2011.     

Toward the three-dimensional structure and lysophosphatidic acid binding characteristics of the LPA4/p2y9/GPR23 receptor: A homology modeling study

Lysophosphatidic acid (LPA) is a naturally occurring phospholipid that initiates a broad array of biological processes, including those involved in cell proliferation, survival and migration via activation of specific G protein-coupled receptors located on the cell surface. To date, at least five receptor subtypes (LPA_1–5) have been identified. The LPA_1–3 receptors are members of the endothelial cell differentiation gene (Edg) family. LPA₄, a member of the purinergic receptor family, and the recently identified LPA₅ are structurally distant from the canonical Edg LPA_1–3 receptors. LPA₄ and LPA₅ are linked to G_q, G_12/13 and G_s but not G_i, while LPA_1–3 all couple to G_i in addition to G_q and G_12/13. There is also evidence that LPA₄ and LPA₅ are functionally different from the Edg LPA receptors. Computational modeling has provided useful information on the structure–activity relationship (SAR) of the Edg LPA receptors. In this work, we focus on the initial analysis of the structural and ligand-binding properties of LPA₄, a prototype non-Edg LPA receptor. Three homology models of the LPA₄ receptor were developed based on the X-ray crystal structures of the ground state and photoactivated bovine rhodopsin and the recently determined human β₂-adrenergic receptor. Docking studies of LPA in the homology models were then conducted, and plausible LPA binding loci were explored. Based on these analyses, LPA is predicted to bind to LPA₄ in an orientation similar to that reported for LPA_1–3, but through a different network of hydrogen bonds. In LPA_1–3, the ligand polar head group is reported to interact with residues at positions 3.28, 3.29 and 7.36, whereas three non-conserved amino acid residues, S114(3.28), T187(EL2) and Y265(6.51), are predicted to interact with the polar head group in the LPA₄ receptor models.