PubMed Health⌕ Search

Biomedical subjects

Stefano Piana

Publications and source records attributed to Stefano Piana.

13 recordsLinked to original sources
PubMed2006

The nature of the adsorption of nucleobases on the gold [111] surface.

Biochip technology is based on the immobilization of biological macromolecules on the surface of electronic devices. The biochemical properties of the immobilized molecules can be influenced to a large extent by the interaction with the inorganic surface. The interaction of DNA with the surface of gold, a metal commonly used in biochip technologies, is sequence dependent as the four nucleobases, adenine, thymine, cytosine, and guanine, interact to a different extent with the gold surface. The nature of nucleobase adsorption on the gold [111] surface has been investigated by performing density functional theory and post-Hartree-Fock calculations. It turns out that the interaction is dominated by dispersion forces and an appreciable degree of chemisorption is observed for adenine only. A set of Lennard-Jones parameters that describe the interaction was derived from the post-Hartree-Fock calculations. Classical molecular dynamics simulations of nucleobase monolayers based on these parameters are in remarkable agreement with the experiment and show that the interaction of the nucleobases with the gold surface is strongly modulated by base-base interactions and reaches a maximum when a full monolayer is formed.

Adenine↗
PubMed2006

Assisted desolvation as a key kinetic step for crystal growth.

The crystallization of materials from a supersaturated solution is a fundamental chemical process. Although several very successful models that provide a qualitative understanding of the crystal growth process exist, in most cases the atomistic detail of crystal growth is not fully understood. In this work, molecular dynamics simulations of the morphologically most important surfaces of barite in contact with a supersaturated solution have been performed. The simulations show that an ordered and tightly bound layer of water molecules is present on the crystal surface. The approach of an ion to the surface requires desolvation of both the surface and the ion itself leading to an activated process that is rate limiting for two-dimensional nucleation to occur. However, desolvation on specific surfaces can be assisted by anions adsorbed on the crystal surface. This hypothesis, corroborated by crystallization and scanning electron microscopy studies, allows the rationalization of the morphology of barite crystals grown at different supersaturations.

Journal Article↗
PubMed2006

Deuterium isotope effects in A:T and A:U base pairs: a computational NMR study.

Recent measurements of trans-hydrogen bond deuterium isotope effects (DIEs) on 13C chemical shifts in nucleic acids (Vakonakis, I.; LiWang, A. C. J. Biomol. NMR 2004, 29, 65; J. Am. Chem. Soc. 2004, 126, 5688) have led to intriguing results: (i) the DIEs of A:T pairs in DNA are about 5 ppb smaller than those of A:U in RNA and (ii) A:T DIEs vary by as much as 13 ppb among the oligonucleotides. The first observation suggests that inter-base H-bonds in RNA may be stronger than those in DNA, while the second indicates that the conformation of the base pair modulates the transmission of the isotope effect across the hydrogen bond. In an effort at providing a rationale--so far unknown--for the observed DIEs in nucleic acids, density functional theory and hybrid Car-Parrinello/molecular mechanical calculations of DIEs on nucleosides and nucleotides in the gas phase and in aqueous solution have been performed. The calculations suggest that (i) the DIE in an isolated A:T base pair differs from that in an A:U base pair because of the changes in the magnetic properties caused by the replacement of a methyl group on passing from U to T, (ii) the DIEs depend crucially on the conformation of the base pairs, and (iii) the DIEs are strongly affected by magnetic and electrostatic interactions with the surrounding environment.

Journal Article↗
PubMed2005

Structure and energy of a DNA dodecamer under tensile load.

