Octadecylamine (ODA) was used to intercalate a fine-grained and a coarse-grained fraction of natural Mg-vermiculite (VER) using a low-temperature melting procedure. Mixtures of Mg-vermiculite fractions and powdered ODA in the molar ratios of 2:1, 1:1, 1:2 and 1:6 were homogenized and heated for 1, 3, 15 and 30 h at 80°C to prepare intercalated samples. X-ray powder diffraction analysis of intercalated samples was combined with molecular modeling to investigate their interlayer structure. Significant amounts of non-intercalated vermiculite and diffuse peaks with very low intensity and basal spacings close to 29 Å were identified when the lowest concentration (molar ratio VER:ODA = 2:1) was used. According to molecular modeling, this indicates the initial stage of a one-layer arrangement of distorted ODA molecules in the interlayer. If the concentration of ODA molecules and treatment time were increased, a two-layer arrangement of ODA molecules with a different ODA chain-disorder and interlayer-space saturation was identified. Interlayer ODA molecules were inclined to the vermiculite basal plane with an inclination angle for two-layer arrangements that ranged from 76 to 95°. Experimental basal spacings with these two-layer arrangements varied from 52 to 58 Å and were in agreement with molecular modeling results. A fully-saturated 58 Å two-layer ODA arrangement was identified when higher ODA concentrations (VER:ODA = 1:2 and 1:6) and 15 and 30 h treatment times were used. There was no significant difference between ODA-intercalated samples prepared using fine-grained and coarse-grained Mg-vermiculite fractions. A grafted ODA-chain nano-layer with a 49.6(2.1) Å average height was observed on the surface of thin ODA-intercalated micro-flakes using atomic force microscopy. Grafted ODA chains not only created an homogeneous surface nano-layer, but also variable-width channels between the ODA molecules.

The sugarcane giant borer, Telchin licus licus, is an insect pest that causes significant losses in sugarcane crops and in the sugar-alcohol sector. Chemical and manual control methods are not effective. As an alternative, in the current study, we have screened Bacillus thuringiensis (Bt) Cry toxins with high toxicity against this insect. Bioassays were conducted to determine the activity of four Cry toxins (Cry1A (a, b, and c) and Cry2Aa) against neonate T. licus licus larvae. Notably, the Cry1A family toxins had the lowest LC50 values, in which Cry1Ac presented 2.1-fold higher activity than Cry1Aa, 1.7-fold larger than Cry1Ab, and 9.7-fold larger than Cry2Aa toxins. In silico analyses were performed as a perspective to understand putative interactions between T. licus licus receptors and Cry1A toxins. The molecular dynamics and docking analyses for three putative aminopeptidase N (APN) receptors (TlAPN1, TlAPN3, and TlAPN4) revealed evidence for the amino acids that may be involved in the toxin–receptor interactions. Notably, the properties of Cry1Ac point to an interaction site that increases the toxin's affinity for the receptor and likely potentiate toxicity. The interacting amino acid residues predicted for Cry1Ac in this work are probably those shared by the other Cry1A toxins for the same region of APNs. Thus, the presented data extend the existing knowledge of the effects of Cry toxins on T. licus licus and should be considered in further development of transgenic sugarcane plants resistant to this major occurring insect pest in sugarcane fields.

Herbicides targeting grass plastidic acetyl coenzyme A carboxylase (ACC) are effective selective graminicides. Their intensive use worldwide has selected for resistance genes in a number of grass weed species. Biochemistry and molecular biology have been the means of determining the herbicidal activity and selectivity toward crop plants of ACC-inhibiting herbicides. In recent years, elucidation of the tridimensional structure of ACC and identification of five amino acid residues within the ACC carboxyl transferase domain that are critical determinants for herbicide sensitivity shed light on the basis of ACC-based resistance to herbicides. However, metabolism-based resistance to ACC-inhibiting herbicides is much less well known, although this type of resistance seems to be widespread. A number of genes thus endow resistance to ACC-inhibiting herbicides, with the possibility for various resistance genes that confer dominant resistance at the herbicide field rate to accumulate within a single weed population or plant. This, together with a poor knowledge of the genetic parameters driving resistance, renders the evolution of resistance to ACC-inhibiting herbicides unpredictable. Future research should consider developing tactics to slow the spread of resistance. For this purpose, it is crucial that our understanding of metabolism-based resistance improves rapidly because this mechanism is complex and can confer resistance to herbicides with different target sites.

