Interactions of aryloxyphenoxypropionic acids with sensitive and resistant acetyl-coenzyme a carboxylase by homology modeling and molecular dynamic simulations

Acetyl-coenzyme A carboxylase (ACCase) has been identified as one of the most important targets of herbicides. In the present study, we constructed homology models of the carboxyl-transferase (CT) domain of ACCase from sensitive and resistant foxtail and used these models as templates to study the molecular mechanism of herbicide resistance and stereochemistry-activity relationships of aryloxyphenoxypropionates (APPs). In the homology modeling structures, the dimer of the CT domain was formed by the side-to-side arrangement of the two monomers, in such a way that the N domain of one molecule is placed next to the C domain of the other. The dimeric association of sensitive foxtail CT was found to differ from that of resistant foxtail CT, and the spatial orientation of two key residues, Leu-695 and Ile-695, in these dimers also differed. The mutation of Ile to Leu may perturb the conformation of the dimeric interface, which may account for the molecular mechanism of herbicide resistance. Further docking analysis indicated that the binding model of high-active compounds is similar to that in the crystal structure of the enzyme-ligand complex. The different spatial orientations of ester groups of the isomers of APPs may explain the stereochemistry-activity relationship. Ser-698 formed a H-bonding interaction with all of the docked ligands, while Tyr-728 formed a pi-pi stacking interaction with some of the APPs. These findings may enhance our understanding of the molecular mechanism of herbicide resistance and stereochemistry-activity relationships, which may provide a new starting point for the identification of more potent inhibitors against both sensitive and resistant ACCase.     

分子动力学模拟研究点突变(Met108→Leu108)对树胶醛糖结合蛋白与配体作用的影响

为了研究点突变(Met108→Leu108)对树胶醛糖结合蛋白(ABP)与配体结合能力的影响,对ABP、ABP结合树胶醛糖复合物及ABP结合半乳糖复合物以及它们各自的突变体分别进行60 ns的分子动力学模拟.模拟结果表明,108号残基突变前后,电子等排体的两个氨基酸残基,使蛋白与配体间的范德华相互作用发生明显变化,同时导致蛋白的内部运动也发生变化,进而影响蛋白与配体的相互作用.进一步分析表明,突变前后的蛋白构象变化都趋向于两个结构域张开,而与配体的结合可减缓张开程度.     

Interactions between cycloguanil derivatives and wild type and resistance-associated mutant Plasmodium falciparum dihydrofolate reductases

Abstract  Comparative molecular field analysis (CoMFA) and quantum chemical calculations were performed on cycloguanil (Cyc) derivatives of the wild type and the quadruple mutant (Asn51Ile, Cys59Arg, Ser108Asn, Ile164Leu) of Plasmodium falciparum dihydrofolate reductase (PfDHFR). The represented CoMFA models of wild type ( and r ² = 0.985) and mutant type ( and r ² = 0.979) can describe the differences of the Cyc structural requirements for the two types of PfDHFR enzymes and can be useful to guide the design of new inhibitors. Moreover, the obtained particular interaction energies between the Cyc and the surrounding residues in the binding pocket indicated that Asn108 of mutant enzyme was the cause of Cyc resistance by producing steric clash with p-Cl of Cyc. Consequently, comparing the energy contributions with the potent flexible WR99210 inhibitor, it was found that the key mutant residue, Asn108, demonstrates attractive interaction with this inhibitor and some residues, Leu46, Ile112, Pro113, Phe116, and Leu119, seem to perform as second binding site with WR99210. Therefore, quantum chemical calculations can be useful for investigating residue interactions to clarify the cause of drug resistance. Graphical Abstract  CoMFA steric contour maps of Cyc derivatives against the quadruple mutant PfDHFR. Electrostatic van der Waals surface interaction of Cyc and the key resistance residue Asn108.      

Importance of individual side chains for the stability of a protein fold: computational alanine scanning of the insulin monomer

