Theoretical and computational strategies for rational molecularly imprinted polymer design

The further evolution of molecularly imprinted polymer science and technology necessitates the development of robust predictive tools capable of handling the complexity of molecular imprinting systems. A combination of the rapid growth in computer power over the past decade and significant software developments have opened new possibilities for simulating aspects of the complex molecular imprinting process. We present here a survey of the current status of the use of in silico-based approaches to aspects of molecular imprinting. Finally, we highlight areas where ongoing and future efforts should yield information critical to our understanding of the underlying mechanisms sufficient to permit the rational design of molecularly imprinted polymers.     

Study of the Chelating Capacity of Nucleobase Analogs with Biological Interest: Ab initio Molecular Orbital Calculations and Single-Crystal X-ray Structural Study of 6-Amino-5-formyl-1,3-dimethyluracil-benzoylhydrazone

Ab initio molecular orbital calculations at the RHF/6–31G* level using the GAUSSIAN94 program package have allowed us to simulate the molecular structures for different conformations of 6-amino-5-formyl-1,3-dimethyluracil-benzoylhydrazone. The contribution of the atomic orbitals of the potential donor atoms to the higher occupied molecular orbitals allows us to propose theoretical arguments to justify the different chelating behavior found for this compound in several metal complexes. Further, the molecular structure of 6-amino-5-formyl-1,3-dimethyluracil-benzoylhydrazone has been determined by single-crystal X-ray diffraction methods. The compound crystallizes in the monoclinic system (space group P2₁/n) with cell dimensions: a = 12.111(5), b = 5.743(5), c = 22.636(5) Å, = 98.60(5)°. The structure was solved from 1719 reflections with I>2(I). The final R [I>2(I)] was 0.0506 for 217 parameters. The azomethine double bond substituents are in the E conformation and the N51 atom is in the cis position with respect to the N6 atom, due to the formation of an intramolecular hydrogen bond N6-H···N51. The geometrical data are in good agreement with those calculated by means of ab initio methods.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089400060630     

DFT and MP2 studies on the C2C2α bond cleavage in thiamin catalysis

Thiamin diphosphate (ThDP), the biologically active derivative of vitamin B₁, is an important cofactor of several enzymes that catalyze the oxidative and non-oxidative conversion of α-keto acids. The final step of non-oxidative decarboxylation of pyruvate by pyruvate decarboxylase – the liberation of acetaldehyde – requires deprotonation of the α-hydroxyl group and cleavage of the C2–C2α bond of the transitory 2-(1-hydroxyethyl)-ThDP intermediate. It has been proposed that the cofactor 4′-amino/imino function is essentially involved in the deprotonation of the α-hydroxyl group. Proton transfer and C2–C2α cleavage may occur in a stepwise manner, or, alternatively in a concerted mechanism. Here, density functional theory (DFT) calculations as well as second order Møller–Plesset perturbation theory (MP2) studies were performed on a simple model for the enzyme using the program package Gaussian 03. Calculations favor a stepwise mechanism with initial formation of the C2α alkoxide, followed by C2–C2α bond cleavage.     

Calculations of the C<Subscript>2</Subscript> fragmentation energies of higher fullerenes C<Subscript>80</Subscript> and C<Subscript>82</Subscript>