Herbicides that target the enzyme acetolactate synthase (ALS) are among the most widely used in the world. Unfortunately, these herbicides are also notorious for their ability to select resistant (R) weed populations. Now, there are more weed species that are resistant to ALS-inhibiting herbicides than to any other herbicide group. In most cases, resistance to ALS-inhibiting herbicides is caused by an altered ALS enzyme. The frequent occurrence of weed populations resistant to ALS inhibitors can be attributed to the widespread usage of these herbicides, how they have been used, the strong selection pressure they exert, and the resistance mechanism. In several cropping systems, ALS-inhibiting herbicides were used repeatedly as the primary mechanism of weed control. These herbicides exert strong selection pressure because of their high activity on sensitive biotypes at the rates used and because of their soil residual activity. Several point mutations within the gene encoding ALS can result in a herbicide-resistant ALS. From investigations of numerous R weed biotypes, five conserved amino acids have been identified in ALS that, on substitution, can confer resistance to ALS inhibitors. Substitutions of at least 12 additional ALS amino acids can also confer herbicide resistance in plants and other organisms but, to date, have not been found in R weed populations. Mutations in ALS conferring herbicide resistance are at least partially dominant, and because the gene is nuclear inherited, it is transmitted by both seed and pollen. Furthermore, in many cases there is apparently a negligible fitness cost of the resistance gene in the absence of herbicide selection. Although resistance to ALS-inhibiting herbicides has been a bane to weed management, it has spurred many advances within and beyond the weed science discipline. As examples, resistance to ALS-inhibiting herbicides has been exploited in the development of herbicide-resistant crops, studies of weed population dynamics, and in developing protocols for targeted gene modification. Resistance to ALS-inhibiting herbicides has greatly affected weed science by influencing how we view the sustainability of our weed management practices, what we consider when developing and marketing new herbicides, and how we train new weed scientists.

Understanding protein–inorganic surface interactions is central to the rational design of new tools in biomaterial sciences, nanobiotechnology and nanomedicine. Although a significant amount of experimental research on protein adsorption onto solid substrates has been reported, many aspects of the recognition and interaction mechanisms of biomolecules and inorganic surfaces are still unclear. Theoretical modeling and simulations provide complementary approaches for experimental studies, and they have been applied for exploring protein–surface binding mechanisms, the determinants of binding specificity towards different surfaces, as well as the thermodynamics and kinetics of adsorption. Although the general computational approaches employed to study the dynamics of proteins and materials are similar, the models and force-fields (FFs) used for describing the physical properties and interactions of material surfaces and biological molecules differ. In particular, FF and water models designed for use in biomolecular simulations are often not directly transferable to surface simulations and vice versa. The adsorption events span a wide range of time- and length-scales that vary from nanoseconds to days, and from nanometers to micrometers, respectively, rendering the use of multi-scale approaches unavoidable. Further, changes in the atomic structure of material surfaces that can lead to surface reconstruction, and in the structure of proteins that can result in complete denaturation of the adsorbed molecules, can create many intermediate structural and energetic states that complicate sampling. In this review, we address the challenges posed to theoretical and computational methods in achieving accurate descriptions of the physical, chemical and mechanical properties of protein-surface systems. In this context, we discuss the applicability of different modeling and simulation techniques ranging from quantum mechanics through all-atom molecular mechanics to coarse-grained approaches. We examine uses of different sampling methods, as well as free energy calculations. Furthermore, we review computational studies of protein–surface interactions and discuss the successes and limitations of current approaches.

