Transition State Conformations Research Articles

Carbon Nanomaterials and DNA: from Molecular Recognition to Applications.

DNA is polymorphic. Increasing evidence has indicated that many biologically important processes are related to DNA's conformational transition and assembly states. In particular, noncanonical DNA structures, such as the right-handed A-form, the left-handed Z-form, the triplex, the G-quadruplex, the i-motif, and so forth, have been specific targets for the diagnosis and therapy of human diseases. Meanwhile, they have been widely used in the construction of smart DNA nanomaterials and nanoarchitectures. As rising stars in materials science, the family of carbon nanomaterials (CNMs), including two-dimensional graphene, one-dimensional carbon nanotubes (CNTs), and zero-dimensional graphene or carbon quantum dots (GQDs or CQDs), interact with DNA and are able to regulate the conformational transitions of DNA. The interaction of DNA with CNMs not only opens new opportunities for specific molecular recognition, but it also expands the promising applications of CNMs from materials science to biotechnology and biomedicine. In this Account, we focus on our contributions to the field of interactions between CNMs and DNA in which we have explored their promising applications in nanodevices, sensing, materials synthesis, and biomedicine. For one-dimensional CNTs, two-dimensional graphene, and zero-dimensional GQDs and CQDs, the basic principles, binding modes, and applications of the interactions between CNMs and DNA are reviewed. We aim to give prominence to the important status of CNMs in the field of molecular recognition for DNA. First, we summarized our discovery of the interactions between single-walled carbon nanotubes (SWNTs) with duplex, triplex, and human telomeric i-motif DNA and their interesting applications. For example, SWNTs are the first chemical agents that can selectively stabilize human telomeric i-motif DNA and induce its formation under physiological conditions. On the basis of this principle, two types of nanodevices were designed. One was used for highly sensitive detection of ppm levels of SWNTs in cells, and the other monitored i-motif DNA formation. Further studies indicated that SWNTs could inhibit telomerase activity in living cells and cause telomere dysfunction, providing new insight into the biological effects of SWNTs. Then, some applications that are based on the interactions between graphene and DNA are also summarized. Combined with other nanomaterials, such as metal and upconversion nanoparticles, several hybrid nanomaterials were successfully constructed, and a series of DNA logic gates were successfully developed. Afterwards, the newcomer of the carbon nanomaterials family, carbon quantum dots (CQDs), were found to be capable of modulating right-handed B-form DNA to left-handed Z-form DNA. These were further used to design FRET logic gates that were based on the CQD-derived DNA conformational transition. Taking into account the remaining challenges and promising aspects, CNM-based DNA nanotechnology and its biomedical applications will attract more attention and produce new breakthroughs in the near future.

SsDNA-Functionalized Nanoceria: A Redox-Active Aptaswitch for Biomolecular Recognition.

Quantification of biomolecular binding events is a critical step for the development of biorecognition assays for diagnostics and therapeutic applications. This paper reports the design of redox-active switches based on aptamer conjugated nanoceria for detection and quantification of biomolecular recognition. It is shown that the conformational transition state of the aptamer on nanoceria, combined with the redox properties of these particles can be used to create surface based structure switchable aptasensing platforms. Changes in the redox properties at the nanoceria surface upon binding of the ssDNA and its target analyte enables rapid and highly sensitive measurement of biomolecular interactions. This concept is demonstrated as a general applicable method to the colorimetric detection of DNA binding events. An example of a nanoceria aptaswitch for the colorimetric sensing of Ochratoxin A (OTA) and applicability to other targets is provided. The system can sensitively and selectivity detect as low as 0.15 × 10(-9) m OTA. This novel assay is simple in design and does not involve oligonucleotide labeling or elaborate nanoparticle modification steps. The proposed mechanism discovered here opens up a new way of designing optical sensing methods based on aptamer recognition. This approach can be broadly applicable to many bimolecular recognition processes and related applications.

Low-Mode Conformational Search Method with Semiempirical Quantum Mechanical Calculations: Application to Enantioselective Organocatalysis.

