Docking and free energy simulations to predict conformational domains involved in hCG-LH receptor interactions using recombinant antibodies

Single chain fragment variables (ScFvs) have been extensively employed in studying the protein-protein interactions. ScFvs derived from phage display libraries have an additional advantage of being generated against a native antigen, circumventing loss of information on conformational epitopes. In the present study, an attempt has been made to elucidate human chorionic gonadotropin (hCG)-luteinizing hormone (LH) receptor interactions by using a neutral and two inhibitory ScFvs against hCG. The objective was to dock a computationally derived model of these ScFvs onto the crystal structure of hCG and understand the differential roles of the mapped epitopes in hCG-LH receptor interactions. An anti-hCG ScFv, whose epitope was mapped previously using biochemical tools, served as the positive control for assessing the quality of docking analysis. To evaluate the role of specific side chains at the hCG-ScFv interface, binding free energy as well as residue interaction energies of complexes in solution were calculated using molecular mechanics Poisson-Boltzmann/surface area method after performing the molecular dynamic simulations on the selected hCG-ScFv models and validated using biochemical and SPR analysis. The robustness of these calculations was demonstrated by comparing the theoretically determined binding energies with the experimentally obtained kinetic parameters for hCG-ScFv complexes. Superimposition of hCG-ScFv model onto a model of hCG complexed with the 51-266 residues of LH receptor revealed importance of the residues previously thought to be unimportant for hormone binding and response. This analysis provides an alternate tool for understanding the structure-function analysis of ligand-receptor interactions. Proteins 2011;79:3108-3122. (C) 2011 Wiley-Liss, Inc.

Self-Assembly of Metallamacrocycles Using a Dinuclear Organometallic Acceptor: Synthesis, Characterization, and Sensing Study

A dinuclear organometallic acceptor 4,4'-bis[trans-Pt(PEt(3))(2)(O(3)SCF(3))(ethynyl)]biphenyl (1) containing Pt-ethynyl functionality is synthesized. Multinuclear NMR ((1)H, (31)P, and (13)C), infrared (IR), and electrospray ionization mass spectrometry (ESI-MS) including single-crystal X-ray diffraction analysis established the formation of 1. Equimolar treatment of acceptor 1 separately with three different ``clip'' type ditopic donors (L(a)-L(c)) yielded [2 + 2] self-assembled three metallamacrocycles 2a-2c, respectively. These macrocycles were characterized by various spectroscopic techniques, and their sizes/shapes were obtained through geometry optimization using molecular mechanics universal force field (MMUFF) simulations. Attachment of unsaturated ethynyl functionality to biphenyl building unit helped to make the macrocycles (2a-2c) pi-electron rich and thereby fluorescent in nature. Furthermore, 2c in solution has been examined to be suitable for sensing electron-deficient nitroaromatic like picric acid, which is often considered as a secondary chemical explosive. The fluorescence study of 2c showed a marked quenching of initial emission intensity upon titrating with picric acid (PA), and it exhibited the largest fluorescence quenching response with high selectivity among various other electron deficient aromatic compounds tested.

Pt-6(II) nanoscopic cages with an organometallic backbone as sensors for picric acid

An organometallic building block 1,3,5-tris(4-trans-Pt(PEt3)(2)I(ethynyl)phenyl)benzene (1) incorporating Pt-ethynyl functionality has been synthesized and characterized. 2 + 3] self-assembly of its nitrate analogue 1,3,5-tris(4-trans-Pt(PEt3)(2)(ONO2)(ethynyl)phenyl)benzene (2) with ``clip'' type bidentate donors (L-1-L-3) separately afforded three trigonal prismatic architectures (3a-3c), respectively. All these prisms were characterized and their shapes/sizes are predicted through geometry optimization employing molecular mechanics universal force field (MMUFF) simulation. The extended p-conjugation including the presence of Pt-ethynyl functionality makes them electron rich as well as luminescent in nature. Macrocycles 3b and 3c exhibit fluorescence quenching in solution upon addition of picric acid PA], which is a common constituent of many explosives. Interestingly, the non-responsive nature of fluorescent intensity towards other electron-deficient nitro-aromatic explosives (NAEs) makes them promising selective sensors for PA with a detection limit predicted to be ppb level. Furthermore, solid-state quenching of fluorescent intensity of the thin film of 3b upon exposure to saturated vapor of picric acid has drawn special attention for infield applications.

