Modeling the interactions between MC2R and ACTH models from human

Melanocortin system is composed of four peptide hormones namely α-, β-, -γ, and adrenocorticotropic hormone (ACTH), derived from post-translational cleavage of a polypeptide precursor ‘proopiomelanocortin (POMC).’ Among these hormones, ACTH, a 38 amino acid residue peptide fragment is an important hormone as it is involved in steroid secretion. In addition to this, to cite a few, this hormone is also known to induce variety of other effects, such as alterations in motor/sexual behavior, improvement in memory, and anti-inflammatory effects. To date, five melanocortin receptors (MC1R–MC5R) have been characterized with tissue-specific expression patterns and different binding affinities for each of the melanocortin hormones to regulate various biological functions. In the present work, three-dimensional (3D) models of MC2R and ACTH from human have been predicted, followed by docking and molecular dynamics simulation. While the 3D model of MC2R receptor has been predicted through threading approach, structure of ACTH was built based on ab initio technique. The MC2R model was later successfully docked onto the ACTH structure. Molecular dynamics (MD) simulation for 20?ns was used to compute the binding free energy of MC2R with ACTH model under implicit solvent conditions.     

Effect of Solvation on Ozonolysis Reaction Intermediates and Transition States

Electrostatic solvent effects on the ozonolysis of ethylene have been investigated using correlated ab initio and density functional approaches. We use a simple polarizable continuum model for the solvent. It allows us to evaluate the medium effect on both the electronic and nuclear structure of the chemical species involved in the reaction. The computations confirm that basically the reaction proceeds through the Criegee mechanism. However, formation of the van der Waals complexes ethyl-ene/ozone and carbonyl oxide/formaldehyde also appears to play a role. All the calculated species are stabilized with respect to the reactants except the transition state corresponding to the primary ozonide formation. In general, electrostatic solvent effects are relatively small for activation barriers of single reaction steps and more substantial for the corresponding reaction energies. Moreover, the medium significantly modifies the structure of some species for which polarization effects are crucial.     

Modelling of Structural and Vibrational Properties of Poly(p-Phenylene) and Polypyrrole Using Molecular Orbital Methods

Abstract

This paper concentrates on two very important conducting polymers poly(p-phenylene) and polypyrrole. Detailed atomistic molecular models have been developed with the help of ab initio and semi-empirical quantum mechanical calculations using the Cerius² and WinMOPAC (version 6.0) programs.

Their optimised geometry had been calculated and compared with experimental X-ray diffraction data. The simulated and experimental vibrational spectra of biphenyl as well as isolated pyrrole monomers and oligomers from n = 1 and 2, where n is the number of structural repeat units used, have been computed using the ab initio 3–21G basis set. The results obtained are compared with experimental data for the case of biphenyl and for oligomers with n = 2 to 5 for both neutral benzenoid and quinonoid oligopyrroles, from semi-empirical predictions obtained by AM1 and PM3. The trends in the computed harmonic force fields, vibrational frequencies and intensities are monitored as a function of the chain length. The data are analyzed in conjunction with the trends in computed equilibrium geometries.     

Applicability of PM3 to transphosphorylation reaction path: Toward designing a minimal ribozyme

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at theab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on the system under investigation. In the preent study we assess the validity of the MNDO/PM3 method on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar toab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.     

Coordination and Thermodynamics of Stable Zn(II) Complexes in the Gas Phase

Abstract

Ab initio molecular orbital methods in combination with DFT calculations were used to study the structural and thermodynamic properties of 17 complexes containing zinc cation and four first-shell ligands as models of active site of metalloenzymes (e.g. angiotensin converting enzyme, thermolysin). The geometry of the complexes was relaxed by complete optimization by ab initio molecular orbital methods at Hertree-Fock level with 3–21G basis set. Following single point calculation with tight SCF criteria at the B3LYP level with 6–311+G(2d,p) basis set was used to calculate accurate interaction enthalpies. The structure and thermodynamics of optimized complexes are discussed from the point of view of their biological importance.     