A conformational search program for finding low-energy conformations of large noncovalent complexes has been developed. A quantitatively reliable semiempirical quantum mechanical PM6-DH+ method, which is able to accurately describe noncovalent interactions at a low computational cost, was employed in contrast to conventional conformational search programs in which molecular mechanical methods are usually adopted. Our approach is based on the low-mode method whereby an initial structure is perturbed along one of its low-mode eigenvectors to generate new conformations. This method was applied to determine the most stable conformation of transition state for enantioselective alkylation by the Maruoka and cinchona alkaloid catalysts and Hantzsch ester hydrogenation of imines by chiral phosphoric acid. Besides successfully reproducing the previously reported most stable DFT conformations, the conformational search with the semiempirical quantum mechanical calculations newly discovered a more stable conformation at a low computational cost.

Conformational Interconversions of Amino Acid Derivatives.

Exhaustive conformational interconversions including transition structure analyses of N-acetyl-l-glycine-N-methylamide as well as its alanine, serine, and cysteine analogues have been investigated at the MP2/6-31G** level, yielding a total of 142 transition states. Improved estimates of relative energies were obtained by separately extrapolating the Hartree-Fock and MP2 energies to the basis set limit and adding the difference between CCSD(T) and MP2 results with the cc-pVDZ basis set to the extrapolated MP2 results. The performance of eight empirical force fields (AMBER94, AMBER14SB, MM2, MM3, MMFFs, CHARMM22_CMAP, OPLS_2005, and AMOEBAPRO13) in reproducing ab initio energies of transition states was tested. Our results indicate that commonly used class I force fields employing a fixed partial charge model for the electrostatic interaction provide mean errors in the ∼10 kJ/mol range for energies of conformational transition states for amino acid conformers. Modern reparametrized versions, such as CHARMM22_CMAP, and polarizable force fields, such as AMOEBAPRO13, have slightly lower mean errors, but maximal errors are still in the 35 kJ/mol range. There are differences between the force fields in their ability for reproducing conformational transitions classified according to backbone/side-chain or regions in the Ramachandran angles, but the data set is likely too small to draw any general conclusions. Errors in conformational interconversion barriers by ∼10 kJ/mol suggest that the commonly used force field may bias certain types of transitions by several orders of magnitude in rate and thus lead to incorrect dynamics in simulations. It is therefore suggested that information for conformational transition states should be included in parametrizations of new force fields.

Understanding metastable phase transformation during crystallization of RDX, HMX and CL-20: experimental and DFT studies.

Multiphase growth during crystallization severely affects deliverable output of explosive materials. Appearance and incomplete transformation of metastable phases are a major source of polymorphic impurities. This article presents a methodical and molecular level understanding of the metastable phase transformation mechanism during crystallization of cyclic nitramine explosives, viz. RDX, HMX and CL-20. Instantaneous reverse precipitation yielded metastable γ-HMX and β-CL-20 which undergo solution mediated transformation to the respective thermodynamic forms, β-HMX and ε-CL-20, following 'Ostwald's rule of stages'. However, no metastable phase, anticipated as β-RDX, was evidenced during precipitation of RDX, which rather directly yielded the thermodynamically stable α-phase. The γ→β-HMX and β→ε-CL-20 transformations took 20 and 60 minutes respectively, whereas formation of α-RDX was instantaneous. Density functional calculations were employed to identify the possible transition state conformations and to obtain activation barriers for transformations at wB97XD/6-311++G(d,p)(IEFPCM)//B3LYP/6-311G(d,p) level of theory. The computed activation barriers and lattice energies responsible for transformation of RDX, HMX and CL-20 metastable phases to thermodynamic ones conspicuously supported the experimentally observed order of phase stability. This precise result facilitated an understanding of the occurrence of a relatively more sensitive and less dense β-CL-20 phase in TNT based melt-cast explosive compositions, a persistent and critical problem unanswered in the literature. The crystalline material recovered from such compositions revealed a mixture of β- and ε-CL-20. However, similar compositions of RDX and HMX never showed any metastable phase. The relatively long stability with the highest activation barrier is believed to restrict complete β→ε-CL-20 transformation during processing. Therefore a method is suggested to overcome this issue.

A Single Glycosidase Harnesses Different Pyranoside Ring Transition State Conformations for Hydrolysis of Mannosides and Glucosides