VIBRATIONAL CHARACTERISTICS OF ZIGZAG, ARMCHAIR AND CHIRAL CANTILEVER SINGLE-WALLED CARBON NANOTUBES

Finite element analysis has been performed to study vibrational characteristics of cantilever single walled carbon nanotubes. Finite element models are generated by specifying the C-C bond rigidities, which are estimated by equating energies from molecular mechanics and continuum mechanics. Bending, torsion, and axial modes are identified based on effective mass for armchair, zigzag and chiral cantilever single walled carbon nanotubes, whose Young's modulus is evaluated from the bending frequency. Empirical relations are provided for frequencies of bending, torsion, and axial modes.

A novel method for studying the buckling of nanotubes considering geometrical imperfections

Buckling of nanotubes has been studied using many methods such as molecular dynamics (MD), molecular mechanics, and continuum-based shell theories. In MD, motion of the individual atoms is tracked under applied temperature and pressure, ensuring a reliable estimate of the material response. The response thus simulated varies for individual nanotubes and is only as accurate as the force field used to model the atomic interactions. On the other hand, there exists a rich literature on the understanding of continuum mechanics-based shell theories. Based on the observations on the behavior of nanotubes, there have been a number of shell theory-based approaches to study the buckling of nanotubes. Although some of these methods yield a reasonable estimate of the buckling stress, investigation and comparison of buckled mode shapes obtained from continuum analysis and MD are sparse. Previous studies show that the direct application of shell theories to study nanotube buckling often leads to erroneous results. The present study reveals that a major source of this error can be attributed to the departure of the shape of the nanotube from a perfect cylindrical shell. Analogous to the shell buckling in the macro-scale, in this work, the nanotube is modeled as a thin-shell with initial imperfection. Then, a nonlinear buckling analysis is carried out using the Riks method. It is observed that this proposed approach yields significantly improved estimate of the buckling stress and mode shapes. It is also shown that the present method can account for the variation of buckling stress as a function of the temperature considered. Hence, this can prove to be a robust method for a continuum analysis of nanosystems taking in the effect of variation of temperature as well.

Electronically nonadiabatic decomposition mechanisms of clusters of zinc and dimethylnitramine

Electronically nonadiabatic decomposition mechanisms of dimethylnitramine (DMNA) in presence of zinc metal clusters are explored. Complete active space self-consistent field (CASSCF) calculation is employed for DMNA-Zn and ONIOM (Our own N-layered integrated molecular orbital and molecular mechanics) methodology is coupled with CASSCF methodology for DMNA-Zn-10 cluster. Present computational results show that DMNA-Zn clusters undergo electronically nonadiabatic reactions, rendering nitro-nitrite isomerization followed by NO elimination. The overall reactions are also found to be highly exothermic in nature. This is the first report on electronically nonadiabatic decomposition pathways of DMNA-Zn-n neutral clusters. (C) 2014 Elsevier B.V. All rights reserved.

Template-free coordination-driven self-assembly of discrete hexanuclear prismatic cages employing half-sandwich octahedral Ru-2(II) acceptors and triimidazole donors

Coordination-driven self-assembly of dinuclear half-sandwich p-cymene ruthenium(II) complexes Ru-2(mu-eta(4)-C2O4)(CH3OH)(2)(eta(6)-p-cymene)(2)](O3SCF3)(2) (1a) and Ru-2(mu-eta(4)-C6H2O4)(CH3OH)(2)(eta(6)-p-cymene)(2)](O3SCF3)(2) (1b) separately with imidazole-based tritopic donors (L-1-L-2) in methanol yielded a series of hexanuclear 3+2] trigonal prismatic cages (2-5), respectively L-1 = 1,3,5-tris(imidazole-1-yl) benzene; L-2 = 4,4',4 `'-tris(imidazole-1-yl) triphenylamine]. All the self-assembled cages 2-5 were characterized by various spectroscopic techniques (multinuclear NMR, Infra-red and ESI-MS) and their sizes, shapes were obtained through geometry optimization using molecular mechanics universal force field (MMUFF) computation. Despite the possibility due to the free rotation of donor sites of imidazole ligands, of two different atropoisomeric prismatic cages (C-3h or C-s) and polymeric product, the self-selection of single (C(3)h) conformational isomeric cages as the only product is a noteworthy observation. (C) 2015 Elsevier B.V. All rights reserved.