In the last decade, methods to study single DNA molecules under tensile load have been developed. These experiments measure the force required to stretch and melt the double helix and provide insights into the structural stability of DNA. However, it is not easy to directly relate the shape of the force curve to the structural changes that occur in the double helix under tensile load. Here, state-of-the-art computer simulations of short DNA sequences are preformed to provide an atomistic description of the stretching of the DNA double helix. These calculations show that for extensions larger that approximately 25% the DNA undergoes a structural transformation and a few base pairs are lost from both the terminal and central part of the helix. This locally melted DNA duplex is stable and can be extended up to approximately 50-60% of the equilibrium length at a constant force. It is concluded that melting under tension cannot be modeled as a simple two-state process. Finally, the important role of the cantilever stiffness in determining the shape of the force-extension curve and the most probable rupture force is discussed.

Computer Simulation↗
PubMed2005

Simulating micrometre-scale crystal growth from solution.

Understanding crystal growth is essential for controlling the crystallization used in industrial separation and purification processes. Because solids interact through their surfaces, crystal shape can influence both chemical and physical properties. The thermodynamic morphology can readily be predicted, but most particle shapes are actually controlled by the kinetics of the atomic growth processes through which assembly occurs. Here we study the urea-solvent interface at the nanometre scale and report kinetic Monte Carlo simulations of the micrometre-scale three-dimensional growth of urea crystals. These simulations accurately reproduce experimentally observed crystal growth. Unlike previous models of crystal growth, no assumption is made that the morphology can be constructed from the results for independently growing surfaces or from an a priori specification of surface defect concentration. This approach offers insights into the role of the solvent, the degree of supersaturation, and the contribution that extended defects (such as screw dislocations) make to crystal growth. It also connects observations made at the nanometre scale, through in situ atomic force microscopy, with those made at the macroscopic level. If extended to include additives, the technique could lead to the computer-aided design of crystals.

Journal Article↗
PubMed2005

Folding pathways for initiator and effector procaspases from computer simulations.

The folding pathways of procaspases 3, 7, and 8 have been studied using a Go-like Hamiltonian and molecular dynamics simulations coupled with a parallel tempering scheme. The folding pathways and the overall structures of procaspases 3 and 7 are similar, and are characterized by monomeric as well as dimeric folding intermediates in agreement with the available structural and thermochemical data. The folding pathway of procaspase 8, on the other hand, is characterized by a larger population of monomers and partially folded dimer intermediates, and only a relatively small population of folded dimer species. The most stable structure predicted for procaspase 8 is a dimer, in which the position of the linker is remarkably different from the one observed in procaspases 3 and 7, leading to the fact that all the contacts that stabilize the active site are essentially formed. This novel and unexpected structure provides a rationale for the observed activity of the procaspase 8 dimer, and thus could be highly relevant for the initiation of FAS-mediated apoptosis.

Caspase 3↗
PubMed2005

Understanding the barriers to crystal growth: dynamical simulation of the dissolution and growth of urea from aqueous solution.

Both the dissolution and growth of a molecular crystalline material, urea, has been studied using dynamical atomistic simulation. The kinetic steps of dissolution and growth are clearly identified, and the activation energies for each possible step are calculated. Our molecular dynamics simulations indicate that crystal growth on the [001] face is characterized by a nucleation and growth mechanism. Nucleation on the [001] urea crystal face is predicted to occur at a very high rate, followed by rapid propagation of the steps. The rate-limiting step for crystallization is actually found to be the removal of surface defects, rather than the initial formation of the next surface layer. Through kinetic Monte Carlo modeling of the surface growth, it is found that this crystal face evolves via a rough surface topography, rather than a clean layer-by-layer mechanism.

Journal Article↗
PubMed2004

Molecular dynamics simulations of structural changes during procaspase 3 activation.

Molecular dynamics (MD) simulations of the structural rearrangements on the pathway leading to procaspase 3 activation are presented. A retrostructural approach is used to build procaspase 3 from mature caspase 3. The peptide bond that is cleaved during enzyme maturation is gradually reformed during the MD simulation and the most relevant structural changes that occur as a consequence are analyzed. The main structural features that characterize this procaspase 3 model are compared with the available X-ray structure of procaspase 7 as the only zymogen structure that has been crystallised so far. The MD simulations indicate that in the free caspase 3, the flexible selectivity loop is already preorganized to accomodate the substrate. Such a preorganization is not present in either monomeric caspase 3 or in the procaspase 3 dimer, indicating that the structure of the selectivity loop is highly sensitive to perturbations.

