Abstract

In this paper, we have computationally described the interaction of yDNA bases with small gold clusters through a variety of aspects including geometrical, spectroscopic and energetic one. To get insights in to essential elements of this interaction, the contribution of typical energy interaction components have also been investigated.

Abstract

Controlling dynamics at the atomic level is a primary goal of chemists. In the last decade, Optimal control theory (OCT) has emerged as one comprehensive means of designing laser pulses for achieving prescribed dynamical goals. We have used OCT to obtain infrared laser pulses for the selective vibrational excitation of several different systems, from a simple diatomic to triatomic systems located in complex molecular environments. We will present our results obtained for the diatomic system HF applying OCT using different algorithmic methods, namely an iterative method, conjugate gradient method and genetic algorithms. The results of a new dynamic method for applying the penalty term needed to keep the power in the electric field of the laser within reasonable bounds will be discussed. For the triatomic FeCO, we use a two mathematical dimensional model of carboxy-myoglobin. Density functional theory is used to obtain the potential energy and dipole moment surfaces of the active site model. The Conjugate gradient method is employed to optimize the cost functional and to obtain the optimized laser pulses. Optimized laser fields are found which give virtually 100% excitation probability to preselected vibrational levels. Similar application to hydrogen bonding in a complex system will also be presented.

Abstract

The synthesis of modified versions of DNA is an area that is receiving much attention because of their potential applications in biotechnology and nanotechnology. Recently, Eric Kool et al have re-engineered an expanded version of DNA by pairing size expanded x-bases with normal bases in a complementary manner, which they called as xDNA. xDNA bases are analogues of natural bases that are formed by insertion of benzene ring into the natural bases. Recently, there was arisen a huge interest in conducting properties of metal bound DNA structures. We use quantum chemistry methods to investigate the structural and electronic properties of x-bases and their interaction with gold atoms, in order to explore develop insight into the nature of binding between gold atoms and size expanded DNA.

Abstract

The problem of predicting the 3D native conformation of a polypeptide given the primary sequence is of fundamental importance in contemporary biology. Here we present a novel Genetic algorithm method to obtain the optimal geometry of a pentapeptide. The pentapeptide moiety is found in wide set of configurations in protein databases. A genetic algorithm utilizes optimization procedure similar to natural genetic evolution. Unlike other genetic algorithms, which are constrained to operate on bit strings, the algorithm used here operates on real values. A random population of conformations is generated taking into account the constraints of Ramachandran plot. This population advances through successive generations in which the solutions evolve via genetic operators (Mutate, Variate and Crossover); in such a way that it gives both better survival chance to fitter soluctions and a larger diversity within the population. More attention will be given to favorable local structures while unfavorable local structures will be rapidly abandoned on the basis of their fitness functions (potential energy), which are calculated using the GROMOS-96 force field of Spdbv version-3.7. The SCWRL3.0 program developed at Dunbrack lab is used for placing side chains to a fixed backbone coordinates. This algorithm proved to be very efficient for the penta-peptide Met-Enkephalin, a natural Opioid polypeptide produced in the brain. The algorithm converged only in 41 minutes and in 23 generations keeping the population size fixed (20 individual per population). The converged structure is energetically as stable as the native conformation of the protein. Its shows a RMSD of 5.9 from the original native conformation.

Abstract

With availability of fast computers, electronic structure and properties calculations of small to medium sized systems using quantum chemical methods are now commonplace. Several options in addition to obtaining structures with energy optimisation are now easily realisable, e.g., carrying out calculations of electrostatic potential interations, atoms in molecule approach of Bader, Natural Bonding orbital method, etc. We discuss results obtained in specific cases using these options. We have recently carried out studies on several such systems including amino-methyl furan carboxylic acid system and binding of gold nanoclusters with size expanded DNA bases, xA, xC, xG, and xT. We will present the motivation behind such works and discuss the results obtained using different methods. Geometries of these molecular systems were fully optimized using the Hartree Fock SCF method and using B3LYP density functional method (DFT). The calculations provide a theoretical understanding beyond a phenomenological observation of the preferential cyclotrimerization of 5- (aminomethyl)-2-furancarboxylic acid (AMFC). Other than calculations of thermodynamic free energy, a detailed analysis of the molecular properties is carried out. The electrostatic potentials (ESP) in the region around molecular skeletons are investigated to provide justication for the preferential formation of the cyclic tripeptide over the cyclic dipeptide from the building block, AMFC. The ESP are calculated and computed on the molecular surfaces of the oligopeptides. There are repulsions observed in the dimer due to the proximity of the ethereal oxygen atoms of the furan rings. Natural bond orbital (NBO) analysis is done to find the second-order interactions, and atoms in molecules (AIM) calculations are performed to explore the interactions between the donor and acceptor moieties in the molecule. Our results on the gold complexed x-bases show that the gold clusters around xA and xT adopt triangular geometries, whereas irregular structures are obtained in the case of gold clusters complexed around xC and xG. The lengths of the bonds between atoms in the x-bases increase on gold complexation. The aromaticcharacter of the x-bases also increases on gold complexation except for the five member rings. A significant charge transfer from the x-base to gold atoms is seen in these complexes. Second-order interactions are observed in addition to direct covalent bonds between gold atoms and x bases.