Although there are many glutamate receptors in the retina, 2-amino-4-phosphonobutyrate (L-AP4) is an agonist that acts selectively at metabotropic glutamate receptors (mGluR6) of ON bipolar cells. We explored the properties of agonists that activate this receptor. The effects of various glutamate analogs on the b-wave of the electroretinogram (ERG) were used as a measure of their activity. Conformational comparisons among agonists suggest that ligands in an extended conformation preferentially bind to the ON bipolar synaptic receptor. But this property is insufficient to explain the selectivity of mGluR6 because some inactive glutamate analogs could also match this extended conformation. Comparative molecular field analysis (CoMFA) was used to compare the electrostatic and steric potentials of agonists with their action at the ON bipolar synapse. Steric potentials beneath a plane defined by the three putative binding sites plays a key role in determining agonist activity. The CoMFA model was used to predict the activity of glutamate analogs and correlations between predicted and measured activity support the model.

The active sites and substrate bindings of Rhizobium trifolii molonyl-CoA synthetase (MCS) catalyzing the malonyl-CoA formation from malonate and CoA have been determined based on NMR spectroscopy, site-directed mutagenesis, and comparative modeling methods. The MCS-bound conformation of malonyl-CoA was determined from two-dimensional–transferred nuclear Overhauser effect spectroscopy data. MCS protein folds into two structural domains and consists of 16 α-helices, 24 β-strands, and several long loops. The core active site was determined as a wide cleft close to the end of the small C-terminal domain. The catalytic substrate malonate is placed between ATP and His206 in the MCS enzyme, supporting His206 in its catalytic role as it generates reaction intermediate, malonyl-AMP. These findings are strongly supported by previous biochemical data, as well as by the site-directed mutagenesis data reported here. This structure reveals the biochemical role as well as the substrate specificity that conservative residues of adenylate-forming enzymes have.

The ecdysone receptor (ECR), a nuclear transcription factor controlling insect development, is a novel target for insecticides such as dibenzoylhydrazines with low environmental and toxicological impacts. To understand the high selectivity of such synthetic molecules toward ECR, two homology models of the Chironomus tentans ECR ligand-binding domain (LDB) have been constructed by taking as templates the known LBD crystal structures of the retinoic acid and vitamin D receptors. Docking of 20-hydroxyecdysone (20E) and dibenzoylhydrazines to the receptor suggests a novel superposition of the natural and synthetic molecules; the N-tert-butyl substituent of the dibenzoylhydrazines extends significantly beyond the 20E volume. Our ECR–LBD protein models rationalize how 20E and dibenzoylhydrazines interact with the ligand-binding pocket. The homology model complexes provide new insights that can be exploited in the rational design of new environmentally safe insecticides.

In a general approach to the understanding of protein adaptation to high temperature, molecular models of the closely related mesophilic Streptomyces sp. S38 Xyl1 and thermophilic Thermomonospora fusca TfxA family 11 xylanases were built and compared with the three-dimensional (3D) structures of homologous enzymes. Some of the structural features identified as potential contributors to the higher thermostability of TfxA were introduced in Xyl1 by site-directed mutagenesis in an attempt to improve its thermostability and thermophilicity. A new Y11–Y16 aromatic interaction, similar to that present in TfxA and created in Xyl1 by the T11Y mutation, improved both the thermophilicity and thermostability. Indeed, the optimum activity temperature (70 vs. 60 °C) and the apparent Tm were increased by about 9 °C, and the mutant was sixfold more stable at 57 °C. The combined mutations A82R/F168H/N169D/Δ170 potentially creating a R82–D169 salt bridge homologous to that present in TfxA improved the thermostability but not the thermophilicity. Mutations R82/D170 and S33P seemed to be slightly destabilizing and devoid of influence on the optimal activity temperature of Xyl1. Structural analysis revealed that residues Y11 and Y16 were located on β-strands B1 and B2, respectively. This interaction should increase the stability of the N-terminal part of Xyl1. Moreover, Y11 and Y16 seem to form an aromatic continuum with five other residues forming putative subsites involved in the binding of xylan (+3, +2, +1, −1, −2). Y11 and Y16 might represent two additional binding subsites (−3, −4) and the T11Y mutation could thus improve substrate binding to the enzyme at higher temperature and thus the thermophilicity of Xyl1.