Caspase 3↗
PubMed2003

Protonation states of methionine aminopeptidase and their relevance for inhibitor binding and catalytic activity.

We have performed a computational study of different protomeric states of the methionine aminopeptidase active site using a combined quantum-mechanical/molecular mechanical simulation approach. The aim of this study was to clarify the native protonation state of the enzyme, which is needed for the development of novel irreversible inhibitors that can possibly be used as antiangiogenic and antibiotic drugs by virtual screening and other drug design methods. The results of the simulations indicated that two protonation states are possible without disturbing the overall geometry of the active site. We then verified experimentally the presence of the two protonation states by studying the substrate hydrolysis and inhibitor binding reactions at different pH values and come to the conclusion that one of the protomeric states is relevant for inhibitor binding, whereas the other is relevant for substrate hydrolysis. This result has implications for the development of other inhibitors of this class of enzymes and adds a new perspective to the pharmacological properties of the antiangiogenic drug fumagillin, which is an irreversible inhibitor of the human methionine aminopeptidase type II.

Aminopeptidases↗
PubMed2003

Structure-based thermodynamic analysis of caspases reveals key residues for dimerization and activity.

Cysteine-dependent aspartic proteases (caspases) are a family of enzymes which play a crucial role in apoptosis. Caspases accumulate in eukaryotic cells in the form of low-activity proenzyme precursors. Proteolytic cleavage of specific sites triggers conformational changes that lead to full activation and thus to the initiation of the apoptotic cascade. Several experimental observations suggest that dimerization is crucial for activity and regulation, but the underlying molecular mechanisms have not yet been completely resolved. In this work, we have used a structure-based thermodynamic analysis method [Edgcomb, S. P., and Murphy, K. P. (2000) Curr. Opin. Biotechnol. 11, 62-66] to calculate the free energy of association and folding for all the caspases and procaspases whose structures are known at present. In all cases, analysis of the single-residue contributions to the dimerization energy shows that 30-50% of the dimer stability originates from the highly specific interaction of 12-14 residues located at the N- and C-termini of the large and small subunits, respectively. Moreover, our calculations indicate that these residues are also critical for stabilizing the conformation of the active site loops, which in turn is crucial for the binding of substrates and inhibitors. Thus, our results help to rationalize the relation between dimerization and activity in this important class of enzymes and can be used as a starting point for an active manipulation of the monomer-dimer equilibrium for preparatory and regulatory purposes.

Amino Acid Sequence↗
PubMed2002

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease.

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.

Binding Sites↗
PubMed2002

Drug resistance in HIV-1 protease: Flexibility-assisted mechanism of compensatory mutations.

The emergence of drug-resistant variants is a serious side effect associated with acquired immune deficiency syndrome therapies based on inhibition of human immunodeficiency virus type 1 protease (HIV-1 PR). In these variants, compensatory mutations, usually located far from the active site, are able to affect the enzymatic activity via molecular mechanisms that have been related to differences in the conformational flexibility, although the detailed mechanistic aspects have not been clarified so far. Here, we perform multinanosecond molecular dynamics simulations on L63P HIV-1 PR, corresponding to the wild type, and one of its most frequently occurring compensatory mutations, M46I, complexed with the substrate and an enzymatic intermediate. The quality of the calculations is established by comparison with the available nuclear magnetic resonance data. Our calculations indicate that the dynamical fluctuations of the mutated enzyme differ from those in the wild type. These differences in the dynamic properties of the adducts with the substrate and with the gem-diol intermediate might be directly related to variations in the enzymatic activity and therefore offer an explanation of the observed changes in catalytic rate between wild type and mutated enzyme. We anticipate that this "flexibility-assisted" mechanism might be effective in the vast majority of compensatory mutations, which do not change the electrostatic properties of the enzyme.

Drug Resistance, Viral↗