Two-level parallelization of Ehrenfest force calculations in ab initio molecular dynamics simulation

A two-level parallel code for Ehrenfest force calculations in ab initio molecular dynamics simulations was developed for a shared memory multiprocessor cluster. Coarse-grain parallelism was implemented by atomic decomposition and a fine-grained parallelism was exploited to perform matrix multiplications. This two-level parallelism efficiently enhances the speed of computations.     

Recent advances in implicit solvent-based methods for biomolecular simulations

Implicit solvent-based methods play an increasingly important role in molecular modeling of biomolecular structure and dynamics. Recent methodological developments have mainly focused on the extension of the generalized Born (GB) formalism for variable dielectric environments and accurate treatment of nonpolar solvation. Extensive efforts in parameterization of GB models and implicit solvent force fields have enabled ab initio simulation of protein folding to native or near-native structures. Another exciting area that has benefited from the advances in implicit solvent models is the development of constant pH molecular dynamics methods, which have recently been applied to the calculations of protein pK(a) values and the studies of pH-dependent peptide and protein folding.     

Quadrupole Relaxation of the 7Li+ Ion in Dilute Aqueous Solution Determined by Experimental and Theoretical Methods

A combination of molecular dynamics simulations (MD), ab initio selfconsistent field (SCF) calculations and nuclear magnetic resonance relaxation time experiments (NMR) is a powerful battery of techniques to investigate the molecular origins of the nuclear quadrupole relaxation mechanism for ⁷Li⁺ ions in dilute aqueous solution.     

Towards the Understanding of the Folding of Methylene Units in the Glutamine Residue

The conformational preferences of the methylenic sequence in the side chain of the glutamine residue were investigated by ab initio and semi-empirical quantum mechanical calculations and examination of both the Brookhaven Protein Databank and Cambridge Structural Data Base. The results were analysed on the basis of our previous findings about the folding of methylene groups in aliphatic segments. Both energy calculations and the crystallographic structure of small peptides indicate that methylene units of the glutamine residue tend to fold in a gauche conformation. In contrast, such groups usually adopt an all-trans conformation in proteins due basically to the entropic and solvent contributions. These results have been demonstrated by computing the entropic correction to the free energy and evaluating the solvent effects through SCRF calculations     

Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.     

Local conformational variations observed in B-DNA crystals do not improve base stacking: computational analysis of base stacking in a d(CATGGGCCCATG)(2) B<-->A intermediate crystal structure

The crystal structure of d(CATGGGCCCATG)₂ shows unique stacking patterns of a stable B↔A-DNA intermediate. We evaluated intrinsic base stacking energies in this crystal structure using an ab initio quantum mechanical method. We found that all crystal base pair steps have stacking energies close to their values in the standard and crystal B-DNA geometries. Thus, naturally occurring stacking geometries were essentially isoenergetic while individual base pair steps differed substantially in the balance of intra-strand and inter-strand stacking terms. Also, relative dispersion, electrostatic and polarization contributions to the stability of different base pair steps were very sensitive to base composition and sequence context. A large stacking flexibility is most apparent for the CpA step, while the GpG step is characterized by weak intra-strand stacking. Hydration effects were estimated using the Langevin dipoles solvation model. These calculations showed that an aqueous environment efficiently compensates for electrostatic stacking contributions. Finally, we have carried out explicit solvent molecular dynamics simulation of the d(CATGGGCCCATG)₂ duplex in water. Here the DNA conformation did not retain the initial crystal geometry, but moved from the B↔A intermediate towards the B-DNA structure. The base stacking energy improved in the course of this simulation. Our findings indicate that intrinsic base stacking interactions are not sufficient to stabilize the local conformational variations in crystals.     

Structure and dynamics of human vimentin intermediate filament dimer and tetramer in explicit and implicit solvent models