We have determined the three-dimensional structure of the potassium channel inhibitor HsTX1, using nuclear magnetic resonance and molecular modeling. This protein belongs to the scorpion short toxin family, which essentially contains potassium channel blockers of 29 to 39 amino acids and three disulfide bridges. It is highly active on voltage-gated Kv1.3 potassium channels. Furthermore, it has the particularity to possess a fourth disulfide bridge. We show that HsTX1 has a fold similar to that of the three-disulfide–bridged toxins and conserves the hydrophobic core found in the scorpion short toxins. Thus, the fourth bridge has no influence on the global conformation of HsTX1. Most residues spatially analogous to those interacting with voltage-gated potassium channels in the three-disulfide–bridged toxins are conserved in HsTX1. Thus, we propose that Tyr21, Lys23, Met25, and Asn26 are involved in the biological activity of HsTX1. As an additional positively charged residue is always spatially close to the aromatic residue in toxins blocking the voltage-gated potassium channels, and as previous mutagenesis experiments have shown the critical role played by the C-terminus in HsTX1, we suggest that Arg33 is also important for the activity of the four disulfide-bridged toxin. Docking calculations confirm that, if Lys23 and Met25 interact with the GYGDMH motif of Kv1.3, Arg33 can contact Asp386 and, thus, play the role of the additional positively charged residue of the toxin functional site. This original configuration of the binding site of HsTX1 for Kv1.3, if confirmed experimentally, offers new structural possibilities for the construction of a molecule blocking the voltage-gated potassium channels.

Comparative sequence analysis complements experimental methods for the determination of RNA three-dimensional structure. This approach is based on the concept that different sequences within the same gene family form similar higher-order structures. The large number of rRNA sequences with sufficient variation, along with improved covariation algorithms, are providing us with the opportunity to identify new base triples in 16S rRNA. The three-dimensional conformations for one of our strongest candidates involving U121 (C124:G237) and/or U121 (U125:A236) (Escherichia coli sequence and numbering) are analyzed here with different molecular modeling tools. Molecular modeling shows that U121 interacts with C124 in the U121 (C124:G237) base triple. This arrangement maintains isomorphic structures for the three most frequent sequence motifs (approximately 93% of known bacterial and archaeal sequences), is consistent with chemical reactivity of U121 in E. coli ribosomes, and is geometrically favorable. Further, the restricted set of observed canonical (GU, AU, GC) base-pair types at positions 124:237 and 125:236 is consistent with the fact that the canonical base-pair sets (for both base pairs) that are not observed in nature prevent the formation of the 121(124:237) base triple. The analysis described here serves as a general scheme for the prediction of specific secondary and tertiary structure base pairing where there is a network of correlated base changes.

The human C3a anaphylatoxin receptor (C3aR) is a G protein-coupled receptor (GPCR) composed of seven transmembrane α-helices connected by hydrophilic loops. Previous studies of chimeric C3aR/C5aR and loop deletions in C3aR demonstrated that the large extracellular loop2 plays an important role in noneffector ligand binding; however, the effector binding site for C3a has not been identified. In this study, selected charged residues in the transmembrane regions of C3aR were replaced by Ala using site-directed mutagenesis, and mutant receptors were stably expressed in the RBL-2H3 cell line. Ligand binding studies demonstrated that R161A (helix IV), R340A (helix V), and D417A (helix VII) showed no binding activity, although full expression of these receptors was established by flow cytometric analysis. C3a induced very weak intracellular calcium flux in cells expressing these three mutant receptors. H81A (helix II) and K96A (helix III) showed decreased ligand binding activity. The calcium flux induced by C3a in H81A and K96A cells was also consistently reduced. These findings suggest that the charged transmembrane residues Arg161, Arg340, and Asp417 in C3aR are essential for ligand effector binding and/or signal coupling, and that residues His81 and Lys96 may contribute less directly to the overall free energy of ligand binding. These transmembrane residues in C3aR identify specific molecular contacts for ligand interactions that account for C3a-induced receptor activation.