Hydrolysis of β-d-mannosides by β-mannosidases typically proceeds via a B2,5 transition state conformation for the pyranoside ring, while that of β-d-glucosides by β-glucosidases proceeds through a distinct 4H3 transition state conformation. However, rice Os7BGlu26 β-glycosidase hydrolyzes 4-nitrophenyl β-d-glucoside and β-d-mannoside with similar efficiencies. The origin of this dual substrate specificity was investigated by kinetic, structural, and computational approaches. The glycosidase inhibitors glucoimidazole and mannoimidazole inhibited Os7BGlu26 with Ki values of 2.7 nM and 10.4 μM, respectively. In X-ray crystal structures of complexes with Os7BGlu26, glucoimidazole bound to the active site in a 4E conformation, while mannoimidazole bound in a B2,5 conformation, suggesting different transition state conformations. Moreover, calculation of quantum mechanics/molecular mechanics (QM/MM) free energy landscapes showed that 4-nitrophenyl β-d-glucoside adopts a 1S3/4E conformation in the Michaelis com...

Contributions of unique active site residues of eukaryotic UDP-galactopyranose mutases to substrate recognition and active site dynamics.

UDP-galactopyranose mutase (UGM) catalyzes the interconversion between UDP-galactopyranose and UDP-galactofuranose. Absent in humans, galactofuranose is found in bacterial and fungal cell walls and is a cell surface virulence factor in protozoan parasites. For these reasons, UGMs are targets for drug discovery. Here, we report a mutagenesis and structural study of the UGMs from Aspergillus fumigatus and Trypanosoma cruzi focused on active site residues that are conserved in eukaryotic UGMs but are absent or different in bacterial UGMs. Kinetic analysis of the variants F66A, Y104A, Q107A, N207A, and Y317A (A. fumigatus numbering) show decreases in kcat/KM values of 200–1000-fold for the mutase reaction. In contrast, none of the mutations significantly affect the kinetics of enzyme activation by NADPH. These results indicate that the targeted residues are important for promoting the transition state conformation for UDP-galactofuranose formation. Crystal structures of the A. fumigatus mutant enzymes were determined in the presence and absence of UDP to understand the structural consequences of the mutations. The structures suggest important roles for Asn207 in stabilizing the closed active site, and Tyr317 in positioning of the uridine ring. Phe66 and the corresponding residue in Mycobacterium tuberculosis UGM (His68) play a role as the backstop, stabilizing the galactopyranose group for nucleophilic attack. Together, these results provide insight into the essentiality of the targeted residues for realizing maximal catalytic activity and a proposal for how conformational changes that close the active site are temporally related and coupled together.

Capturing Transition Paths and Transition States for Conformational Rearrangements in the Ribosome

To reveal the molecular determinants of biological function, one seeks to characterize the interactions that are formed in conformational and chemical transition states. In other words, what interactions govern the molecule’s energy landscape? To accomplish this, it is necessary to determine which degrees of freedom can unambiguously identify each transition state. Here, we perform simulations of large-scale aminoacyl-transfer RNA (aa-tRNA) rearrangements during accommodation on the ribosome and project the dynamics along experimentally accessible atomic distances. From this analysis, we obtain evidence for which coordinates capture the correct number of barrier-crossing events and accurately indicate when the aa-tRNA is on a transition path. Although a commonly used coordinate in single-molecule experiments performs poorly, this study implicates alternative coordinates along which rearrangements are accurately described as diffusive movements across a one-dimensional free-energy profile. From this, we provide the theoretical foundation required for single-molecule techniques to uncover the energy landscape governing aa-tRNA selection by the ribosome.

Simulations and experiments in protein folding.

The interplay between simulations and experiments of protein folding has largely contributed to the elucidation of many important aspects of the phenomenon. In this chapter, I briefly describe the experiments which provide information on the kinetics of the protein folding process, and help to characterize the folding transition state. Then, I show how to probe the kinetics of protein folding using molecular dynamics simulations, how to compare the simulations with the experiments and how to help and rationalize the latter, ultimately offering a molecular picture of the process. After the production of suitable molecular dynamics simulation data in the form of trajectories, the procedure involves sequentially the identification of the stable states of the protein, the identification of the transition pathways connecting the stable states, the identification of the transition state conformations, comparison with experimental results, and finally, the identification of the molecular determinants or reaction coordinates of the folding process, that is, the features that clearly help distinguishing the transition state from the stable states.