Evaluating boundary dependent errors in QM/MM simulations.

Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations provide a powerful tool for studying chemical reactions, especially in complex biochemical systems. In most works to date, the quantum region is kept fixed throughout the simulation and is defined in an ad hoc way based on chemical intuition and available computational resources. The simulation errors associated with a given choice of the quantum region are, however, rarely assessed in a systematic manner. Here we study the dependence of two relevant quantities on the QM region size: the force error at the center of the QM region and the free energy of a proton transfer reaction. Taking lysozyme as our model system, we find that in an apolar region the average force error rapidly decreases with increasing QM region size. In contrast, the average force error at the polar active site is considerably higher, exhibits large oscillations and decreases more slowly, and may not fall below acceptable limits even for a quantum region radius of 9.0 A. Although computation of free energies could only be afforded until 6.0 A, results were found to change considerably within these limits. These errors demonstrate that the results of QM/MM calculations are heavily affected by the definition of the QM region (not only its size), and a convergence test is proposed to be a part of setting up QM/MM simulations.

Theoretical studies on conformation comparison of braid-like and triplex DNA

Based on the experimental data of scanning tunneling microscopy (STM), models of three-stranded braid-like DNAs composed by three kinds of base triplets AAA, TAT and GCA were constructed. We investigated the braid-like DNAs and their comparative tripler DNAs using a molecular mechanics method. The three strands of braid-like DNAs are proven equivalent, while those of tripler DNAs are not. The conformational energies for braid-like DNAs were found to be higher than that for tripler DNAs. Each period in one strand of braid-like DNA has 18 nucleotides, half of which are right-handed, while the other half are left-handed. Additional discussions concerning sugar puckering modes and the H-bonds are also included. (C) 1999 Elsevier Science B.V. All rights reserved.

Tests of an adaptive QM/MM calculation on free energy profiles of chemical reactions in solution.

We present reaction free energy calculations using the adaptive buffered force mixing quantum mechanics/molecular mechanics (bf-QM/MM) method. The bf-QM/MM method combines nonadaptive electrostatic embedding QM/MM calculations with extended and reduced QM regions to calculate accurate forces on all atoms, which can be used in free energy calculation methods that require only the forces and not the energy. We calculate the free energy profiles of two reactions in aqueous solution: the nucleophilic substitution reaction of methyl chloride with a chloride anion and the deprotonation reaction of the tyrosine side chain. We validate the bf-QM/MM method against a full QM simulation, and show that it correctly reproduces both geometrical properties and free energy profiles of the QM model, while the electrostatic embedding QM/MM method using a static QM region comprising only the solute is unable to do so. The bf-QM/MM method is not explicitly dependent on the details of the QM and MM methods, so long as it is possible to compute QM forces in a small region and MM forces in the rest of the system, as in a conventional QM/MM calculation. It is simple, with only a few parameters needed to control the QM calculation sizes, and allows (but does not require) a varying and adapting QM region which is necessary for simulating solutions.

Interface energy and its influence on interface fracture between metal and ceramic thin films in nanoscale

A theoretical model about the size-dependent interface energy between two thin films with different materials is developed by considering the chemical bonding contribution based on the thermodynamic expressions and the structure strain contribution based on the mechanical characteristics. The interface energy decreases with reducing thickness of thin films, and is determined by such available thermodynamic and mechanical parameters as the melting entropy, the melting enthalpy, the shear modulus of two materials, etc. The predicted interface energies of some metal/MgO and metal/Al2O3 interfaces based on the model are consistent with the results based on the molecular mechanics calculation. Furthermore, the interface fracture properties of Ag/MgO and Ni/Al2O3 based on the atomistic simulation are further compared with each other. The fracture strength and the toughness of the interface with the smaller structure interface energy are both found to be lower. The intrinsic relations among the interface energy, the interface strength, and the fracture toughness are discussed by introducing the related interface potential and the interface stress. The microscopic interface fracture toughness is found to equal the structure interface energy in nanoscale, and the microscopic fracture strength is proportional to the fracture toughness. (C) 2010 American Institute of Physics. [doi:10.1063/1.3501090]

Study of adsorption of methylene blue and new methylene blue in liquid-solid interface by slab optical waveguide spectroscopy