A new computational approach is proposed to probe the importance of residue side chains for the stability of a protein fold. Computational mutations to estimate protein stability (CMEPS) is based on the notion that the binding free energy corresponding to the complexation of a given side chain, considered as a "pseudo-ligand" of the wild type protein, reflects the importance of this side chain to the thermodynamic stability of the protein. The contribution of a particular side chain to the folding energy is estimated according to the molecular mechanics-generalized born surface area MM-GBSA approach, using a single molecular dynamics simulation trajectory of the wild type protein. CMEPS is a first principles method which does not contain any adjustable parameter that could be fitted to experimental data. The approach is first validated for Barnase and the B1 domain of protein L, for which a correlation coefficient R = 0.73, between experimental and CMEPS calculated DeltaDeltaG values, is found and then applied to the insulin monomer. In the present application, CMEPS replaces each amino acid by an alanine residue. Therefore, most mutations lead to cavities in the protein. From this the change in stability can be correlated with increased cavity volume. For insulin, this correlation is very similar compared with data previously analyzed for T4 lysozyme from an experiment for buried apolar side chains. There, the increased cavity volume has been related to the hydrophobic effect. However, since CMEPS uses the energetics in terms of electrostatic and van der Waals interactions (and not the hydrophobic effect which is difficult to relate to physical interactions), it is possible to study the effect of mutations of polar and solvent accessible side chains. According to CMEPS, residues Leu A16, Tyr A19, Leu B11, Leu B15, and Arg B22 are most important for the stability of the monomeric insulin fold. This is in agreement with experimental data. As a consequence, mutation of these residues may lead to misfolded and inactive insulin analogues.     

Insight into herbicide resistance of W574L mutant Arabidopsis thaliana acetohydroxyacid synthase: molecular dynamics simulations and binding free energy calculations

Acetohydroxyacid synthase (AHAS) is the target enzyme of several classes of herbicides, such as sulfonylureas and imidazolinones. Now many mutant AHASs with herbicide resistance have emerged along with extensive use of herbicides, therefore it is imperative to understand the detailed interaction mechanism and resistance mechanism so as to develop new potent inhibitors for wild-type or resistant AHAS. With the aid of available crystal structures of the Arabidopsis thaliana (At) AHAS-inhibitor complex, molecular dynamics (MD) simulations were used to investigate the interaction and resistance mechanism directly and dynamically at the atomic level. Nanosecond-level MD simulations were performed on six systems consisting of wild-type or W574L mutant AtAHAS in the complex with three sulfonylurea inhibitors, separately, and binding free energy was calculated for each system using the MM-GBSA method. Comprehensive analyses from structural and energetic aspects confirmed the importance of residue W574, and also indicated that W574L mutation might alert the structural charactersistic of the substrate access channel and decrease the binding affinity of inhibitors, which cooperatively weaken the effective channel-blocked effect and finally result in weaker inhibitory effect of inhibitor and corresponding herbicide resistance of W574L mutant. To our knowledge, it is the first report about MD simulations study on the AHAS-related system, which will pave the way to study the interactions between herbicides and wild-type or mutant AHAS dynamically, and decipher the resistance mechanism at the atomic level for better designing new potent anti-resistance herbicides.     

The relationship between Van der Waals volume of 4-substituent of phenol and sulfation in the rat.

The influence of 4-substituent of phenol on the sulfate conjugation was studied in the rat. 4-Substituents used were 4-ethyl, 4-tert-butyl and 4-phenyl. The in vivo and in vitro sulfations decreased in the order of phenol, 4-ethylphenol, 4-tert-butylphenol and 4-phenylphenol. The activity of sulfotransferase was correlated with the van der Waals volume of 4-substituent. The modification of thiol and arginyl residues of cytosol protein showed a similar degree of inactivation of sulfotransferase to all substituents used. Inactivation by the modification of histidyl residue of the protein increased with the increase of the van der Waals volume of 4-substituent. This result suggests that histidyl residue recognized van der Waals volume of 4-substituent.     

Molecular dynamics simulation and free energy calculation studies of the binding mechanism of allosteric inhibitors with p38α MAP kinase