α-Fluorophosphonates reveal how a phosphomutase conserves transition state conformation over hexose recognition in its two-step reaction.

β-Phosphoglucomutase (βPGM) catalyzes isomerization of β-D-glucose 1-phosphate (βG1P) into D-glucose 6-phosphate (G6P) via sequential phosphoryl transfer steps using a β-D-glucose 1,6-bisphosphate (βG16BP) intermediate. Synthetic fluoromethylenephosphonate and methylenephosphonate analogs of βG1P deliver novel step 1 transition state analog (TSA) complexes for βPGM, incorporating trifluoromagnesate and tetrafluoroaluminate surrogates of the phosphoryl group. Within an invariant protein conformation, the β-D-glucopyranose ring in the βG1P TSA complexes (step 1) is flipped over and shifted relative to the G6P TSA complexes (step 2). Its equatorial hydroxyl groups are hydrogen-bonded directly to the enzyme rather than indirectly via water molecules as in step 2. The (C)O-P bond orientation for binding the phosphate in the inert phosphate site differs by ∼ 30° between steps 1 and 2. By contrast, the orientations for the axial O-Mg-O alignment for the TSA of the phosphoryl group in the catalytic site differ by only ∼ 5°, and the atoms representing the five phosphorus-bonded oxygens in the two transition states (TSs) are virtually superimposable. The conformation of βG16BP in step 1 does not fit into the same invariant active site for step 2 by simple positional interchange of the phosphates: the TS alignment is achieved by conformational change of the hexose rather than the protein.

To be or not to be butterfly: New mechanistic insights in the Aza-Michael asymmetric addition of lithium (R)-N-benzyl-N-(α-methylbenzyl)amide.

The asymmetric Aza-Michael addition of homochiral lithium benzylamides to α,β-unsaturated esters represents an extended protocol to obtain enantioenriched β-amino esters. An exhaustive mechanistic revision of the originally proposed mechanism is reported, developing a quantum mechanics/molecular mechanics protocol for the asymmetric Aza-Michael reaction of homochiral lithium benzylamides. Explicit and implicit solvent schemes were considered, together with a proper account of long-range dispersion forces, evaluated through a density functional theory benchmark of different functionals. Theoretical results showed that the diastereoselectivity is mainly controlled by the N-α-methylbenzyl moiety placing, deriving a Si/Re 99:1 diastereoselective ratio, in good agreement with reported experimental results. The main transition state geometries are two transition state conformers in a "V-stacked" orientation of the amide's phenyl rings, differing in the tetrahydrofuran molecule arrangement coordinated to the metal center. Extensive conformational sampling and quantum-level refinement give reasonable good speed/accuracy results, allowing this protocol to be extended to other similar Aza-Michael reaction systems.

Molecular mechanism of a hotdog-fold acyl-CoA thioesterase.

Thioesterases are enzymes that hydrolyze thioester bonds between a carbonyl group and a sulfur atom. They catalyze key steps in fatty acid biosynthesis and metabolism, as well as polyketide biosynthesis. The reaction molecular mechanism of most hotdog-fold acyl-CoA thioesterases remains unknown, but several hypotheses have been put forward in structural and biochemical investigations. The reaction of a human thioesterase (hTHEM2), representing a thioesterase family with a hotdog fold where a coenzyme A moiety is cleaved, was simulated by quantum mechanics/molecular mechanics metadynamics techniques to elucidate atomic and electronic details of its mechanism, its transition-state conformation, and the free energy landscape of the process. A single-displacement acid-base-like mechanism, in which a nucleophilic water molecule is activated by an aspartate residue acting as a base, was found, confirming previous experimental proposals. The results provide unambiguous evidence of the formation of a tetrahedral-like transition state. They also explain the roles of other conserved active-site residues during the reaction, especially that of a nearby histidine/serine pair that protonates the thioester sulfur atom, the participation of which could not be elucidated from mutation analyses alone.

ChemInform Abstract: Supramolecular Catalysis. Part 2. Artificial Enzyme Mimics

AbstractReview: 852 refs.

Stretching of single poly-ubiquitin molecules revisited: Dynamic disorder in the non-exponential unfolding kinetics

