Biomedical subjects

P Carloni

Publications and source records attributed to P Carloni.

At least 19 recordsLinked to original sources

PubMed2006

Using metadynamics to understand the mechanism of calmodulin/target recognition at atomic detail.

The ability of calcium-bound calmodulin (CaM) to recognize most of its target peptides is caused by its binding to two hydrophobic residues ('anchors'). In most of the CaM complexes, the anchors pack against the hydrophobic pockets of the CaM domains and are surrounded by fully conserved Met side chains. Here, by using metadynamics simulations, we investigate quantitatively the energetics of the final step of this process using the M13 peptide, which has a high affinity and spans the sequence of the skeletal myosin light chain kinase, an important natural CaM target. We established the accuracy of our calculations by a comparison between calculated and NMR-derived structural and dynamical properties. Our calculations provide novel insights into the mechanism of protein/peptide recognition: we show that the process is associated with a free energy gain similar to that experimentally measured for the CaM complex with the homologous smooth muscle MLCK peptide (Ehrhardt et al., 1995, Biochemistry 34, 2731). We suggest that binding is dominated by the entropic effect, in agreement with previous proposals. Furthermore, we explain the role of conserved methionines by showing that the large flexibility of these side chains is a key feature of the binding mechanism. Finally, we provide a rationale for the experimental observation that in all CaM complexes the C-terminal domain seems to be hierarchically more important in establishing the interaction.

Amino Acid Sequence↗

PubMed2005

Unwinding the helical linker of calcium-loaded calmodulin: a molecular dynamics study.

The fold of calmodulin (CaM) consists of two globular domains connected by a helical segment (the linker), whose conformational properties play a crucial role for the protein's molecular recognition processes. Here we investigate the structural properties of the linker by performing a 11.5 ns molecular dynamics (MD) simulation of calcium-loaded human CaM in aqueous solution. The calculations are based on the AMBER force field. The calculated S2 order parameters are in good accord with NMR data: The structure of the linker in our simulations is much more flexible than that emerging from the Homo sapiens X-ray structure, consistently with the helix unwinding observed experimentally in solution. This process occurs spontaneously in a nanosecond timescale, as observed also in a very recent simulation based on the GROMOS force field. A detailed description of the mechanism that determines the linker unwinding is provided, in which electrostatic contacts between the two globular domains play a critical role. The orientation of the domains emerging from our MD calculations is consistent both with former X-ray scattering data and a recent NMR work. Based on our findings, a rationale for the experimentally measured entropy cost associated to binding to the protein's cellular partners is also given.

Binding Sites↗

PubMed2005

A homology model of the pore region of HCN channels.

HCN channels are activated by membrane hyperpolarization and regulated by cyclic nucleotides, such as cyclic adenosine-mono-phosphate (cAMP). Here we present structural models of the pore region of these channels obtained by using homology modeling and validated against spatial constraints derived from electrophysiological experiments. For the construction of the models we make two major assumptions, justified by electrophysiological observations: i), in the closed state, the topology of the inner pore of HCN channels is similar to that of K(+) channels. In particular, the orientation of the S5 and S6 helices of HCN channels is very similar to that of the corresponding helices of the K(+) KcsA and K(+) KirBac1.1 channels. Thus, we use as templates the x-ray structure of these K(+) channels. ii), In the open state, the S6 helix is bent further than it is in the closed state, as suggested (but not proven) by experimental data. For this reason, the template of the open conformation is the x-ray structure of the MthK channel. The structural models of the closed state turn out to be consistent with all the available electrophysiological data. The model of the open state turned out to be consistent with all the available electrophysiological data in the filter region, including additional experimental data performed in this work. However, it required the introduction of an appropriate, experimentally derived constraint for the S6 helix. Our modeling provides a structural framework for understanding several functional properties of HCN channels: i), the cysteine ring at the inner mouth of the pore may act as a sensor of the intracellular oxidizing/reducing conditions; ii), the bending amplitude of the S6 helix upon gating appears to be significantly smaller than that found in MthK channels; iii), the reduced ionic selectivity of HCN channels, relative to that of K(+) channels, may be caused, at least in part, by the larger flexibility of the inner pore of HCN channels.

Amino Acid Sequence↗

PubMed2003

Reaction mechanism of caspases: insights from QM/MM Car-Parrinello simulations.