p38 MAP kinase is a promising target for anti-inflammatory treatment. The classical kinase inhibitors imatinib and sorafenib as well as BI-1 and BIRB-796 were reported to bind in the DFG-out form of human p38α, known as type II or allosteric kinase inhibitors. Although DFG-out conformation has attracted great interest in the design of type II kinase inhibitors, the structural requirements for binding and mechanism of stabilization of DFG-out conformation remain unclear. As allosteric inhibition is important to the selectivity of kinase inhibitor, herein the binding modes of imatinib, sorafenib, BI-1 and BIRB-796 to p38α were investigated by molecular dynamics simulation. Binding free energies were calculated by molecular mechanics/Poisson-Boltzmann surface area method. The predicted binding affinities can give a good explanation of the activity difference of the studied inhibitors. Furthermore, binding free energies decomposition analysis and further structural analysis indicate that the dominating effect of van der Waals interaction drives the binding process, and key residues, such as Lys53, Gly71, Leu75, Ile84, Thr106, Met109, Leu167, Asp168, and Phe169, play important roles by forming hydrogen bond, salt bridge, and hydrophobic interactions with the DFG-out conformation of p38α. Finally, we also conducted a detailed analysis of BI-1, imatinib, and sorafenib binding to p38α in comparison with BIRB-796 exploited for gaining potency as well as selectivity of p38 inhibitors. These results are expected to be useful for future rational design of novel type II p38 inhibitors.     

An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis.

The Src-homology-3 (SH3) domain of the Caenorhabditis elegans protein Sem-5 binds proline-rich sequences. It is reported that the SH3 domains broadly accept amide N-substituted residues instead of only recognizing prolines on the basis of side chain shape or rigidity. We have studied the interactions between Sem-5 and its ligands using molecular dynamics (MD), free energy calculations, and sequence analysis. Relative binding free energies, estimated by a method called MM/PBSA, between different substitutions at sites -1, 0, and +2 of the peptide are consistent with the experimental data. A new method to calculate atomic partial charges, AM1-BCC method, is also used in the binding free energy calculations for different N-substitutions at site -1. The results are very similar to those obtained from widely used RESP charges in the AMBER force field. AM1-BCC charges can be calculated more rapidly for any organic molecule than can the RESP charges. Therefore, their use can enable a broader and more efficient application of the MM/PBSA method in drug design. Examination of each component of the free energy leads to the construction of van der Waals interaction energy profiles for each ligand as well as for wild-type and mutant Sem-5 proteins. The profiles and free energy calculations indicate that the van der Waals interactions between the ligands and the receptor determine whether an N- or a Calpha-substituted residue is favored at each site. A VC value (defined as a product of the conservation percentage of each residue and its van der Waals interaction energy with the ligand) is used to identify several residues on the receptor that are critical for specificity and binding affinity. This VC value may have a potential use in identifying crucial residues for any ligand-protein or protein-protein system. Mutations at two of those crucial residues, N190 and N206, are examined. One mutation, N190I, is predicted to reduce the selectivity of the N-substituted residue at site -1 of the ligand and is shown to bind similarly with N- and Calpha-substituted residues at that site.     

Van der Waals surface graphs and molecular shape

A van der Waals surface graph is the graph defined on a van der Waals surface by the intersections of the atomic van der Waals spheres. A van der Waals shape graph has a vertex for each atom with a visible face on the van der Waals surface, and edges between vertices representing atoms with adjacent faces on the van der Waals surface. These are discrete invariants of three‐dimensional molecular shape. Some basic properties of van der Waals surface graphs are studied, including their relationship with the Voronoi diagram of the atom centres, and a class of molecular embeddings is identified for which the dual of the van der Waals surface graph coincides with the van der Waals shape graph. This revised version was published online in July 2006 with corrections to the Cover Date.     

Multiple states of the Tyr318Leu mutant of dihydroorotate dehydrogenase revealed by single-molecule kinetics