Neurotrophins are a family of proteins with pleiotropic effects mediated by two distinct receptor types, namely the Trk family, and the common neurotrophin receptor p75NTR. Binding of four mammalian neurotrophins, nerve growth factor (NGF), brain-derived neurotrophic factor (BDNF), neurotrophin-3 (NT-3), and neurotrophin-4/5 (NT-4/5), to p75NTR is studied by molecular modeling based on X-ray structures of the neurotrophins and the extracellular domain of p55TNFR, a homologue of p75NTR. The model of neurotrophin/receptor interactions suggests that the receptor binding domains of neurotrophins (loops I and IV) are geometrically and electrostatically complementary to a putative binding site of p75NTR, formed by the second and part of the third cysteine-rich domains. Geometric match of neurotrophin/receptor binding domains in the complexes, as characterized by shape complementarity statistic Sc, is comparable to known protein/protein complexes. All charged residues within the loops I and IV of the neurotrophins, previously determined as being critical for p75NTR binding, directly participate in receptor binding in the framework of the model. Principal residues of the binding site of p75NTR include Asp47, Lys56, Asp75, Asp76, Asp88, and Glu89. The additional involvement of Arg80 and Glu53 is specific for NGF and BDNF, respectively, and Glu73 participates in binding with NT-3 and NT-4/5. Neurotrophins are likely to induce similar, but not identical, conformational changes within the p75NTR binding site.

B*2701 differs from all other HLA-B27 subtypes of known peptide specificity in that, among its natural peptide ligands, arginine is not the only allowed residue at peptide position 2. Indeed, B*2701 is unique in binding many peptides with Gln2 in vivo. However, the mutation (Asp74Tyr) responsible for altered selectivity is far away from the B pocket of the peptide binding site to which Gln/Arg2 binds. Here, we present a model that explains this effect. It is proposed that a new rotameric state of the conserved Lys70 is responsible for the unique B*2701 binding motif. This side chain should be either kept away from pocket B through its interaction with Asp74 in most HLA-B27 subtypes, or switched to this pocket if residue 74 is Tyr as in B*2701. Involvement of Lys70 in pocket B would thus allow binding of peptides with Gln2. Binding of Arg2-containing peptides to B*2701 is also possible because Lys70 could adopt another conformation, H-bonded to Asn97, which preserves the same binding mode of Arg2 as in B*2705. This model was experimentally validated by mutating Lys70 into Ala in B*2701. Edman sequencing of the B*2701(K70A) peptide pool showed only Arg2, characteristic of HLA-B27-bound peptides, and no evidence for Gln2. This supports the computational model and demonstrates that allowance of B*2701 for peptides with Gln2 is due to the long-range effect of the polymorphic residue 74 of HLA-B27, by inducing a conformational switch of the conserved Lys70.

Thermotoga maritima (Tm) expresses a 7 kDa monomeric protein whose 18 N-terminal amino acids show 81% identity to N-terminal sequences of cold shock proteins (Csps) from Bacillus caldolyticus and Bacillus stearothermophilus. There were only trace amounts of the protein in Thermotoga cells grown at 80 °C. Therefore, to perform physicochemical experiments, the gene was cloned in Escherichia coli. A DNA probe was produced by PCR from genomic Tm DNA with degenerated primers developed from the known N-terminus of TmCsp and the known C-terminus of CspB from Bacillus subtilis. Southern blot analysis of genomic Tm DNA allowed to produce a partial gene library, which was used as a template for PCRs with gene- and vector-specific primers to identify the complete DNA sequence. As reported for other csp genes, the 5′ untranslated region of the mRNA was anomalously long; it contained the putative Shine–Dalgarno sequence. The coding part of the gene contained 198 bp, i.e., 66 amino acids. The sequence showed 61% identity to CspB from B. caldolyticus and high similarity to all other known Csps. Computer-based homology modeling allowed the conclusion that TmCsp represents a β-barrel similar to CspB from B. subtilis and CspA from E. coli.