A slab optical waveguide (SOWG) has been used for study of adsorption of both methylene blue (MB) and new methylene blue (NMB) in liquid-solid interface. Adsorption characteristics of MB and NMB on both bare SOWG and silanized SOWG by octadecyltrichlorosilane (ODS) were compared. Effect of pH on adsorption on MB and NMB was investigated. Binding rate constant analysis showed that both MB and NMB on bare SOWG demonstrates larger association constants than those on ODS-SOWG. Interactions of NIB and NMB on bare SOWG and ODS-SOWG were analyzed by molecular mechanics calculation method. The binding energy change was in the following order: ENMB-bare > EMB-bare > ENMB-ODS > EMB-ODS. (c) 2004 Elsevier B.V. All rights reserved.

Novel inclusion complex of ibuprofen tromethamine with cyclodextrins: Physico-chemical characterization

Guest-host interactions of ibuprofen tromethamine salt (Ibu.T) with native and modified cyclodextrins (CyDs) have been investigated using several techniques, namely phase solubility diagrams (PSDs), proton nuclear magnetic resonance (H-1 NMR), differential scanning calorimetry (DSC), Fourier transform infrared spectroscopy (FT-IR), X-ray powder diffractometry (XRPD). scanning-electron microscopy (SEM) and molecular mechanics (MM). From the analysis of PSD data (A(L)-type) it is concluded that the anionic tromethamine salt of ibuprofen (pK(a) = 4.55) forms 1: 1 soluble complexes with all CyDs investigated in buffered water at pH 7.0, while the neutral form of Ibu forms an insoluble complex with beta-CyD (B-S-type) in buffered water at pH 2.0. Ibu.T has a lower tendency to complex with beta-CyD (K-11 = 58 M-1 at pH 7.0) compared with the neutral Ibu (K-11 = 4200 M (1)) in water. Complex formation of Ibu.T with beta-CyD (Delta G degrees = -20.4 kJ/mol) is enthalpy driven (Delta H degrees = -22.9 kJ/mol) and is accompanied by a small unfavorable entropy (Delta S degrees = -8.4 J/mol K) change. H-1 NMR studies and MM computations revealed that, on complexation, the hydrophobic central benzene ring of lbu.T and part of the isobutyl group reside within the beta-CyD cavity leaving the peripheral groups (carboxylate, tromethamine and methyl groups) located near the hydroxyl group networks at either rim of beta-CyD. PSD, H-1 NMR, DSC, FT-IR, XRPD, SEM and MM studies confirmed the formation of Ibu.T/beta-CyD inclusion complex in solution and the solid state. (C) 2009 Elsevier B.V. All rights reserved.

The Effect of Tacticity on the Conformational Properties of Poly(1-olefin sulfone)s

The effect of tacticity on the conformational properties of poly(olefin sulfone)s was studied. Tactic polymers, prepared from racemic thiirane monomers using chiral inititators were compared with atactic polymers prepared by free radical co-polymerisation of the 1-olefin with sulfur dioxide. Analysis of the XRD patterns showed that the tactic polymers formed more ordered structures in the bulk with longer layer spacings, consistent with a model in which their side chains meet at the tips in contrast with the atactic polymers whose side chains interdigitate. 13C MAS nmr experiments suggest that as tacticity increases so too does the proportion of C-S bonds in the gauche conformation, however the proportion of S-C bonds in the trans conformation falls, in contrast to a reported molecular mechanics study. Finally, DSC measurements on the polymers with longer side chains showed the presence of two endotherms on heating, illustrating definite liquid crystalline behaviour.

Suggested Binding Mechanism of the HIV-gp120 to its CD4 Receptor

The molecular recognition and attachment of the CD4 molecule and the HIV envelope glycoprotein (gp120) might be described as a consecutive three-step molecular recognition process. 1. (a) Long range interaction: electrostatic pre-orientation, 2. (b) short range interaction: electronic attachment followed by a ‘Locking-in’ (via aromatic ring orientation) and 3. (c) internal interaction (induced fit): conformational readjustment of the protein molecules. On the basis of the preliminary investigations (X-ray structures of CD4 and biological studies of CD4 and gp120 point mutants) we described a computational model. This approach consists of empirical calculations as well as ab initio level of quantum chemistry. The conformational analysis of the wild type and mutant CD4 molecules was supported by molecular mechanics and dynamics (Amber force field). The latter analysis involves the application of a novel method, the Amino Acid Conformation Assignment of Proteins (ACAP) software, developed for the notation of secondary protein structures. According to the cardinal role of the electrostatic factors during this interaction, several ab initio investigations were performed for better understanding of the recognition process on submolecular level. Using the above mentioned computational model, we could interpret the basic behaviours and predict some additional features of CD4-gp120 interaction, in spite of the missing gp120 X-ray structure.