Dihydroorotate dehydrogenase (DHOD) from Escherichia coli is a monomeric membrane-associated flavoprotein that catalyzes the oxidation of dihydroorotate to orotate. By using confocal fluorescence spectroscopy on the highly fluorescent Tyr318Leu DHOD mutant, we studied the catalytic turnover of single enzyme molecules through the characteristic on-off fluorescence signal, which corresponds to flavin mononucleotide (FMN) interconverting between the oxidized and reduced states during turnover. Our single-molecule data provide evidence of a distinct static heterogeneity in the enzymatic activity, with some molecules going through the on-off cycles 5-fold faster than others, however, there is no detectable dynamic disorder in DHOD turnover. When 0.1% reduced Triton X-100, a detergent that more closely simulates the natural membrane environment, is added, our data suggest the degree of static molecular heterogeneity is reduced. The observation of static heterogeneity suggests that the enzyme, which associates with the membrane in vivo, is present in distinct conformations that result in different catalytic efficiencies. The alternate conformations are most likely the result of the loss of van der Waals or other interactions between tyrosine 318 and FMN in the catalytic site with the mutation of Tyr318Leu, which disrupts the native structure of wild-type DHOD.     

小麦ACCase CT功能域基因在大肠杆菌中的表达及与除草剂的相互作用

以中国春小麦幼苗为材料,克隆构建了小麦质体乙酰辅酶A羧化酶(ACCase)的羧基转移酶(CT)重组质粒( RCP18-5),并实现了重组质粒在大肠杆菌中的可溶性高表达.对重组蛋白的性质研究表明,该蛋白具有较强的疏水性,稳定性不高.为改善这种状况,对CT功能域基因进行了截短和延长,同样于大肠杆菌中进行表达.结果表明,仅长...     

A molecular model of a point mutation (Val297Met) in the serine protease domain of protein C

A heterozygous GTG to ATG (Val297Met) mutation was detected in a patient with inherited protein C deficiency and deep vein thrombosis. Cosegregation of the mutation with protein C deficiency was observed through a family pedigree study. Molecular models of the serine protease domains of wild type and mutant protein C were constructed by standard comparative method. Val 297 was found to be located in the hydrophobic core of the protein. Although the substitution of Met for Val does not greatly alter the hydrophobicity of the protein, it introduces a bulkier side chain, which yields steric hindrance between this residue and adjacent residues, such as Met364, Tyr393, Ile321, Ile323, and Val378. It seems that the Met can not fit into the tight packing into which it is trapped, thereby probably inducing misfolding and/or greater instability of the protein. Such misfolding and/or instability thereby eventually disturbs the catalytic triad, in consistent with the observed type I deficiency state.     

Tunneling through weak interactions: comparison of through-space-, H-bond-, and through-bond-mediated tunneling

Results from ab initio electronic structure theory calculations on model systems allow for the detailed comparison of tunneling through covalently bonded contacts, hydrogen bonds, and van der Waals contacts. Considerable geometrical sensitivity as well as an exponential distance dependence of the tunneling is observed for tunneling through various nonbonded contacts. However, the fundamental result from the present study is that at most a modest difference is observed between tunneling mediated by H-bonds and tunneling mediated by van der Waals contacts at typical distances for each type of interaction. These results are considered in relation to the pathways model of Beratan and Onuchic, and implications for understanding long-range tunneling in biological systems are discussed.     

聚合物流体的范德华凝聚行为

研究了由聚合物的范德华作用导致的凝聚行为. 研究发现, 尽管聚合物同小分子的相行为的形成原因不同(聚合物体系的相行为是由动能、构象熵项和范德华作用能三项相互竞争的结果, 而小分子的相行为是由动能和范德华作用能相互竞争的结果), 但是它们表现出了极为相似的相行为.     

Extended electron distributions applied to the molecular mechanics of some intermolecular interactions. II. Organic complexes

Summary Extended electron distributions (XEDs) have been used to simulate the formation of complexes by intermolecular interaction via: (i) aromatic stacking; and (ii) hydrogen bonding. The results qualitatively reproduce experimental observations. In contrast, atom-centred partial charges fail to reproduce highly hydrogen-bonded systems, but make little difference in cases where interactions are driven largely by van der Waals forces. The dielectric constant used in the Coulombic term has been shown to be significant in defining the type and properties of these interactions when XEDs are employed. Some consideration has been given to solvation and entropy effects.     