The C₂ fragmentation energies of the most stable isolated-pentagon-rule (IPR) isomers of the C₈₀ and C₈₂ fullerenes were evaluated with second-order Møller-Plesset (MP2) theory, density-functional theory (DFT) and the semiempirical self-consistent charge density-functional tight-binding (SCC-DFTB) method. Zero-point energy, ionization energy and empirical C₂ corrections were included in the calculation of fragmentation energies for comparison with experimental C₂ fragmentation energies of the fullerene cations. In the case of the most probable Stone-Wales pathway of C₂ fragmentation of C₈₀, the calculated \(D_{0} {\left( {{\text{C}}_{{{\text{80}}}} ^{ + } } \right)}\) agree well with experimental data, whereas in the case of C₈₂ fragmentation, the calculated \(D_{0} {\left( {{\text{C}}_{{{\text{82}}}} ^{ + } } \right)}\) exceed by up to 1.2 eV the experimental ones, which suggests that other IPR isomers may be present in sufficient amounts in experimental samples. Computer-intensive MP2 calculations and DFT calculations with larger basis sets do not yield much improved C₂ fragmentation energies, compared to those reported earlier with B3LYP/3-21G. On the other hand, semiempirical approaches such as SCC-DFTB, which are orders of magnitude less intensive, yield satisfactory fragmentation energies for higher fullerenes and may become a method of choice for routine calculations of fullerenes and carbon nanotubes.

Open image in new window

Figure C₂ fragmentation energies of C₈₀ and C₈₂ fullerenes have been calculated with B3LYP/6-31G* model chemistry, with semiempirical self-consistent-charge density-functional tight-binding (SCC-DFTB) method and with the more rigorous MP2 method. The influence of basis set extension and level of theory on the resulting fragmentation energies is discussed

     

An update view on the substrate recognition mechanism of phosphodiesterases: a computational study of PDE10 and PDE4 bound with cyclic nucleotides

Early studies strongly implied that the specificity of cyclic nucleotide phosphodiesterases (PDEs) toward its endogenous substrates can be uniquely determined by the amido orientation of the invariant glutamine locating in the binding pocket of the enzyme. However, recently solved crystal structures of PDE4 (cAMP specific) and PDE10 (dual specific) in the presence of endogenous substrates have revealed that their invariant glutamine orientations are very similar despite exhibiting different substrate specificities proven physiologically. To understand this subtle specificity issue in the PDE family, here several experimentally inaccessible PDE-substrate complex models have been studied computationally, and the results are juxtaposed and compared in detail. Modeling results show that PDE10 in fact favors cAMP energetically but still can bind to cGMP owing to the robust hydrogen-bond network in the vicinity of the invariant glutamine side chain. PDE4 fails to accommodate cGMP is correlated to the weakening of this same hydrogen-bond network but not owing to any steric strain in the binding pocket. An Asn residue in the binding pocket of PDE4 has enhanced the specificity of the binding to cAMP sideway as observed in our computer simulation. Further to the previously studied syn- versus anti-conformational specificity of cAMP in PDE10, the unexpected substrate-binding mode in PDE10 versus PDE4 as reported here strongly suggested that there are remaining uncertainties in the substrate orientation and recognition mechanism in the PDE families. The molecular details of the binding pocket observed in this study provide hints for more optimal PDE4 and PDE10 inhibitor design.     

A computational investigation of sulfur-containing heterocyclic components from the anal sac secretions of Mustela species