Caspases are fundamental targets for pharmaceutical interventions in a variety of diseases involving disregulated apoptosis. Here, we present a quantum mechanics/molecular mechanics Car-Parrinello study of key steps of the enzymatic reaction for a representative member of this family, caspase-3. The hydrolysis of the acyl-enzyme complex is described at the density functional (BLYP) level of theory while the protein frame and solvent are treated using the GROMOS96 force field. These calculations show that the attack of the hydrolytic water molecule implies an activation free energy of ca. DeltaF(A) approximately equal 19 +/- 4 kcal/mol in good agreement with experimental data and leads to a previously unrecognized gem-diol intermediate that can readily (DeltaF(A) approximately equal 5 +/- 3 kcal/mol) evolve to the enzyme products. Our findings assist in elucidating the striking difference in catalytic activity between caspases and other structurally well-characterized cysteine proteases (papains and cathepsins) and may help design novel transition-state analog inhibitors.

Caspase 3↗

PubMed2003

Molecular dynamics studies of caspase-3.

Caspase-3 is a fundamental target for pharmaceutical interventions against a variety of diseases involving disregulated apoptosis. The enzyme is active as a dimer with two symmetry-related active sites, each featuring a Cys-His catalytic dyad and a selectivity loop, which recognizes the characteristic DEVD pattern of the substrate. Here, a molecular dynamics study of the enzyme in complex with two pentapeptide substrates DEVDG is presented, which provides a characterization of the dynamic properties of the active form in aqueous solution. The mobility of the substrate and that of the catalytic residues are rather low indicating a distinct preorganization effect of the Michaelis complex. An essential mode analysis permits us to identify coupled motions between the two monomers. In particular, it is found that the motions of the two active site loops are correlated and tend to steer the substrate toward the reactive center, suggesting that dimerization has a distinct effect on the dynamic properties of the active site regions. The selectivity loop of one monomer turns out to be correlated with the N-terminal region of the p12 subunit of the other monomer, an interaction that is also found to play a fundamental role in the electrostatic stabilization of the quaternary structure. To further characterize the specific influence of dimerization on the enzyme essential motions, a molecular dynamics analysis is also performed on the isolated monomer.

Binding Sites↗

PubMed2001

Ab initio molecular dynamics-based assignment of the protonation state of pepstatin A/HIV-1 protease cleavage site.

A recent 13C NMR experiment (Smith et al. Nature Struct. Biol. 1996, 3, 946-950) on the Asp 25-Asp25' dyad in pepstatin A/HIV-1 protease measured two separate resonance lines, which were interpreted as being a singly protonated dyad. We address this issue by performing ab initio molecular dynamics calculations on models for this site accompanied by calculations of 13C NMR chemical shifts and isotopic shifts. We find that already on the picosecond time-scale the model proposed by Smith et al. is not stable and evolves toward a different monoprotonated form whose NMR pattern differs from the experimental one. We suggest, instead, a different protonation state in which both aspartic groups are protonated. Despite the symmetric protonation state, the calculated 13C NMR properties are in good agreement with the experiment. We rationalize this result using a simple valence bond model, which explains the chemical inequality of the two C sites. The model calculations, together with our calculations on the complex, allow also the rationalization of 13C NMR properties on other HIV-1 PR/inhibitor complexes. Both putative binding of the substrate to the free enzyme, which has the dyad singly protonated (Piana, S.; Carloni, P. Proteins: Struct., Funct., Genet. 2000, 39, 26-36), and pepstatin A binding to the diprotonated form are consistent with the inverse solvent isotope effect on the onset of inhibition of pepsin by pepstatin and the kinetic iso-mechanism proposed for aspartic proteases (Cho, T.-K.; Rebholz, K.; Northrop, D.B. Biochemistry 1994, 33, 9637-9642).

Aspartic Acid↗

PubMed2001

Molecular dynamics studies on HIV-1 protease drug resistance and folding pathways.

Drug resistance to HIV-1 protease involves the accumulation of multiple mutations in the protein. We investigate the role of these mutations by using molecular dynamics simulations that exploit the influence of the native-state topology in the folding process. Our calculations show that sites contributing to phenotypic resistance of FDA-approved drugs are among the most sensitive positions for the stability of partially folded states and should play a relevant role in the folding process. Furthermore, associations between amino acid sites mutating under drug treatment are shown to be statistically correlated. The striking correlation between clinical data and our calculations suggest a novel approach to the design of drugs tailored to bind regions crucial not only for protein function, but for folding as well.

