The role of lysine-41 in RNase A catalysis - A quantum chemical study on the active site ligand complex

Several mechanisms have been proposed to explain the action of enzymes at the atomic level. Among them, the recent proposals involving short hydrogen bonds as a step in catalysis by Gerlt and Gassman [1] and proton transfer through low barrier hydrogen bonds (LBHBs) [2, 3] have attracted attention. There are several limitations to experimentally testing such hypotheses, Recent developments in computational methods facilitate the study of active site-ligand complexes to high levels of accuracy, Our previous studies, which involved the docking of the dinucleotide substrate UpA to the active site of RNase A [4, 5], enabled us to obtain a realistic model of the ligand-bound active site of RNase A. From these studies, based on empirical potential functions, we were able to obtain the molecular dynamics averaged coordinates of RNase A, bound to the ligand UpA. A quantum mechanical study is required to investigate the catalytic process which involves the cleavage and formation of covalent bonds. In the present study, we have investigated the strengths of some of the hydrogen bonds between the active site residues of RNase A and UpA at the ab initio quantum chemical level using the molecular dynamics averaged coordinates as the starting point. The 49 atom system and other model systems were optimized at the 3-21G level and the energies of the optimized systems were obtained at the 6-31G* level. The results clearly indicate the strengthening of hydrogen bonds between neutral residues due to the presence of charged species at appropriate positions. Such a strengthening manifests itself in the form of short hydrogen bonds and a low barrier for proton transfer. In the present study, the proton transfer between the 2'-OH of ribose (from the substrate) and the imidazole group from the H12 of RNase A is influenced by K41, which plays a crucial role in strengthening the neutral hydrogen bond, reducing the barrier for proton transfer.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

Effect of Structural Modification on the Quantum-Size Effect in II-VI Semiconducting Nanocrystals

The electronic structure of group II-VI semiconductors in the stable wurtzite form is analyzed using state-of-the-art ab initio approaches to extract a simple and chemically transparent tight-binding model. This model can be used to understand the variation in the bandgap with size, for nanoclusters of these compounds. Results complement similar information already available for same systems in the zinc blende structure. A comparison with all available experimental data on quantum size effects in group II-VI semiconductor nanoclusters establishes a remarkable agreement between theory and experiment in both structure types, thereby verifying the predictive ability of our approach. The significant dependence of the quantum size effect on the structure type suggests that the experimental bandgap change at a given size compared to the bulk bandgap, may be used to indicate the structural form of the nanoclusters, particularly in the small size limit, where broadening of diffraction features often make it difficult to unambiguously determine the structure.

Studies involving electrostatic potential of valinomycin

The electrostatic potential of valinomycin in various conformations as obtained by the crystal structures (uncomplexed, complexed) and theoretical considerations have been evaluated and compared. The potential energy profiles along the æ axis of the bracelet-like structures show a systematic variation from the uncomplexed to the complexed structure. This type of conformational change and the potential variation are probably associated with different states of ion transport, like the capture and release of ions by the ionophore. Also, the asymmetry of the molecule due to D-HyIV on one side and L-Lac on the other side is reflected in the potential values along the Z-axis, the magnitude of which, is considerable in the uncomplexed structure. The evaluation of the potential at the ab-initio level on smaller fragments indicate that the order of liganding capacity of oxygen is amide ether ester. Also, the inductive effects due to alkyl substitution is negligible as evidenced by the potential studies on the substituted amides and esters.

Systematic variation in bond lengths of peptides

The lengths of the carbonyl as well as of the adjacent C-N and C-C bonds in peptides are shown to vary systematically with the central C-N bond length. Results of ab initio calculations on N-methylacetamide and its Li+, Na+ and Mg2+ complexes are also discussed.

Design, Synthesis, and DNA Binding Properties of Photoisomerizable Azobenzene−Distamycin Conjugates: An Experimental and Computational Study