Cloning, Expression and Purification of Wheat Acetyl-CoA Carboxylases CT Domain in E.coil

The entire gene of carboxyltransferase(CT) domain of acetyl-CoA carboxylase(ACCase) from Chinese Spring wheat(CSW) plastid was cloned firstly,and the 2.3 kb gene was inserted into PET28a+ vector and expressed in E.coil in a soluble state.The (His)6 fusion protein was identified by SDS-PAGE and Western blot.The recombinant protein was purified by affinity chromatography,and the calculated molecular mass(Mr) was 88000.The results of the sequence analysis indicate that the cloned gene(GeneBank accession No.EU124675) was a supplement and revision of the reported ACCase CT partial cDNA from Chinese Spring wheat plastid.The recombinant protein will be significant for us to investigate the recognizing mechanism between ACCase and herbicides,and further to screen new herbicides.     

An Extension of the Consistent Valence Force Field (CVFF) with the Aim to Simulate the Structures of Vanadium Phosphorus Oxides and the Adsorption of n‐Butane and of 1‐Butene on their Crystal Planes

A specific force field of Consistent Valence Force Field type was developed with the aim to simulate the structures of catalysts of vanadium phosphorus oxide type and the reversible adsorption of organic compounds on specific crystallographic planes of such catalysts by molecular modeling. The appropriate parameters were derived for the bonded (stretching, bending, and torsional deformations) and nonbonded (attractive and repulsive van der Waals and Coulomb forces) atomic interactions for V—O and P—O bonds in typical fragments of these catalysts with the vanadium atom in the oxidation state IV. The parameters for bonded interactions were computed from Hessian matrices, supplied by the program DMol for performing Density Functional Theory, by means of a program for non‐linear regression. The DMol program was applied to energy minimize structures of known vanadium phosphorus oxides, which were compared with X‐ray structures, and to obtain their Hessian matrices as a basis for the force constants needed. Some hypothetical structural models had to be added. The van der Waals parameters were estimated by means of correlations between van der Waals radii and the repulsive parameters and between polarizabilities and the dispersive parameters from the literature. The force field obtained was applied to simulate the crystal structure of vanadyl pyrophosphate and to compute the heat of adsorption of n‐butane and of 1‐butene on its (100) plane (computer codes of company Biosym/MSI/Accelrys). The experimental crystal structure and the adsorption energies were fairly well reproduced, except that the a lattice constant proves somewhat too large.     

Electrostatic complementarity between proteins and ligands. 3. Structural basis

Summary Electrostatic potential complementarity between ligands and their receptor sites is evaluated by the superposition of the electrostatic potential, generated by the receptor, onto the ligand potential over the ligand van der Waals surface. We would like to examine which structural factors generate this pattern of superposition. Example studies suggest that in many ligand-protein pairs, there exist principal formal charges on each molecule, largely responsible for the electrostatic potential complementarity observed. Electrostatic potential complementarity depends on the relative disposition of these principal charges and the ligand van der Waals surface. Simple mathematical models were constructed to predict the complementarity solely from structural considerations. The essential conditions for electrostatic potential complementarity were elucidated. These can be used in ligand design strategies to obtain an electrostatically optimal ligand.     

On the decisive role of finite chain extensibility and global interactions in networks

The stress-strain dependence of dry networks at unidirectional extension and compression is studied. The phenomenological van der Waals equation of state is compared with different molecular models in order to provide an interpretation of the van der Waals corrections. It is shown that the stress-strain behavior predicted by the phantom, Langevin, and constrained junction fluctuation models are altogether covered by the van der Waals approach. The relationship between the suppression of junction fluctuation parameter introduced by Dossin and Graessley and the van der Waals corrections has been worked out. The effect of junction functionality on the small strain modulus as well as on the second Mooney-Rivlin coefficient is also presented.