A computational investigation of the sulfur-containing heterocyclic components (substituted thietanes and 1,2-dithiolanes) of Mustela anal sac secretions has been carried out. A cluster analysis of the chemical compositions of Mustela anal sac volatiles reveals little similarity with established phylogenetic relationships between members of the genus. Ab initio calculations [MP2/6–311++G(2df,2p)//B3LYP/6–311++G**] show the lowest-energy C₅H₁₀S isomeric thietane to be 2,2-dimethylthietane, which is also the most abundant of the Mustela thietanes. Similarly, 3,3-dimethyl-1,2-dithiolane is the lowest-energy C₅H₁₀S₂ compound. 2-n-Propylthietane is the highest-energy C₆H₁₂S compound, but the most abundant Mustela C₆H₁₂S compound produced, whereas cis-2-ethyl-4-methylthietane, the lowest-energy C₆H₁₂S thietane, has never been observed in Mustela anal sac secretions. A molecular docking analysis of the Mustela sulfur-containing heterocycles into both porcine and bovine odorant binding proteins reveals the interactions of the docked ligands with the proteins to be largely hydrophobic, and have binding energies generally lower than typical odorant molecules such as linalool or eugenol. Figure Mustela anal sac volatile components, 2,2-dimethylthietane and cis-3,4-dimethyl–1,2-dithiolane.     

Exploring the structural basis of the selective inhibition of monoamine oxidase A by dicarbonitrile aminoheterocycles: Role of Asn181 and Ile335 validated by spectroscopic and computational studies

Since cyanide potentiates the inhibitory activity of several monoamine oxidase (MAO) inhibitors, a series of carbonitrile-containing aminoheterocycles was examined to explore the role of nitriles in determining the inhibitory activity against MAO. Dicarbonitrile aminofurans were found to be potent, selective inhibitors against MAO A. The origin of the MAO A selectivity was identified by combining spectroscopic and computational methods. Spectroscopic changes induced in MAO A by mono- and dicarbonitrile inhibitors were different, providing experimental evidence for distinct binding modes to the enzyme. Similar differences were also found between the binding of dicarbonitrile compounds to MAO A and to MAO B. Stabilization of the flavin anionic semiquinone by monocarbonitrile compounds, but destabilization by dicarbonitriles, provided further support to the distinct binding modes of these compounds and their interaction with the flavin ring. Molecular modeling studies supported the role played by the nitrile and amino groups in anchoring the inhibitor to the binding cavity. In particular, the results highlight the role of Asn181 and Ile335 in assisting the interaction of the nitrile-containing aminofuran ring. The network of interactions afforded by the specific attachment of these functional groups provides useful guidelines for the design of selective, reversible MAO A inhibitors.     

Engineering Strontium Binding Affinity in an EF-hand Motif: A Quantum Chemical and Molecular Dynamics Study

Abstract

Proteins with the ability to specifically bind strontium would potentially be of great use in the field of nuclear waste management. Unfortunately, no such peptides or proteins are known—indeed, it is uncertain whether they exist under natural conditions due to low environmental concentrations of strontium. To investigate the possibility of devising such molecules, one of us (CV), in a previous experimental study [J. Biol. Inorg. Chem. 8, 33440 (2003)], proposed starting from an EF-hand motif of the protein calmodulin and mutating some residues to change the motif's specificity for calcium into one for strontium. In this paper, which represents a theoretical complement to the experimental work, we analyzed small-molecule crystallographic structures and performed quantum chemical calculations to identify possible mutations. We then constructed seven mutant sequences of the EF-hand motif and analyzed their dynamical and binding behaviors using molecular dynamics simulations and free-energy calculations (using the MM/PBSA method). As a result of these analyzes we were able to isolate some characteristics that could lead to mutant peptides with enhanced strontium affinity.     

The structure-function relationship of the Aspergillus fumigatuscyp51A L98H conversion by site-directed mutagenesis: the mechanism of L98H azole resistance

Since 1998, the rapid emergence of multi-azole-resistance (MAR) was observed in Aspergillus fumigatus in the Netherlands. Two dominant mutations were found in the cyp51A gene, a 34 bp tandem repeat (TR) in the promoter region combined with a leucine to histidine substitution at codon 98 (L98H). In this study, we show that molecular dynamics simulations combined with site-directed mutagenesis of amino acid substitutions in the cyp51A gene, correlate to the structure–function relationship of the L98H substitution conferring to MAR in A. fumigatus. Because of a L98H directed change in the flexibility of the loops, that comprise a gate-like structure in the protein, the capacity of the two ligand entry channels is modified by narrowing the diameter and thereby binding of azoles is obstructed. Moreover, the L98H induced relocation of tyrosine 121 and tyrosine 107 seems to be related to the MAR phenotype, without affecting the biological activity of the CYP51A protein. Site-directed mutagenesis showed that both the 34 bp TR and the L98H mutation are required to obtain the MAR phenotype. Furthermore, the amino acid leucine in codon 98 in A. fumigatus is highly conserved and important for maintaining the structure of the CYP51A protein that is essential for azole docking.     

Mutability of DNA polymerase I: implications for the creation of mutant DNA polymerases

DNA polymerases of the Family A catalyze the addition of deoxynucleotides to a primer with high efficiency, processivity, and selectivity-properties that are critical to their function both in nature and in the laboratory. These polymerases tolerate many amino acid substitutions, even in regions that are evolutionarily conserved. This tolerance can be exploited to create DNA polymerases with novel properties and altered substrate specificities, using rational design and molecular evolution. These efforts have focused mainly on the Family A DNA polymerises -Taq, E. coli Pol I, and T7 - because they are widely utilized in biotechnology today. The redesign of polymerases often requires knowledge of the function of specific residues in the protein, including those located in six evolutionarily conserved regions. The most well characterized of these are motifs A and B, which regulate the fidelity of replication and the incorporation of nucleotide analogs such as dideoxynucleotides. Regions that remain to be more thoroughly characterized are motif C, which is critical for catalysis, and motifs 1, 2 and 6, all of which bind to DNA primer or template. Several recently identified mutants with abilities to incorporate nucleotides with bulky adducts have mutations that are not located within conserved regions and warrant further study. Analysis of these mutants will help advance our understanding of how DNA polymerases select bases with high fidelity.     

H2 production by mixed cultures

Production of H₂ from glucose by an anoxygenic phototrophic bacterium (Rhodobacter sphaeroides), a cyanobacterium (Synechococcus cedrorum) and a heterotrophic bacterium (Pseudomonas fluorescens) was tested individually and in mixed cultures of various combinations in light. H₂ production was maximal with a mixed culture of R. sphaeroides and P. fluorescens, which could be further enhanced by immobilization of the bacteria in alginate gel. Inhibition of H₂ photoproduction was observed in a mixture of S. cedrorum and P. fluorescens and a co-culture of all the three organisms.Ch. Sasikala and Ch. V. Ramana are and G. S. Prasad was with the Microbial Biotechnology Laboratory, Department of Botany, Osmania University, Hyderabad-500 007, India. G. S. Prasad is now with the Microbial Type Culture Collection Centre (MTCC), IMTECH, Chandigar, India.     

Prebiotic Synthesis of Simple Sugars by an Interstellar Formose Reaction

The prebiotic possibilities for the synthesis of interstellar carbohydrates through a protic variant of the formose reaction under gas phase conditions were studied. Ab initio calculations were used to evaluate potential mechanisms. Based on considerations of barrier heights and temperature variations in the Interstellar Medium the plausibility of extended sugar synthesis will be discussed.     

Chiral-at-metal tetrahydrosalen complexes of resolved titanium(IV) sec-butoxides: Ligand wrapping and multiple asymmetric catalytic induction

The reaction of the racemic and resolved tetrahydrosalen derivative LH₂ (LH₂ = N,N’-bis(3,5-dichloro-2-hydroxybenzyl)-trans-1,2-diaminocyclohexane) with the resolved titanium(IV) sec-butoxides Ti(O^R-2Bu)₄ or Ti(O^S-2Bu)₄ yielded a series of four compounds, LTi(O²Bu)₂ (1-4), which have been characterized by IR, elemental analysis, ¹H and ¹³C NMR and X-ray crystallography. X-ray crystallography revealed the co-crystallization of two pseudo-C₂-symmetric products from racemic LH₂, whereas a perfect chiral induction of the ligand to the metal occurred when resolved (R,R)-LH₂ was used, resulting in a Δ fac-fac wrapping mode of the tetradentate ligand about the titanium center. Ab initio electronic structure calculations (DFT) are in agreement that the lowest energy isomer is that which is experimentally observed. Catalysis screenings show that Ti(O^S-2Bu)₄, in conjunction with (R,R)-LH₂, forms a matched pair that catalyzes the addition of dimethyl zinc to benzaldehyde with higher enantioselectivity than that observed for resolved (R,R)-LH₂ with Ti(O^R-2Bu)₄ or achiral Ti(OⁱPr)₄. Increasing the temperature of the system results in slightly increased enantiomeric excess.     

Screening of the active site from water by the incoming ligand triggers catalysis and inhibition in serine proteases

The pKa of the catalytic His57 N(epsilon)H in the tetrahedral complex (TC) of chymotrypsin with trifluoromethyl ketone inhibitors is 4-5 units higher relative to the free enzyme (FE). Such stable TC's, formed with transition state (TS) analog inhibitors, are topologically similar to the catalytic TS. Thus, analysis of this pKa shift may shed light on the role of water solvation in the general base catalysis by histidine. We applied our QM/SCRF(VS) approach to study this shift. The method enables explicit quantum mechanical DFT calculations of large molecular clusters that simulate chemical reactions at the active site (AS) of water solvated enzymes. We derived an analytical expression for the pKa dependence on the degree of water exposure of the ionizable group, and on the total charge in the enzyme AS, Q(A) and Q(B), when the target ionizable functional group (His57 in this study) is in the acidic (A) and basic (B) forms, respectively. Q2(B) > Q2(A) both in the FE and in the TC of chymotrypsin. Therefore, water solvation decreases the relative stability of the protonated histidine in both. Ligand binding reduces the degree of water solvation of the imidazole ring, and consequently elevates the histidine pKa. Thus, the binding of the ligand plays a triggering role that switches on the cascade of catalytic reactions in serine proteases.     

TCNE-aniline charge transfer complex: ab initio and TDDFT investigations in gas phase

The geometric and electronic structure of tetracyanoethylene (TCNE)-aniline (donor-acceptor type) complex has been investigated in gas phase using ab initio and time dependent density functional theory calculations. Both the above calculations predict a composed structure for the complex, in which the interacting site is a C≡N and C=C bond center in the TCNE and, –NH₂ and π-electrons of aniline. The N atom of aniline is oriented toward the TCNE molecule. The charge transfer transition energy, estimated by calculating the ground-to-excited state transition electric dipole moments of the complex, agree well with the reported experimental value in chloroform medium. TCNE-aniline at ground state. TCNE-aniline at excited state     

Structural and computational basis for potent inhibition of glutamate carboxypeptidase II by carbamate-based inhibitors

A series of carbamate-based inhibitors of glutamate carboxypeptidase II (GCPII) were designed and synthesized using ZJ-43, N-[[[(1S)-1-carboxy-3-methylbutyl]amino]carbonyl]-l-glutamic acid, as a molecular template in order to better understand the impact of replacing one of the two nitrogen atoms in the urea-based GCPII inhibitor with an oxygen atom. Compound 7 containing a C-terminal 2-oxypentanedioic acid was more potent than compound 5 containing a C-terminal glutamic acid (2-aminopentanedioic acid) despite GCPII’s preference for peptides containing an N-terminal glutamate as substrates. Subsequent crystallographic analysis revealed that ZJ-43 and its two carbamate analogs 5 and 7 with the same (S,S)-stereochemical configuration adopt a nearly identical binding mode while (R,S)-carbamate analog 8 containing a d-leucine forms a less extensive hydrogen bonding network. QM and QM/MM calculations have identified no specific interactions in the GCPII active site that would distinguish ZJ-43 from compounds 5 and 7 and attributed the higher potency of ZJ-43 and compound 7 to the free energy changes associated with the transfer of the ligand from bulk solvent to the protein active site as a result of the lower ligand strain energy and solvation/desolvation energy. Our findings underscore a broader range of factors that need to be taken into account in predicting ligand-protein binding affinity. These insights should be of particular importance in future efforts to design and develop GCPII inhibitors for optimal inhibitory potency.     

Ribonucleotide reductase inhibition by p-alkoxyphenols studied by molecular docking and molecular dynamics simulations

Ribonucleotide reductase (RNR) is necessary for production of the precursor deoxyribonucleotides for DNA synthesis. Class Ia RNR functions via a stable free radical in one of the two components protein R2. The enzyme mechanism involves long range (proton coupled) electron transfer between protein R1 and the tyrosyl radical in protein R2. Earlier experimental studies showed that p-alkoxyphenols inhibit RNR. Here, molecular docking and molecular dynamics simulations involving protein R2 suggest an inhibition mechanism for p-alkoxyphenols . A low energy binding pocket is identified in protein R2. The preferred configuration provides a structural basis explaining their specific binding to the Escherichia coli and mouse R2 proteins. Trp48 (E. coli numbering), on the electron transfer pathway, is involved in the interactions with the inhibitors. The relative order of the binding energies calculated for the phenol derivatives to protein R2 is correlated with earlier experimental data on inhibition efficiency, in turn related to increasing size of the hydrophobic alkyl substituents. Using the configuration identified by molecular docking as a starting point for molecular dynamics simulations, we find that the p-allyloxyphenol interrupts the catalytic electron transfer pathway of the R2 protein by forming hydrogen bonds with Trp48 and Asp237, thus explaining the inhibitory activity of p-alkoxyphenols.     

Determinative factors in inhibition of aquaporin by different pharmaceuticals: Atomic scale overview by molecular dynamics simulation

The inhibition of water permeation through aquaporins by ligands of pharmaceutical compounds is considered as a method to control the cell lifetime. The inhibition of aquaporin 1 (AQP1) by bacopaside-I and torsemide, was explored and its atomistic nature was elucidated by molecular docking and molecular dynamics (MD) simulation collectively along with Poisson-Boltzmann surface area (PBSA) method. Docking results revealed that torsemide has a lower level of docking energy in comparison with bacopaside-I at the cytoplasmic side. Furthermore, the effect of steric constraints on water permeation was accentuated. Bacopaside-I inhibits the channel properly due to the strong interaction with the channel and larger spatial volume, whereas torsemide blocks the cytoplasmic side of the channel imperfectly. The most probable active sites of AQP1 for the formation of hydrogen bonds between the inhibitor and the channel were identified by numerical analysis of the bonds. Eventually, free energy assessments indicate that binding of both inhibitors is favorable in complex with AQP1, and van der Waals interaction has an important contribution in stabilizing the complexes.     

Experimental and computational analyses of the energetic basis for dual recognition of immunity proteins by colicin endonucleases

Colicin endonucleases (DNases) are bound and inactivated by immunity (Im) proteins. Im proteins are broadly cross-reactive yet specific inhibitors binding cognate and non-cognate DNases with K_d values that vary between 10^− 4 and 10^− 14 M, characteristics that are explained by a ‘dual-recognition’ mechanism. In this work, we addressed for the first time the energetics of Im protein recognition by colicin DNases through a combination of E9 DNase alanine scanning and double-mutant cycles (DMCs) coupled with kinetic and calorimetric analyses of cognate Im9 and non-cognate Im2 binding, as well as computational analysis of alanine scanning and DMC data. We show that differential ΔΔGs observed for four E9 DNase residues cumulatively distinguish cognate Im9 association from non-cognate Im2 association. E9 DNase Phe86 is the primary specificity hotspot residue in the centre of the interface, which is coordinated by conserved and variable hotspot residues of the cognate Im protein. Experimental DMC analysis reveals that only modest coupling energies to Im9 residues are observed, in agreement with calculated DMCs using the program ROSETTA and consistent with the largely hydrophobic nature of E9 DNase-Im9 specificity contacts. Computed values for the 12 E9 DNase alanine mutants showed reasonable agreement with experimental ΔΔG data, particularly for interactions not mediated by interfacial water molecules. ΔΔG predictions for residues that contact buried water molecules calculated using solvated rotamer models met with mixed success; however, we were able to predict with a high degree of accuracy the location and energetic contribution of one such contact. Our study highlights how colicin DNases are able to utilise both conserved and variable amino acids to distinguish cognate from non-cognate Im proteins, with the energetic contributions of the conserved residues modulated by neighbouring specificity sites.