Intermediate filaments, in addition to microtubules and microfilaments, are one of the three major components of the cytoskeleton in eukaryotic cells, and play an important role in mechanotransduction as well as in providing mechanical stability to cells at large stretch. The molecular structures, mechanical and dynamical properties of the intermediate filament basic building blocks, the dimer and the tetramer, however, have remained elusive due to persistent experimental challenges owing to the large size and fibrillar geometry of this protein. We have recently reported an atomistic-level model of the human vimentin dimer and tetramer, obtained through a bottom-up approach based on structural optimization via molecular simulation based on an implicit solvent model (Qin et al. in PLoS ONE 2009 4(10):e7294, 9). Here we present extensive simulations and structural analyses of the model based on ultra large-scale atomistic-level simulations in an explicit solvent model, with system sizes exceeding 500,000 atoms and simulations carried out at 20 ns time-scales. We report a detailed comparison of the structural and dynamical behavior of this large biomolecular model with implicit and explicit solvent models. Our simulations confirm the stability of the molecular model and provide insight into the dynamical properties of the dimer and tetramer. Specifically, our simulations reveal a heterogeneous distribution of the bending stiffness along the molecular axis with the formation of rather soft and highly flexible hinge-like regions defined by non-alpha-helical linker domains. We report a comparison of Ramachandran maps and the solvent accessible surface area between implicit and explicit solvent models, and compute the persistence length of the dimer and tetramer structure of vimentin intermediate filaments for various subdomains of the protein. Our simulations provide detailed insight into the dynamical properties of the vimentin dimer and tetramer intermediate filament building blocks, which may guide the development of novel coarse-grained models of intermediate filaments, and could also help in understanding assembly mechanisms.     

Methane in Water: An Ab Initio Study

Abstract

In preceding publications we discussed some properties of pure water in condensed phases using an ab initio approach. Here this study is used as a basis of comparison for analysing the behaviour of water as a solvent in the presence of an apolar molecule. Our analysis is focused on the process of organization of the hydrogen bonding network around the solute. For this purpose we perform some ab initio calculations for a system of 32 water molecules and one methane molecule at 300 K; in particular, the average molecular dipole moment of water is determined and the result is compared with that of pure water. Next the attention is switched to the methane molecule; related properties such as excluded volume and sphericity of its shape are illustrated and discussed. A comparison with results obtained using classical approaches suggests that some classical models of water can be considered to be still valid when they are used to analyse the water-methane system.     

A force field for naproxen

A united-atom potential model for naproxen suitable for molecular dynamics (MD) simulation has been developed. The charge distribution is approximated by point charges obtained from ab initio calculations using the CHELPG method. Also the intramolecular interactions such as bond and angle vibration, and the torsion potential are obtained from ab initio calculations. The dispersive interaction contribution is taken from the literature. By MD simulation using a naproxen film in slap geometry, the temperature dependence of the density, surface tension and self-diffusion coefficient as well as the melting temperature for the developed potential model are obtained.     

Hydroxy-1,2,5-oxadiazolyl moiety as bioisoster of the carboxy function. A computational study on γ-aminobutyric acid (GABA) related compounds

Recently, our research group has proposed the hydroxyfurazanyl (4-hydroxy-1,2,5-oxadiazole-3-yl) moiety as a new non-classical isoster of the carboxy function in the design of γ-aminobutyric acid (GABA) analogues. Some compounds showed significant activity at the GABA_A receptor, representing the only examples of pentatomic heterocycles bearing an ω-aminoalkyl flexible side chain in the position vicinal to the hydroxy group displaying agonist activity at this receptor subtype. In this work, an ab initio analysis of the structural and electronic features of furazan-3-ol is presented, in order to provide a theoretical basis to the claimed bioisosterism with the carboxy function. An ab initio conformational study with the C-PCM implicit solvent model was carried out to elucidate the reasons of the peculiar behaviour of the furazan models. Alongside, another conformational search through molecular dynamics in explicit solvent was accomplished, in order to validate the first method. The electronic features of the 4-hydroxy-1,2,5-oxadiazole-3-yl substructure seem to account for a marked stabilising effect of the putative bioactive conformation at the GABA_A receptor subtype. The 1,2,5-thiadiazole analogue, which shares the same conformational preference of its oxygenated counterpart, was identified as a potential candidate for synthesis and pharmacological testing. Figure 4-(ω-aminoalkyl)-1,2,5-oxadiazole-3-ol analogues of GABA Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorised users.     