Amino Acid Sequence↗

PubMed2001

The rational of catalytic activity of herpes simplex virus thymidine kinase. a combined biochemical and quantum chemical study.

Most antiherpes therapies exploit the large substrate acceptance of herpes simplex virus type 1 thymidine kinase (TK(HSV1)) relative to the human isoenzyme. The enzyme selectively phosphorylates nucleoside analogs that can either inhibit viral DNA polymerase or cause toxic effects when incorporated into viral DNA. To relate structural properties of TK(HSV1) ligands to their chemical reactivity we have carried out ab initio quantum chemistry calculations within the density functional theory framework in combination with biochemical studies. Calculations have focused on a set of ligands carrying a representative set of the large spectrum of sugar-mimicking moieties and for which structural information of the TK(HSV1)-ligand complex is available. The k(cat) values of these ligands have been measured under the same experimental conditions using an UV spectrophotometric assay. The calculations point to the crucial role of electric dipole moment of ligands and its interaction with the negatively charged residue Glu(225). A striking correlation is found between the energetics associated with this interaction and the k(cat) values measured under homogeneous conditions. This finding uncovers a fundamental aspect of the mechanism governing substrate diversity and catalytic turnover and thus represents a significant step toward the rational design of novel and powerful prodrugs for antiviral and TK(HSV1)-linked suicide gene therapies.

Antiviral Agents↗

PubMed2001

Ser133 phosphate-KIX interactions in the CREB-CBP complex: an ab initio molecular dynamics study.

Cyclic AMP response element binding protein (CREB) is involved in activation of transcriptional DNA machinery by binding to the coactivator CREB-binding protein (CBP). The interactions between CREB serine phosphate (pSer133) and specific CBP residues (Tyr658 and Lys662) play a crucial role for the thermodynamic stability of the CREB-CBP complex. Here we use ab initio methods to investigate the dynamics and energetics of a relatively large, fully hydrated model complex representing pSer133 and its counterparts of the CBP domain. The calculations suggest that: (1) key contributions to the stabilization of the complex arise not only from electrostatics (as previously proposed) but also from a previously unrecognized "low-barrier hydrogen bond" between pSer133 and Lys662; (2) hydration plays a crucial role for the stabilization of the phosphate charge; (3) formation of the complex involves a significant degree of reorganization of the electronic charge density.

Binding Sites↗

PubMed2000

Water and potassium dynamics inside the KcsA K(+) channel.

Molecular dynamics simulations and electrostatic modeling are used to investigate structural and dynamical properties of the potassium ions and of water molecules inside the KcsA channel immersed in a membrane-mimetic environment. Two potassium ions, initially located in the selectivity filter binding sites, maintain their position during 2 ns of dynamics. A third potassium ion is very mobile in the water-filled cavity. The protein appears engineered so as to polarize water molecules inside the channel cavity. The resulting water induced dipole and the positively charged potassium ion within the cavity are the key ingredients for stabilizing the two K(+) ions in the binding sites. These two ions experience single file movements upon removal of the potassium in the cavity, confirming the role of the latter in ion transport through the channel.

Bacterial Proteins↗

PubMed2000

Conformational flexibility of the catalytic Asp dyad in HIV-1 protease: An ab initio study on the free enzyme.

The enzyme protease from the human immunodeficiency virus type 1 (HIV-1 PR) is one of the main targets for therapeutic intervention in AIDS. Computer modeling is useful for probing the binding of novel ligands, yet empirical force field-based methods have encountered problems in adequately describing interactions of the catalytic aspartyl pair. In this work we use ab initio dynamic methods to study the molecular interactions and the conformational flexibility of the Asp dyad in the free enzyme. Calculations are performed on model complexes that include, besides the Asp dyad, the conserved Thr26 and Gly27 residues and water molecules present in the active site channel. Our calculations provide proton location and binding mode of the active-site water molecule, which turn out to be different from those of the eukariotic isoenzyme. Furthermore, the calculations reproduce well the structural features of the aspartyl dyad in the protein. Finally, they allow the identification of both dipole/charge interactions and a low-barrier hydrogen bond as important stabilizing factors for the peculiar conformation of the active site. These findings are consistent with site-directed mutagenesis experiments on the 27, 27; positions (Bagossi et al., Protein Eng 1996;9:997-1003). The electric field of the protein frame (included in some of the calculations) does not affect significantly the chemical bonding at the cleavage site. Proteins 2000;39:26-36.