Abstract

The last two decades have seen increasing interest in the area of controlling dynamics at the atomic level. Among the several approaches used, Optimal Control Theory (OCT) has emerged as a powerful method for designing laser pulses to achieve prescribed dynamical goals. We have used OCT to obtain infrared laser pulses for the selective vibrational excitation of several different systems, from a simple diatomic to triatomics located in complex molecular environments. We use time dependent quantum mechanics to monitor the evolution of the system subject to a varying laser field. We have obtained fields for vibrational control of the diatomic system HF applying OCT using different algorithms, namely an iterative method, conjugate gradient method and genetic algorithms. The time-dependent quantum mechanical propagation of the system is achieved using the split operator form of the time evolution operator together with Fast Fourier transforms and grid methods. For the modelling of FeCO in carboxy-myoglobin, we use a two mathematical dimensional model. Density functional theory is used to obtain the potential energy and dipole moment surfaces of the active site model, while the conjugate gradient method is employed within the OCT optimization procedure to obtain the laser pulses. Optimized laser fields are found which give virtually 100% excitation probability to preselected vibrational levels. A similar application to a hydrogen bonded system of biological significance will also be presented. In all cases, detailed analyses of the field are carried out in terms of its power spectrum and the time variation of the contributions from the frequencies involved.

Abstract

The backbone conformation of amino acids in proteins is represented mainly by two torsion angles; phi and psi. Ramachandran map is a way to visualize these torsion angles in protein structures. Here we analyze the nature and distribution of disallowed Ramachandran conformations of amino acid residues observed in protein structures. A non redundant data set consisting of 405 high resolution protein crystal structures from the Brookhaven Protein Data Bank was examined. The data set consisted of a total of 78,213 non-Gly residues, which were characterized on the basis of their backbone dihedral angles phi and psi. Residues falling outside the defined "broad allowed limits" on the Ramachandran map were analyzed. The results suggested that very small fraction of residues (approximately 0.4%) lie in appreciably disallowed regions. Interestingly, small polar/charged residues showed significantly greater probability of adopting unusual backbone with the largest number of examples being observed for Asn, Asp, Ser and Thr residues. The bulky charged residues Glu, Lys, Gln and Arg have a relatively low propensity for backbone distortions. Bulky hydrophobic residues which are found predominantly in well packed interiors of proteins like Ile, Val and aromatic amino acids like Phe, Trp and Tyr do not generally adopt disallowed conformations. In the allowed region irrespective of the type of amino acid except Gly, there is a high propensity for phi value of the range -70 to -60 and for psi value 140 to 150.

Abstract

Depth movies provide 2-D representations of a time-varying 3D scene. Multistream depth movies arise in many structure-capturing setups and have been used for image based rendering. They are bulky and carry redundant information. We propose a proxy-based compression scheme for multistream depth movies of a scene involving dynamic human actors. The input to our system is depth movies from different viewpoints and the calibration parameters to relate depths to 3D points. We use an articulated human model as a proxy to represent the common structure of the scene. The proxy is parametrized by various bone angles.

Abstract

Graph cuts has become a powerful and popular opti-mization tool for energies defined over an MRF and havefound applications in image segmentation, stereo vision,image restoration, etc. The maxflow/mincut algorithm tocompute graph-cuts is computationally heavy. The best-reported implementation of graph cuts takes over 100 mil-liseconds even on images of size640×480and cannot beused for real-time applications or when iterated applica-tions are needed. The commodity Graphics Processor Unit(GPU) has emerged as an economical and fast computationco-processor recently. In this paper, we present an imple-mentation of the push-relabel algorithm for graph cuts onthe GPU. We can perform over 60 graph cuts per secondon1024×1024images and over 150 graph cuts per secondon640×480images on an Nvidia 8800 GTX. The time foreach complete graph-cut is about 1 millisecond when only afew weights change from the previous graph, as on dynamicgraphs resulting from videos. The CUDA code with a well-defined interface can be downloaded for anyone’s use

Abstract

Capturing dynamic scenes in 3D using a suitable depth or structure recovery mechanism has become popular today for image based modelling, 3D telelconferencing, etc. The 2 1 2 D geometric structure and the aligned texture – together called a depth image – are captured by such setups. Time varying sequences of depth images are called depth movies. Depth movies of humans performing interesting actions are being captured using specialized setups today in different research labs. The data so captured is often streamed to another location for immersive viewing using standard depth-image rendering techniques. The depth maps are bulky and need innovative algorithms for efficient representation and compression. In this paper, we present a compression scheme that uses parametric proxy models for the underlying action. We use a generic articulated human model as the proxy and the various joint angles as its parameters for each time instant. The proxy model represents common prediction of the scene structure for that time instant. It can be projected to each view and the difference or residue between the depth due to the proxy and the original depth can be used to represent the scene. The residues compress better and can exploit temporal and spatial relationship between multiple depth movies. We show results on several synthetic actions to demonstrate the usefulness of the scheme