Molecular dynamics simulations to determine the optimal loop length in the helix-loop-helix motif

Novel GABA receptor pesticide targets

The γ-aminobutyric acid (GABA) receptor has four distinct but overlapping and coupled targets of pesticide action importantly associated with little or no cross-resistance. The target sites are differentiated by binding assays with specific radioligands, resistant strains, site-directed mutagenesis and molecular modeling. Three of the targets are for non-competitive antagonists (NCAs) or channel blockers of widely varied chemotypes. The target of the first generation (20th century) NCAs differs between the larger or elongated compounds (NCA-IA) including many important insecticides of the past (cyclodienes and polychlorocycloalkanes) or present (fiproles) and the smaller or compact compounds (NCA-IB) highly toxic to mammals and known as cage convulsants, rodenticides or chemical threat agents. The target of greatest current interest is designated NCA-II for the second generation (21st century) of NCAs consisting for now of isoxazolines and meta-diamides. This new and uniquely different NCA-II site apparently differs enough between insects and mammals to confer selective toxicity. The fourth target is the avermectin site (AVE) for allosteric modulators of the chloride channel. NCA pesticides vary in molecular surface area and solvent accessible volume relative to avermectin with NCA-IBs at 20-22%, NCA-IAs at 40-45% and NCA-IIs at 57-60%. The same type of relationship relative to ligand-docked length is 27-43% for NCA-IBs, 63-71% for NCA-IAs and 85-105% for NCA-IIs. The four targets are compared by molecular modeling for the Drosophila melanogaster GABA-R. The principal sites of interaction are proposed to be: pore V1' and A2' for NCA-IB compounds; pore A2', L6' and T9' for NCA-IA compounds; pore T9' to S15' in proximity to M1/M3 subunit interface (or alternatively an interstitial site) for NCA-II compounds; and M1/M3, M2 interfaces for AVE. Understanding the relationships of these four binding sites is important in resistance management and in the discovery and use of safe and effective pest control agents.

Meta-diamide insecticides acting on distinct sites of RDL GABA receptor from those for conventional noncompetitive antagonists