A theoretical framework based on a generalized Langevin equation (GLE) with fractional Gaussian noise (fGn) and a power-law memory kernel is presented to describe the non-exponential kinetics of the unfolding of a single poly-ubiquitin molecule under a constant force [T.-L. Kuo, S. Garcia-Manyes, J. Li, I. Barel, H. Lu, B. J. Berne, M. Urbakh, J. Klafter, and J. M. Fernández, Proc. Natl. Acad. Sci. U.S.A. 107, 11336 (2010)]. Such a GLE-fGn strategy is made on the basis that the pulling coordinate variable x undergoes subdiffusion, usually resulting from conformational fluctuations, over a one-dimensional force-modified free-energy surface U(x, F). By using the Kramers' rate theory, we have obtained analytical formulae for the time-dependent rate coefficient k(t, F), the survival probability S(t, F) as well as the waiting time distribution function f(t, F) as functions of time t and force F. We find that our results can fit the experimental data of f(t, F) perfectly in the whole time range with a power-law exponent γ = 1/2, the characteristic of typical anomalous subdiffusion. In addition, the fitting of the survival probabilities for different forces facilitates us to reach rather reasonable estimations for intrinsic properties of the system, such as the free-energy barrier and the distance between the native conformation and the transition state conformation along the reaction coordinate, which are in good agreements with molecular dynamics simulations in the literatures. Although static disorder has been implicated in the original work of Kuo et al., our work suggests a sound and plausible alternative interpretation for the non-exponential kinetics in the stretching of poly-ubiquitin molecules, associated with dynamic disorder.

Specific Catalysis of Asparaginyl Deamidation by Carboxylic Acids: Kinetic, Thermodynamic, and Quantitative Structure–Property Relationship Analyses

Asparaginyl (Asn) deamidation could lead to altered potency, safety, and/or pharmacokinetics of therapeutic protein drugs. In this study, we investigated the effects of several different carboxylic acids on Asn deamidation rates using an IgG1 monoclonal antibody (mAb1*) and a model hexapeptide (peptide1) with the sequence YGKNGG. Thermodynamic analyses of the kinetics data revealed that higher deamidation rates are associated with predominantly more negative ΔS and, to a lesser extent, more positive ΔH. The observed differences in deamidation rates were attributed to the unique ability of each type of carboxylic acid to stabilize the energetically unfavorable transition-state conformations required for imide formation. Quantitative structure property relationship (QSPR) analysis using kinetic data demonstrated that molecular descriptors encoding for the geometric spatial distribution of atomic properties on various carboxylic acids are effective determinants for the deamidation reaction. Specifically, the number of O-O and O-H atom pairs on carboxyl and hydroxyl groups with interatomic distances of 4-5 Å on a carboxylic acid buffer appears to determine the rate of deamidation. Collectively, the results from structural and thermodynamic analyses indicate that carboxylic acids presumably form multiple hydrogen bonds and charge-charge interactions with the relevant deamidation site and provide alignment between the reactive atoms on the side chain and backbone. We propose that carboxylic acids catalyze deamidation by stabilizing a specific, energetically unfavorable transition-state conformation of l-asparaginyl intermediate II that readily facilitates bond formation between the γ-carbonyl carbon and the deprotonated backbone nitrogen for cyclic imide formation.

Supramolecular catalysis. Part 2: artificial enzyme mimics

The design of artificial catalysts able to compete with the catalytic proficiency of enzymes is an intense subject of research. Non-covalent interactions are thought to be involved in several properties of enzymatic catalysis, notably (i) the confinement of the substrates and the active site within a catalytic pocket, (ii) the creation of a hydrophobic pocket in water, (iii) self-replication properties and (iv) allosteric properties. The origins of the enhanced rates and high catalytic selectivities associated with these properties are still a matter of debate. Stabilisation of the transition state and favourable conformations of the active site and the product(s) are probably part of the answer. We present here artificial catalysts and biomacromolecule hybrid catalysts which constitute good models towards the development of truly competitive artificial enzymes.

Possible Intervention of Prototropic and Silatropic Ene-type Mechanisms in Asymmetric Catalytic Mukaiyama Aldol Reactions