As indicated by spectroscopic analysis, analytical gel permeation chromatography, and mass spectrometry, overexpression of the recombinant protein yielded authentic TmCsp with a molecular weight of 7,474 Da. This was in agreement with the results of analytical ultracentrifugation confirming the monomeric state of the protein. The temperature-induced equilibrium transition at 87 °C exceeds the maximum growth temperature of Tm and represents the maximal Tm-value reported for Csps so far.

Site-specific photo crosslinking has been used to investigate the RNA neighborhood of 16S rRNA positions U788/U789 in Escherichia coli 30S subunits. For these studies, site-specific psoralen (SSP) which contains a sulfhydryl group on a 17 Å side chain was first added to nucleotides U788/U789 using a complementary guide DNA by annealing and phototransfer. Modified RNA was purified from the DNA and unmodified RNA. For some experiments, the SSP, which normally crosslinks at an 8 Å distance, was derivitized with azidophenacylbromide (APAB) resulting in the photoreactive azido moiety at a maximum of 25 Å from the 4′ position on psoralen (SSP25APA). 16S rRNA containing SSP, SSP25APA or control 16S rRNA were reconstituted and 30S particles were isolated. The reconstituted subunits containing SSP or SSP25APA had normal protein composition, were active in tRNA binding and had the usual pattern of chemical reactivity except for increased kethoxal reactivity at G791 and modest changes in four other regions. Irradiation of the derivatized 30S subunits in activation buffer produced several intramolecular RNA crosslinks that were visualized and separated by gel electrophoresis and characterized by primer extension. Four major crosslink sites made by the SSP reagent were identified at positions U561/U562, U920/U921, C866 and U723; a fifth major crosslink at G693 was identified when the SSP25APA reagent was used. A number of additional crosslinks of lower frequency were seen, particularly with the APA reagent. These data indicate a central location close to the decoding region and central pseudoknot for nucleotides U788/U789 in the activated 30S subunit.

The Pb2+ cleavage of a specific phosphodiester bond in yeast tRNAPhe is the classical model of metal-assisted RNA catalysis. In vitro selection experiments have identified a tRNAPhe variant, the leadzyme, that is very active in cleavage by Pb2+. We present here a three-dimensional modeling protocol that was used to propose a structure for this ribozyme, and is based on the computation of the intersection of conformational space of sequence variants and the use of chemical modification data. Sequence and secondary structure data were used in a first round of computer modeling that allowed identification of conformations compatible with all known leadzyme variants. Common conformations were then tested experimentally by evaluating the activity of analogues containing modified nucleotides in the catalytic core. These experiments led to a new structural hypothesis that was tested in a second round of computer modeling. The resulting proposal for the active conformation of the leadzyme is consistent with all known structural data. The final model suggests an in-line SN2 attack mechanism and predicts two Pb2+ binding sites. The protocol presented here is generally applicable in modeling RNAs whenever the catalytic or binding activity of structural analogues is known.

Two recently published but independently derived structures, namely the X-ray crystallographic structure of ribosomal protein S7 and the “binding pocket” for this protein in a three-dimensional model of the 16S rRNA, have been correlated with one another. The known rRNA–protein interactions for S7 include a minimum binding site, a number of footprint sites, and two RNA–protein crosslink sites on the 16S rRNA, all of which form a compact group in the published 16S rRNA model (despite the fact that these interactions were not used as primary modeling contraints in building that model). The amino acids in protein S7 that are involved in the two crosslinks to 16S rRNA have also been determined in previous studies, and here we have used these sites to orient the crystallographic structure of S7 relative to its rRNA binding pocket. Some minor alterations were made to the rRNA model to improve the fit. In the resulting structure, the principal positively charged surface of the protein is in contact with the 16S rRNA, and all of the RNA–protein interaction data are satisfied. The quality of the fit gives added confidence as to the validity of the 16S rRNA model. Protein S7 is furthermore known to be crosslinked both to P site-bound tRNA and to mRNA at positions upstream of the P site codon; the matched S7-16S rRNA structure makes a prediction as to the location of this crosslink site within the protein molecule.