Here, we present the synthesis, photochemical, and DNA binding properties of three photoisomerizable azobenzene−distamycin conjugates in which two distamycin units were linked via electron-rich alkoxy or electron-withdrawing carboxamido moieties with the azobenzene core. Like parent distamycin A, these molecules also demonstrated AT-specific DNA binding. Duplex DNA binding abilities of these conjugates were found to depend upon the nature and length of the spacer, the location of protonatable residues, and the isomeric state of the conjugate. The changes in the duplex DNA binding efficiency of the individual conjugates in the dark and with their respective photoirradiated forms were examined by circular dichroism, thermal denaturation of DNA, and Hoechst displacement assay with poly[d(A-T).d(T-A)] DNA in 150 mM NaCl buffer. Computational structural analyses of the uncomplexed ligands using ab initio HF and MP2 theory and molecular docking studies involving the conjugates with duplex d[(GC(AT)10CG)]2 DNA were performed to rationalize the nature of binding of these conjugates.

Growth Mechanism of Cadmium Sulfide Nanocrystals

We investigate the growth kinetics of CdS nanocrystals in the quantum confinement regime using time-resolved small-angle X-ray scattering. In contrast to earlier reports for similar systems, we establish that the growth kinetics in this case follows the Lifshitz-Slyozov-Wagner theory, for not only growth of the average diameter of the nanocrystals but also the time dependence of the size distribution and the temperature dependence of the rate constant. This is the first rigorous example of the coarsening process in the quantum confinement (< 5 nm)regime. Ab initio studies for the reaction pathways provide a microscopic understanding of this finding.

Asynchronous behaviour in some three-atom linear concerted radical-exchange reactions

A few simple three-atom thermoneutral radical exchange reactions (i.e. A + BC --> AB + C) are examined by ab initio SCF methods. Emphasis is laid on the detailed analysis of density matrices rather than on energetics. Results reveal that the sum of the bond orders of the breaking and forming bonds is not conserved to unity, due to development of free valence on the migrating atom 'B' in the transition state. Bond orders, free valence and spin densities on the atoms are calculated. The present analysis shows that the bond-cleavage process is always more advanced than the bond-formation process in the transition state. Further analysis shows a development of the negative spin density on the migrating atom 'B' in the transition state. The depletion of the alpha-spin density on the radical site "A" in the reactant during the reaction lags behind the growth of the alpha-spin density on the terminal atom "C" of the reactant bond, 'B-C' in the transition state. But all these processes are completed simultaneously at the end of the reaction. Hence, the reactions are asynchronous but kinetically concerted in most cases.

Geometry changes induced by negative hyperconjugative interactions involving carbonyl and thiocarbonyl groups

MNDO geometry optimizations have been carried out on a series of acyclic and cyclic unsymmetrically disubstituted carbonyl and thiocarbonyl compounds. The C=X unit shows a consistent and often sizeable tilt towards one of the substituents, following the order O > Snot, vert, similarN > C > B. Reference ab initio calculations and available experimental results support the MNDO results. The effect, which is particularly dramatic in small rings, is attributed primarily to favorable negative hyperconjugative interaction between the lone pair on X and a low lying adjacent σ* orbital. Such an interaction can lead to highly distorted structures, including perhaps to a planar molecule with an inverted sp2 carbon center.

Molecular structure-activity relationship study of some non-steroidal antiinflammatory agents using electrostatic potential mapping

A series of 6,11-dihydro-11-oxodibenz[b,e]oxepin-2-acetic acids (DOAA) which are known to be anti-inflammatory agents were studied. The geometries of some of the molecules obtained from X-ray crystallography were used in the calculations as such while the geometries of their derivatives were obtained by local, partial geometry optimization around the Sites of substitution employing the AMI method, keeping the remaining parts of the geometries the same as those in the parent molecules. Molecular electrostatic potential (MEP) mapping was performed for the molecules using optimized hybridization displacement charges (HDC) combined with Lowdin charges, as this charge distribution has been shown earlier to yield near ab initio quality results. A good correlation has been found between the MEP values near the oxygen atoms of the hydroxyl groups of the carboxy groups of the molecules and their anti-inflammatory activities. The result is broadly in agreement with the model proposed earlier by other authors regarding the structure-activity relationship for other similar molecules.

Vibrational analysis of boldine hydrochloride using QM/MM approach