Hybrid MC/QC simulations of water-assisted proton transfer in nucleosides. Guanosine and its analog acyclovir

To provide an in-depth insight into the molecular basis of spontaneous tautomerism in DNA and RNA base pairs, a hybrid Monte Carlo (MC)–quantum chemical (QC) methodology is implemented to map two-dimensional potential energy surfaces along the reaction coordinates of solvent-assisted proton transfer processes in guanosine and its analog acyclovir in aqueous solution. The solvent effects were simulated by explicit inclusion of water molecules that model the relevant part of the first hydration shell around the solute. The position of these water molecules was estimated by carrying out a classical Metropolis Monte Carlo simulation of dilute water solutions of the guanosine (Gs) and acyclovir (ACV) and subsequently analyzing solute–solvent intermolecular interactions in the statistically-independent MC-generated configurations. The solvent-assisted proton transfer processes were further investigated using two different ab initio MP2 quantum chemical approaches. In the first one, potential energy surfaces of the ‘bare’ finite solute–solvent clusters containing Gs/ACV and four water molecules (MP2/6-31+G(d,p) level) were explored, while within the second approach, these clusters were embedded in ‘bulk’ solvent treated as polarizable continuum (C-PCM/MP2/6-31+G(d,p) level of theory). It was found that in the gas phase and in water solution, the most stable tautomer for guanosine and acyclovir is the 1H-2-amino-6-oxo form followed by the 2-amino-6-(sZ)-hydroxy form. The energy barriers of the water-assisted proton transfer reaction in guanosine and in acyclovir are found to be very similar – 11.74 kcal mol^?1 for guanosine and 11.16 kcal mol^?1 for acyclovir, and the respective rate constants (k = 1.5?×?10¹ s^?1, guanosine and k = 4.09?×?10¹ s^?1, acyclovir), are sufficiently large to generate the 2-amino-6-(sZ)-hydroxy tautomer. The analysis of the reaction profiles in both compounds shows that the proton transfer processes occur through the asynchronous concerted mechanism.     

Accurate prediction of protein secondary structure and solvent accessibility by consensus combiners of sequence and structure information

Background  

Structural properties of proteins such as secondary structure and solvent accessibility contribute to three-dimensional structure prediction, not only in the ab initio case but also when homology information to known structures is available. Structural properties are also routinely used in protein analysis even when homology is available, largely because homology modelling is lower throughput than, say, secondary structure prediction. Nonetheless, predictors of secondary structure and solvent accessibility are virtually always ab initio.     

Implicit solvent models for flexible protein-protein docking by molecular dynamics simulation

The suitability of three implicit solvent models for flexible protein-protein docking by procedures using molecular dynamics simulation is investigated. The three models are (i) the generalized Born (GB) model implemented in the program AMBER6.0; (ii) a distance-dependent dielectric (DDD) model; and (iii) a surface area-dependent model that we have parameterized and call the NPSA model. This is a distance-dependent dielectric model modified by neutralizing the ionizable side-chains and adding a surface area-dependent solvation term. These solvent models were first tested in molecular dynamics simulations at 300 K of the native structures of barnase, barstar, segment B1 of protein G, and three WW domains. These protein structures display a range of secondary structure contents and stabilities. Then, to investigate the performance of the implicit solvent models in protein docking, molecular dynamics simulations of barnase/barstar complexation, as well as PIN1 WW domain/peptide complexation, were conducted, starting from separated unbound structures. The simulations show that the NPSA model has significant advantages over the DDD and GB models in maintaining the native structures of the proteins and providing more accurate docked complexes.