The RDL GABA receptor is an attractive target of insecticides. Here we demonstrate that meta-diamides [3-benzamido-N-(4-(perfluoropropan-2-yl)phenyl)benzamides] are a distinct class of RDL GABA receptor antagonists showing high insecticidal activity against Spodoptera litura. We also suggest that the mode of action of the meta-diamides is distinct from that of conventional noncompetitive antagonists (NCAs), such as fipronil, picrotoxin, lindane, dieldrin, and α-endosulfan. Using a membrane potential assay, we examined the effects of the meta-diamide 3-benzamido-N-(2-bromo-4-(perfluoropropan-2-yl)-6-(trifluoromethyl)phenyl)-2-fluorobenzamide (meta-diamide 7) and NCAs on mutant Drosophila RDL GABA receptors expressed in Drosophila Mel-2 cells. NCAs had little or no inhibitory activity against at least one of the three mutant receptors (A2'S, A2'G, and A2'N), which were reported to confer resistance to NCAs. In contrast, meta-diamide 7 inhibited all three A2' mutant receptors, at levels comparable to its activity with the wild-type receptor. Furthermore, the A2'S·T6'V mutation almost abolished the inhibitory effects of all NCAs. However, meta-diamide 7 inhibited the A2'S・T6'S mutant receptor at the same level as its activity with the wild-type receptor. In contrast, a G336M mutation in the third transmembrane domain of the RDL GABA receptor abolished the inhibitory activities of meta-diamide 7, although the G336M mutation had little effect on the inhibitory activities of conventional NCAs. Molecular modeling studies also suggested that the binding site of meta-diamides was different from those of NCAs. Meta-diamide insecticides are expected to be prominent insecticides effective against A2' mutant RDL GABA receptors with a different mode of action.

Mutations of the GABA-A receptor alpha1 subunit M1 domain reveal unexpected complexity for modulation by neuroactive steroids

Neuroactive steroids are among the most efficacious modulators of the mammalian GABA-A receptor. Previous work has proposed that receptor potentiation is mediated by steroid interactions with a site defined by the residues alpha1Asn407/Tyr410 in the M4 transmembrane domain and residue alpha1Gln241 in the M1 domain. We examined the role of residues in the alpha1 subunit M1 domain in the modulation of the rat alpha1beta2gamma2L GABA-A receptor by neuroactive steroids. The data demonstrate that the region is critical to the actions of potentiating neuroactive steroids. Receptors containing the alpha1Q241W or alpha1Q241L mutations were insensitive to (3alpha,5alpha)-3-hydroxypregnan-20-one (3alpha5alphaP), albeit with different underlying mechanisms. The alpha1Q241S mutant was potentiated by 3alpha5alphaP, but the kinetic mode of potentiation was altered by the mutation. It is noteworthy that the alpha1Q241L mutation had no effect on channel potentiation by (3alpha,5alpha)-3-hydroxymethyl-pregnan-20-one, but mutation of the neighboring residue, alpha1Ser240, prevented channel modulation. A steroid lacking an H-bonding group on C3 (5alpha-pregnan-20-one) potentiated the wild-type receptor but not the alpha1Q241L mutant. The findings are consistent with a model in which the alpha1Ser240 and alpha1Gln241 residues shape the surface to which steroid molecules bind.

Structural Analysis of Human Lysozyme Using Molecular Dynamics Simulations

In this study, various molecular dynamics simulations were conducted to investigate the effects of ethanol and temperature on the conformational changes of human lysozyme, which may lead insights into amyloidosis. The analyses of some important structural characteristics, such as backbone root-mean-square deviation, secondary structural stability, radius of gyration, accessible surface area, and hydrophobic contact of the hydrophobic core all show that ethanol tends to destabilize human lysozyme at high temperatures. It can be attributed to that higher temperatures result in the destruction of the native structure of this protein, leading to the exposure of the interior hydrophobic core. At this stage, ethanol plays a role to destroy this region by forming hydrophobic interactions between protein and solvent due to its lower polarity comparing to water. Such newly formed intermolecular interactions accelerate the unfolding of this protein, starting from the core between the alpha- and beta-domains. Our results are in good agreement with the previous hypothesis suggesting that the distortion of the hydrophobic core at the alpha- and beta-interface putatively results in the formation of the initial "seed" for amyloid fibril. Although the present results cannot directly be linked to fibril formation, they still provide valuable insights into amyloidosis of human lysozyme.

Recent developments in structural proteomics for protein structure determination

The major challenges in structural proteomics include identifying all the proteins on the genome-wide scale, determining their structure-function relationships, and outlining the precise three-dimensional structures of the proteins. Protein structures are typically determined by experimental approaches such as X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. However, the knowledge of three-dimensional space by these techniques is still limited. Thus, computational methods such as comparative and de novo approaches and molecular dynamic simulations are intensively used as alternative tools to predict the three-dimensional structures and dynamic behavior of proteins. This review summarizes recent developments in structural proteomics for protein structure determination; including instrumental methods such as X-ray crystallography and NMR spectroscopy, and computational methods such as comparative and de novo structure prediction and molecular dynamics simulations.

Molecular dynamics simulations of various coronavirus main proteinases

In this study, two homology models (denoted as MproST and MproSH) of main proteinase (Mpro) from the novel coronavirus associated with severe acute respiratory syndrome (SARS-CoV) were constructed based on the crystal structures of Mpro from transmissible gastroenteritis coronavirus (TGEV) (MproT) and human coronavirus HcoV-229E (MproH), respectively. Both MproST and MproSH exhibit similar folds as their respective template proteins. These homology models reveal three distinct functional domains as well as an intervening loop connecting domains II and III as found in both template proteins. A catalytic cleft containing the substrate binding sites S1 and S2 between domains I and II are also observed. S2 undergoes more significant structural fluctuation than S1 during the 400 ps molecular dynamics simulations because it is located at the open mouth of the catalytic cleft, while S1 is situated in the very bottom of this cleft. The thermal unfolding of these proteins begins at domain III, where the structure is least conserved among these proteins. Mpro may still maintain its proteolytic activity while it is partially unfolded. The electrostatic interaction between Arg40 and Asp186 plays an important role in maintaining the structural integrity of both S1 and S2.

The Stabilizing Effects ofO-glycosylation on the Secondary Structural Integrity of the Designed α-loop-α motif by Molecular Dynamics Simulations

In this study, various 400 ps molecular dynamics simulations were conducted to determine the stabilizing effect of O-glycosylation on the secondary structural integrity of the design alpha-loop-alpha motif, which has the optimal loop length of 7 Gly residues (denoted as N-A16G7A16-C). In general, O-glycosylation stabilizes the structural integrity of the model peptide regardless of the length and position of glycosylation sites because it decreases the opportunity for water molecules to compete for the intramolecular hydrogen bonds. The designed peptide exhibits the highest helicity when residues 11 and 31 are replaced with Ser residues followed by O-linked with 3 galactose residues, representing the "face-to-face" glycosylation near the loop. In this case, the loop exhibits an extended conformation and several new hydrogen bonds are observed between the main chain of the loop and the galactose residues, resulting in decreasing the fluctuation and increasing the stability of the entire peptide. When the glycosylation are made close to the loop, the secondary structural integrity of the alpha-loop-alpha motif increases with the number of galactose residues. In addition, "face-to-face" glycosylation increases the structural integrity of this motif to a greater extent than "back-to-back" glycosylation. However, when the glycosylation are created away from the loop and near the N- and C-termini, no general rule is found for the stabilizing effect.

Comparative modeling of a GABAA alpha1 receptor using three crystal structures as templates

We built a model of a GABAA alpha1 receptor (GABAAR) that combines the ligand binding (LBD) and the transmembrane domains (TMD). We used six steps: (1) a four-alpha helical bundle in the crystal structure of bovine cytochrome c oxidase (2OCC) was identified as a template for the TMD of a single subunit. (2) The five pore-forming alpha helices of a bacterial mechanosensitive channel (1MSL) served as a template for the pentameric ion channel. (3) Five copies of the tetrameric template from 2OCC were superimposed on 1MSL to produce a homopentamer containing 20 alpha helices arranged around a funnel-shaped central pore. (4) Five copies of the GABAAR sequence were threaded onto the alpha-helical segments of this template and inter-helical loops were generated to produce the TMD model. (5) A model of the LBD was built by threading the aligned sequence of GABAAR onto the crystal structure of the acetylcholine binding protein (1I9B). (6) The models of the LBD and the TMD were aligned along a common five-fold axis, moved together along that axis until in vdW contact, merged, and then optimized with restrained molecular dynamics. Our model corresponds closely with recently published coordinates of the acetylcholine receptor (1OED) but also explains additional features. Our model reveals structures of loops that were not visible in the cryoelectron micrograph and satisfies most labeling and mutagenesis data. It also suggests mechanisms for ligand binding transduction, ion selectivity, and anesthetic binding.