Aspartic Acid↗

PubMed2000

Evaluation of intramolecular equilibria in complexes formed between substituted imidazole ligands and nickel (II), copper (II) or zinc (II).

The metal ion-binding properties of imidazole-4-acetate (ImA-), 4(5)-aminoimidazole-5(4)-carboxamide (AImC), 2,2'biimidazole(BiIm) (I. Török et al., J. Inorg. Biochem. 71 (1998) 7-14), and bis (imidazol-2-yl)methane(BiImM) (K. Várnagy et al., J. Chem. Soc., Dalton Trans. (1994) 2939-2945) have been evaluated by using the recently published stability constants and by applying the recently established log K(ML)M versus pK(HL)H straight-line plots (L. E. Kapinos et al., Inorg. Chim. Acta 280 (1998) 50-56) which hold for simple imidazole-type ligands. The indicated analysis regarding the intramolecular equilibrium between a monodentatally imidazole-nitrogen-coordinated (open) species and a chelated isomer provides helpful insights, e.g., the formation degree of chelates is more favored if six-membered rings can be formed, as in the case with M(BiImM)2+ compared to M(BiIm)2+, though in both instances the formation degree of the chelates is large. The formation degree of chelates in the M(ImA)+ complexes increases in the series Zn(ImA)+ (87%)<Ni(ImA)+ (96%)<Cu(ImA)+ (99.5%). A carbonyl oxygen, if sterically favorably positioned as in the M (AImC)2+ complexes, may also participate well in chelate formation. In this way a carbonyl group can certainly be activated via metal ion coordination and become ready for further reactions. For Ni2+, Cu2+, and Zn2+ the formation degree of the chelated M(AImC)2+ isomers varies between about 30 and 75%. In all of the so-called 'open' species the metal ion is solely coordinated to a pyridine-like nitrogen of the imidazole residue. Some further observations of biological interest are indicated.

Cations, Divalent↗

PubMed2000

Increased oxidative modification of albumin when illuminated in vitro in the presence of a common sunscreen ingredient: protection by nitroxide radicals.

We previously reported on the ability of dibenzoylmethane (DBM) and a relative, Parsol 1789, used as a ultraviolet A (UVA)-absorbing sunscreen, to generate free radicals upon illumination, and as a consequence, to inflict strand breaks in plasmid DNA in vitro. This study has now been extended to determine the effects of Parsol 1789 and DBM on proteins, under UVA illumination, with the sole purpose of gaining more knowledge on the photobiological effects of sunscreen chemicals. Parsol 1789 (100 microM) caused a 2-fold increase in protein carbonyl formation (an index of oxidative damage) in bovine serum albumin (BSA) when exposed to illumination, and this damage was both concentration- and time-dependent. The degree of protein damage was markedly reduced by the presence of free radical scavengers, namely piperidinic and indolinonic nitroxide radicals, in accordance with our previous study. Vitamin E had no effect under the conditions used. The results obtained corroborate the fact that Parsol 1789 generates free radicals upon illumination and that these are, most probably, responsible for the protein damage observed under the conditions used in our system. However, at present, we cannot extrapolate from these results the relevance to human use of sunscreens; therefore, further studies should be necessary to determine the efficacy at the molecular and cellular level of this UVA-absorber in order to ascertain protection against photocarcinogenic risk.

Animals↗

PubMed2000

A comparative study of galactose oxidase and active site analogs based on QM/MM Car-Parrinello simulations.

A parallel study of the radical copper enzyme galactose oxidase (GOase) and a low molecular weight analog of the active site was performed with dynamical density functional and mixed quantum-classical calculations. This combined approach enables a direct comparison of the properties of the biomimetic and the natural systems throughout the course of the catalytic reaction. In both cases, five essential forms of the catalytic cycle have been investigated: the resting state in its semi-reduced (catalytically inactive) and its oxidized (catalytically active) form, A(semi) and A(ox), respectively; a protonated intermediate B; the transition state for the rate-determining hydrogen abstraction step C, and its product D. For A and B the electronic properties of the biomimetic compound are qualitatively very similar to the ones of the natural target. However, in agreement with the experimentally observed difference in catalytic activity, the calculated activation energy for the hydrogen abstraction step is distinctly lower for GOase (16 kcal/mol) than for the mimetic compound (21 kcal/mol). The enzymatic transition state is stabilized by a delocalization of the unpaired spin density over the sulfur-modified equatorial tyrosine Tyr272, an effect that for geometric reasons is essentially absent in the biomimetic compound. Further differences between the mimic and its natural target concern the structure of the product of the abstraction step, which is characterized by a weakly coordinated aldehyde complex for the latter and a tightly bound linear complex for the former.