A complete vibrational analysis was performed on the molecular structure of boldine hydrochloride using QM/MM method. The equilibrium geometry, harmonic vibrational frequencies and infrared intensities were calculated by QM/MM method with B3LYP/6-31G(d) and universal force field (UFF) combination using ONIOM code. We found the geometry obtained by the QM/MM method to be very accurate, and we can use this rapid method in place of time consuming ab initio methods for large molecules. A detailed interpretation of the infrared spectra of boldine hydrochloride is reported. The scaled theoretical wave numbers are in perfect agreement with the experimental values. The FT-IR spectra of boldine hydrochloride in the region 4000-500 cm(-1) were recorded in CsI (solid phase) and in chloroform with concentration 5 and 10 mg/ml.

Size-Dependent Tuning of Mn2+ d Emission in Mn2+-Doped CdS Nanocrystals: Bulk vs Surface

We show that the characteristic Mn2+ d emission color from Mn2+-doped CdS nanocrystals can be tuned over as much as 40 nm, in contrast to what should be expected from such a nearly localized d-d transition. This is achieved surprisingly by a fine-tuning of the host particle diameter from 1.9 to 2.6 nm, thereby changing the overall emission color from red to yellow. Systematic experiments in conjunction with state-of-the-art ab initio calculations with full geometry optimization establish that Mn2+ ions residing at surface/subsurface regions have a distorted tetrahedral coordination resulting in a larger ligand field splitting. Consequently, these near-surface Mn2+ species exhibit a lower Mn2+ d emission energy, compared to those residing at the core of the nanocrystal with an undisturbed tetrahedral coordination. The origin of the tunability of the observed Mn2+ emission is the variation of emission contributions arising from Mn2+ doped at the core, subsurface, and surface of the host. Our findings provide a unique and easy method to identify the location of an emitting Mn2+ ion in the nanocrystal, which would be otherwise very difficult to decipher.

The influence of lone-pair repulsions on C–C bond lengths: a critical evaluation of the experimental and theoretical evidence

X-Ray structural data, as well as semiempirical and ab initio molecular orbital calculations, reveal no systematic and substantial difference between the C–C bond lengths of cis and trans 1,2-diketones. Additional results on various conformations of 1,2-diimines and 1,2-dithiones follow the same pattern. Therefore, lone-pair repulsions cannot be implicated in the observed lengthening of C–C bonds in isatin and several related molecules. Conjugation in these systems occurs peripherally avoiding the participation of the central C–C bond. Negative hyperconjugative interaction between the oxygen lone pairs and the adjacent C–C σ* orbital is suggested to be the principal reason for the relatively long C–C bond in diketones. This effect is found in both the cis and trans conformations.

Patterning Nanoroads and Quantum Dots on Fluorinated Graphene

Using ab initio methods we have investigated the fluorination of graphene and find that different stoichiometric phases can be formed without a nucleation barrier, with the complete “2D-Teflon” CF phase being thermodynamically most stable. The fluorinated graphene is an insulator and turns out to be a perfect matrix-host for patterning nanoroads and quantum dots of pristine graphene. The electronic and magnetic properties of the nanoroads can be tuned by varying the edge orientation and width. The energy gaps between the highest occupied and lowest unoccupied molecular orbitals (HOMO-LUMO) of quantum dots are size-dependent and show a confinement typical of Dirac fermions. Furthermore, we study the effect of different basic coverage of F on graphene (with stoichiometries CF and C4F) on the band gaps, and show the suitability of these materials to host quantum dots of graphene with unique electronic properties.

A parametric method 3 (PM3) study of trans‐stilbene

We report a comparative modified neglect of diatomic overlap (MNDO), Austin method one (AM1), and parametric method 3 (PM3) study of trans‐stilbene (tS) in its ground, excited (singlet and triplet), and ionic (positive and negative polarons and bipolarons) states. We have also calculated the barrier for ring rotation about the backbone single bond. Our results show that PM3 geometries are superior to MNDO and AM1, at least for tS. PM3 predicts, in contrast with MNDO, AM1 and even ab initio 3‐21G, a coplanar structure for tS, in accordance with recent experimental data. Singlet and triplet energies obtained from heats of formation are in surprisingly good agreement with experimental data.