Prototropic and silatropic ene mechanisms were analyzed in the BINOL-TiIV, bisoxazoline (box)-CuII, and BINAP-PdII complex-catalyzed asymmetric Mukaiyama aldol reaction of methyl glyoxylate and acetone-derived 2-trimethylsiloxypropene. Our ab initio calculations indicate that the BINOL-TiIV-catalyzed Mukaiyama aldol reaction proceeds through the concerted prototropic or silatropic ene pathways and that box-CuII- and BINAP-PdII-catalyzed reactions proceed through the stepwise prototropic ene pathways. The transition-state conformations are thus syn-clinal with hydrogen and silane as the migratory groups aligned to the formyl oxygen atom in close syn-clinal proximity, respectively.

Combined Inhibitor Free‐Energy Landscape and Structural Analysis Reports on the Mannosidase Conformational Coordinate

Mannosidases catalyze the hydrolysis of a diverse range of polysaccharides and glycoconjugates, and the various sequence-based mannosidase families have evolved ingenious strategies to overcome the stereoelectronic challenges of mannoside chemistry. Using a combination of computational chemistry, inhibitor design and synthesis, and X-ray crystallography of inhibitor/enzyme complexes, it is demonstrated that mannoimidazole-type inhibitors are energetically poised to report faithfully on mannosidase transition-state conformation, and provide direct evidence for the conformational itinerary used by diverse mannosidases, including β-mannanases from families GH26 and GH113. Isofagomine-type inhibitors are poor mimics of transition-state conformation, owing to the high energy barriers that must be crossed to attain mechanistically relevant conformations, however, these sugar-shaped heterocycles allow the acquisition of ternary complexes that span the active site, thus providing valuable insight into active-site residues involved in substrate recognition.

Combined Inhibitor Free-Energy Landscape and Structural Analysis Reports on the Mannosidase Conformational Coordinate

Mannosidases catalyze the hydrolysis of a diverse range of polysaccharides and glycoconjugates, and the various sequence-based mannosidase families have evolved ingenious strategies to overcome the stereoelectronic challenges of mannoside chemistry. Using a combination of computational chemistry, inhibitor design and synthesis, and X-ray crystallography of inhibitor/enzyme complexes, it is demonstrated that mannoimidazole-type inhibitors are energetically poised to report faithfully on mannosidase transition-state conformation, and provide direct evidence for the conformational itinerary used by diverse mannosidases, including β-mannanases from families GH26 and GH113. Isofagomine-type inhibitors are poor mimics of transition-state conformation, owing to the high energy barriers that must be crossed to attain mechanistically relevant conformations, however, these sugar-shaped heterocycles allow the acquisition of ternary complexes that span the active site, thus providing valuable insight into active-site residues involved in substrate recognition.

Analyses of Conformational States of the Transporter Associated with Antigen Processing (TAP) Protein in a Native Cellular Membrane Environment

The transporter associated with antigen processing (TAP) plays a critical role in the MHC class I antigen presentation pathway. TAP translocates cellular peptides across the endoplasmic reticulum membrane in an ATP hydrolysis-dependent manner. We used FRET spectroscopy in permeabilized cells to delineate different conformational states of TAP in a native subcellular membrane environment. For these studies, we tagged the TAP1 and TAP2 subunits with enhanced cyan fluorescent protein and enhanced yellow fluorescent protein, respectively, C-terminally to their nucleotide binding domains (NBDs), and measured FRET efficiencies under different conditions. Our data indicate that both ATP and ADP enhance the FRET efficiencies but that neither induces a maximally closed NBD conformation. Additionally, peptide binding induces a large and significant increase in NBD proximity with a concentration dependence that is reflective of individual peptide affinities for TAP, revealing the underlying mechanism of peptide-stimulated ATPase activity of TAP. Maximal NBD closure is induced by the combination of peptide and non-hydrolysable ATP analogs. Thus, TAP1-TAP2 NBD dimers are not fully stabilized by nucleotides alone, and substrate binding plays a key role in inducing the transition state conformations of the NBD. Taken together, these findings show that at least three steps are involved in the transport of peptides across the endoplasmic reticulum membrane for antigen presentation, corresponding to three dynamically and structurally distinct conformational states of TAP. Our studies elucidate structural changes in the TAP NBD in response to nucleotides and substrate, providing new insights into the mechanism of ATP-binding cassette transporter function.