Binding Sites↗

PubMed2000

Ab initio molecular dynamics studies on HIV-1 reverse transcriptase triphosphate binding site: implications for nucleoside-analog drug resistance.

Quantum-chemical methods are used to shed light on the functional role of residues involved in the resistance of HIV-1 reverse transcriptase against nucleoside-analog drugs. Ab initio molecular dynamics simulations are carried out for models representing the adduct between the triphosphate substrate and the nucleoside binding site. The triphosphate is considered either deprotonated or protonated at the gamma-position. Although the protonated form already experiences large rearrangements in the ps time scale, the fully deprotonated state exhibits a previously unrecognized low-barrier hydrogen bond between Lys65 and gamma-phosphate. Absence of this interaction in Lys65-->Arg HIV-1 RT might play a prominent role in the resistance of this mutant for nucleoside analogs (Gu Z et al., 1994b, Antimicrob Agents Chemother 38:275-281; Zhang D et al., 1994, Antimicrob Agents Chemother 38:282-287). Water molecules present in the active site, not detected in the X-ray structure, form a complex H-bond network. Among these waters, one may be crucial for substrate recognition as it bridges Gln151 and Arg72 with the beta-phosphate. Absence of this stabilizing interaction in Gln151-->Met HIV-1 RT mutant may be a key factor for the known drug resistance of this mutant toward dideoxy-type drugs and AZT (Shirasaka T et al., 1995, Proc Natl Acad Sci USA 92:2398-2402: Iversen AK et al., 1996, J Virol 70:1086-1090).

Amino Acid Substitution↗

PubMed1999

Serine proteases: an ab initio molecular dynamics study.

In serine proteases (SPs), the H-bond between His57 and Asp102 and that between Gly193 and the transition state intermediate play a crucial role in enzymatic function. To shed light on the nature of these interactions, we have carried out ab initio molecular dynamics simulations on complexes representing adducts between the reaction intermediate and elastase (one protein belonging to the SP family). Our calculations indicate the presence of a low-barrier H-bond between His57 and Asp102, in complete agreement with NMR experiments on enzyme-transition state analogue complexes. Comparison with an ab initio molecular dynamics simulation on a model of the substrate-enzyme adduct indicates that the Gly193-induced strong stabilization of the intermediate is accomplished by charge/dipole interactions and not by H-bonding as previously suggested. Inclusion of the protein electric field in the calculations does not affect significantly the charge distribution.

Enzyme Stability↗

PubMed1999

Protein design is a key factor for subunit-subunit association.

Fundamental questions about role of the quaternary structures are addressed by using a statistical mechanics off-lattice model of a dimer protein. The model, in spite of its simplicity, captures key features of the monomer-monomer interactions revealed by atomic force experiments. Force curves during association and dissociation processes are characterized by sudden jumps followed by smooth behavior and form hysteresis loops. Furthermore, the process is reversible in a finite range of temperature stabilizing the dimer, and the width of the hysteresis loop increases as the design procedure improves: i.e., stabilizes the dimer more. It is shown that, in the interface between the two monomeric subunits, the design procedure naturally favors those amino acids whose mutual interaction is stronger.

Dimerization↗

PubMed1999

Potassium and sodium binding to the outer mouth of the K+ channel.

Molecular dynamics simulations of the K+ channel from Streptomyces lividans (KcsA channel) were performed in a membrane-mimetic environment with Na+ and K+ in different initial locations. The structure of the channel remained stable and well preserved for simulations lasting up to 1.5 ns. Salt bridges between Asp80 and Arg89 of neighboring subunits, not detected in the X-ray structure, enhanced the stability of the tetrameric structure. Na+ or K+ ions located in the channel vestibule lost part of their hydration shell and diffused into the channel inner pore in less than a few hundred picoseconds. This powerful catalytic action was caused by strong electrostatic interactions with Asp80 and Glu71. The hydration state of the metal ions turned out to depend significantly on the conformational flexibility of the channel. Furthermore, Na+ entered the channel inner pore bound to more water molecules than K+. The different hydration state of the two ions may be a determinant factor in the ion selectivity of the channel.

Binding Sites↗