Biomedical subjects

António M Baptista

Publications and source records attributed to António M Baptista.

18 recordsLinked to original sources

PubMed2006

The pH-dependent conformational states of kyotorphin: a constant-pH molecular dynamics study.

An extensive conformational study of the analgesic dipeptide kyotorphin (L-Tyr-L-Arg) at different pH values was performed using a constant-pH molecular dynamics method. This dipeptide showed a remarkable pH-dependent conformational variety. The protonation of the N-terminal amine was identified as a key element in the transition between the more extended and the more packed conformational states, as monitored by the dihedral angle defined by the atoms 1Cbeta-1Calpha-2Calpha-2Cbeta. The principal-component analysis of kyotorphin identified two major conformational populations (the extended trans and the packed cis) together with conformations that occur exclusively at extreme pH values. Other, less stable conformations were also identified, which help us to understand the transitions between the predominant populations. The fitting of kyotorphin's conformational space to the structure of morphine resulted in a set of conformers that were able to fulfill most of the constraints for the mu-receptor. These results suggest that there may be strong similarities between the kyotorphin receptor and the structural family of opioid receptors.

Computer Simulation↗

PubMed2006

Natural domain design: enhanced thermal stability of a zinc-lacking ferredoxin isoform shows that a hydrophobic core efficiently replaces the structural metal site.

Zinc centers play a key role as important structure determinants in a variety of proteins including ferredoxins (Fd). Here, we exploit the availability of two highly similar ferredoxin isoforms from the thermophile Sulfolobus metallicus, which differ in the residues involved in coordinating a His/Asp zinc site that ties together the protein core with its N-terminal extension, to investigate the effect of the absence of this site on ferredoxin folding. The conformational properties of the zinc-containing (FdA) and zinc-lacking (FdB) isoforms were investigated using visible absorption and tryptophan fluorescence emission. Fluorescence quenching studies, together with comparative modeling and molecular dynamics simulations, indicate that the FdB N-terminal extension assumes a fold identical to that of the Zn(2+)-containing isoform. The thermal stability of the isoforms was investigated in a broad pH range (2 < pH < 10), and at physiological pH conditions, both proteins unfold above 100 degrees C. Surprisingly, the Zn(2+)-lacking isoform was always found to be more stable than its Zn(2+)-containing counterpart: a DeltaT(m) approximately 9 degrees C is determined at pH 7, a difference that becomes even more significant at extreme pH values, reaching a DeltaT(m) approximately 24 degrees C at pH 2 and 10. The contribution of the Zn(2+) site to ferredoxin stability was further resolved using selective metal chelators. During thermal unfolding, the zinc scavenger TPEN significantly lowers the T(m) in FdA ( approximately 10 degrees C), whereas it has no effect in FdB. This shows that the Zn(2+) site contributes to ferredoxin stability but that FdB has devised a structural strategy that accounts for an enhanced stability without using a metal cross-linker. An analysis of the FdB sequence and structural model leads us to propose that the higher stability of the zinc-containing ferredoxin results from van der Waals contacts formed between the residues that occupy the same spatial region where the zinc ligands are found in FdA. These favor the formation of a novel local stabilizing hydrophobic core and illustrate a strategy of natural fold design.

Archaeal Proteins↗

PubMed2006

Parametrization of 1-butyl-3-methylimidazolium hexafluorophosphate/nitrate ionic liquid for the GROMOS force field.

A united-atom model of 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF(6)]) and 1-butyl-3-methylimidazolium nitrate ([BMIM][NO(3)]) is developed in the framework of the GROMOS96 43A1(1) force field. These two ionic liquids are parametrized, and their equilibrium properties in the 298-363 K temperature range are subjected to validation against known experimental properties, namely, density, self-diffusion, shear viscosity, and isothermal compressibility. The ionic radial/spatial distributions, pi interaction, gauche/trans populations of the butyl tail, and enthalpies of vaporization are also reported. The properties obtained from the molecular dynamics simulations agree with experimental data and have the same temperature dependence. The strengths and weakness of our model are discussed.

Computer Simulation↗

PubMed2006

Pathways of H2 toward the active site of [NiFe]-hydrogenase.

Hydrogenases catalyze the reversible oxidation of molecular hydrogen (H(2)), but little is known about the diffusion of H(2) toward the active site. Here we analyze pathways for H(2) permeation using molecular dynamics (MD) simulations in explicit solvent. Various MD simulation replicates were done, to improve the sampling of the system states. H(2) easily permeates hydrogenase in every simulation and it moves preferentially in channels. All H(2) molecules that reach the active site made their approach from the side of the Ni ion. H(2) is able to reach distances of <4 A from the active site, although after 6 A permeation is difficult. In this region we mutated Val-67 into alanine and perform new MD simulations. These simulations show an increase of H(2) inside the protein and at lower distances from the active site. This valine can be a control point in the H(2) access to the active center.

Binding Sites↗

PubMed2006

Conformational and orientational guidance of the analgesic dipeptide kyotorphin induced by lipidic membranes: putative correlation toward receptor docking.

The analgesic dipeptide kyotorphin (L-Tyr-L-Arg) and an acylated kyotorphin derivative were studied by a combination of theoretical (molecular dynamics simulation and quantum mechanics methods) and experimental (fluorescence and infrared spectroscopies) approaches both in solution and in model systems of membranes. At biological pH the peptides have a neutral net charge. Nevertheless, their phenolic rings interact with phospholipid molecules (partition coefficient varies from 6 x 10(2) to 2 x 10(4), depending on the lipid and pH used) despite being exposed to the aqueous bulk medium. The lowest energy transition dipole moment is displaced from the normal to the lipid bilayer by 20 degrees on average. The observed extensive interaction, pK(a), precise location, and well-defined orientation in membranes combined with the ability to discriminate rigid raftlike membrane domains suggest that kyotorphin meets the structural constraints needed for receptor-ligand interaction. The acylated kyotorphin derivative mimics kyotorphin properties and represents a promising way for entrapment in a drug carrier and transport across the blood-brain barrier.

Blood-Brain Barrier↗

PubMed2006

Constant-pH molecular dynamics with ionic strength effects: protonation-conformation coupling in decalysine.

A new implementation of the stochastic titration method for constant-pH molecular dynamics is presented, which introduces ionic strength effects in the simulations. In addition, the new implementation uses a faster molecular dynamics algorithm and an improved treatment of protonation events and of their effect on force field parameters. This new methodology is applied to a decalysine peptide, yielding very good quantitative agreement with experiments, both in terms of titration and helix-coil transition. The results show a significant dependence on ionic strength, illustrating the importance of including this parameter in constant-pH molecular dynamics simulations. Overall, the method seems to properly capture the protonation-conformation coupling and its dependence on ionic strength.

Algorithms↗

PubMed2005

Reorganization and conformational changes in the reduction of tetraheme cytochromes.

Molecular dynamics simulation (MD) constitutes an alternative to time-consuming experiments for studying conformational changes. We apply MD on a redox system where experimental information exists for the fully oxidized and fully reduced states: tetraheme cytochrome c3. Instead of doing one simulation for each state, we apply 10 4-ns replicas for both states, which provides robust statistics to characterize the redox changes. Besides these long simulations, we perform 120 short ones (50 ps), where an equilibrated oxidized state is perturbed to a reduced state. This allows the application of a nonequilibrium method, the subtraction technique, which makes it possible to characterize the different timescales of conformational changes. Reduction induces conformational changes in the N-terminus and on the loops spanning residues 36-42 and 88-93, which correlate very well with experiments, demonstrating the applicability of this methodology. We also analyze the effect of reduction on hydrogen bonds, solvent accessible surface and bound water, the changes being found to involve the hemes and propionate groups. Redox-induced protonation is also investigated, by protonating the propionates D from hemes I and IV. Although this change in the former does not have major conformational consequences, it induces in the latter conformational changes beyond the ones obtained with reduction.

Computer Simulation↗

PubMed2005

On the use of different dielectric constants for computing individual and pairwise terms in poisson-boltzmann studies of protein ionization equilibrium.

Poisson-Boltzmann (PB) models are a fast and common tool for studying electrostatic processes in proteins, particularly their ionization equilibrium (protonation and/or reduction), often yielding quite good results when compared with more detailed models. Yet, they are conceptually very simple and necessarily approximate, their empirical character being most evident when it comes to the choice of the dielectric constant assigned to the protein region. The present study analyzes several factors affecting the ability of PB-based methods to model protein ionization equilibrium. We give particular attention to a suggestion made by Warshel and co-workers (e.g., Sham et al. J. Phys. Chem. B 1997, 101, 4458) of using different protein dielectric constants for computing the individual (site) and the pairwise (site-site) terms of the ionization free energies. Our prediction of pK(a) values for several proteins indicates that no advantage is obtained by such a procedure, even for sites that are buried and/or display large pK(a) shifts relative to the solution values. In particular, the present methodology gives the best predictions using a dielectric constant around 20, for shifted/buried and nonshifted/exposed sites alike. The similarities and differences between the PB model and Warshel's PDLD/S model are discussed, as well as the reasons behind their apparently discrepant results. The present PB model is shown to predict also good reduction potentials in redox proteins.

Animals↗

PubMed2005

Water dependent properties of cutinase in nonaqueous solvents: a computational study of enantioselectivity.

The catalytic properties of enzymes in nonaqueous solvents are known to be dependent on the nature of the solvent. Here we present a molecular modeling study of the enantioselective properties of the enzyme cutinase in hexane under varying hydration conditions. Previous simulation studies have shown that for this model enzyme in hexane, the structural and dynamical properties are affected by the amount of water associated with the protein, being more similar to the aqueous simulation at 5-10% of water content. The implications of the hydration levels on the enzyme resolution of (R,S)-1-phenylethanol and (R,S)-2-phenyl-1-propanol are investigated using free energy calculations of the tetrahedral intermediate (TI) model. With this model system we show that the enzyme enantioselective properties are under the control of the amount of water present in the organic media. Maximum enantioselectivity is achieved at 10% water content. The stabilizing effects of the catalytic histidine on the TI are evaluated at different water contents and shown to be correlated. The correlation between the amount of water present in the media and the structural, dynamical, and thermodynamic properties of the enzyme are examined as well as the active site discriminative power.

1-Propanol↗

PubMed2005

Electric-field-induced redox potential shifts of tetraheme cytochromes c3 immobilized on self-assembled monolayers: surface-enhanced resonance Raman spectroscopy and simulation studies.

The tetraheme protein cytochrome c(3) (Cyt-c(3)) from Desulfovibrio gigas, immobilized on a self-assembled monolayer (SAM) of 11-mercaptoundecanoic acid, is studied by theoretical and spectroscopic methods. Molecular dynamics simulations indicate that the protein docks to the negatively charged SAM via its lysine-rich domain around the exposed heme IV. Complex formation is associated with only little protein structural perturbations. This finding is in line with the resonance Raman and surface-enhanced resonance Raman (SERR) spectroscopic results that indicate essentially the same heme pocket structures for the protein in solution and adsorbed on SAM-coated Ag electrodes. Electron- and proton-binding equilibrium calculations reveal substantial negative shifts of the redox potentials compared to the protein in solution. The magnitude of these shifts decreases in the order heme IV (-161 mV) > heme III (-73 mV) > heme II (-57 mV) > heme I (-26 mV), resulting in a change of the order of reduction. These shifts originate from the distance-dependent electrostatic interactions between the SAM headgroups and the individual hemes, leading to a stabilization of the oxidized forms. The results of the potential-dependent SERR spectroscopic analyses are consistent with the theoretical predictions and afford redox potential shifts of -160 mV (heme IV), -90 mV (heme III), -70 mV (heme II), and +20 mV (heme I) relative to the experimental redox potentials for Cyt-c(3) in solution. SERR spectroscopic experiments reveal electric-field-induced changes of the redox potentials also for the structurally very similar Cyt-c(3) from Desulfovibrio vulgaris, although the shifts are somewhat smaller compared to Cyt-c(3) from D. gigas. This study suggests that electric-field-induced redox potential shifts may also occur upon binding to biomembranes or partner proteins and thus may affect biological electron transfer processes.

Biophysical Phenomena↗

PubMed2004

Theoretical identification of proton channels in the quinol oxidase aa3 from Acidianus ambivalens.

Heme-copper oxidases are membrane proteins found in the respiratory chain of aerobic organisms. They are the terminal electron acceptors coupling the translocation of protons across the membrane with the reduction of oxygen to water. Because the catalytic process occurs in the heme cofactors positioned well inside the protein matrix, proton channels must exist. However, due to the high structural divergence among this kind of proteins, the proton channels previously described are not necessarily conserved. In this work we modeled the structure of the quinol oxidase from Acidianus ambivalens using comparative modeling techniques for identifying proton channels. Additionally, given the high importance that water molecules may have in this process, we have developed a methodology, within the context of comparative modeling, to identify high water probability zones and to deconvolute them into chains of ordered water molecules. From our results, and from the existent information from other proteins from the same superfamily, we were able to suggest three possible proton channels: one K-, one D-, and one Q-spatial homologous proton channels. This methodology can be applied to other systems where water molecules are important for their biological function.

Acidianus↗

PubMed2004

Reorganization in apo- and holo-beta-lactoglobulin upon protonation of Glu89: molecular dynamics and pKa calculations.

Molecular dynamics (MD) simulations starting from crystallographic data allowed us to directly account for the effects of the protonation state of Glu89 on the conformational stability of apo- and holo-beta-lactoglobulin (BLG). In apo-BLG simulations starting from the protonated crystal structure, we observe a long-lived H-bond interaction between the protonated Glu89 and Ser116. This interaction, sequestering the proton from the aqueous medium, explains a pK(half) value evaluated at pH 7.3 by continuum electrostatics/Monte Carlo computation on MD data, using linear response approximation. A very large root-mean-square deviation (RMSD; 5.11 A) is observed for the EF loop between protonated and unprotonated apo-BLG. This results from a quite different orientation of the EF loop that acts either as a closed or as an open lid above the protein calyx. Proton exchange by Glu89 in apo- but not in holo-BLG is associated with a reorganization energy of 4.7 kcal/mol. A 3-ns MD simulation starting from the crystal structure of protonated apo-BLG, but considering the Glu89 as unprotonated, shows the progressive opening of the lid giving rise to the Tanford transition. In both holo-BLG forms, the lid is most probably held in place by hydrophobic interactions of amino acid side-chains of the EF loop with the palmitate hydrocarbon tail.

Animals↗

PubMed2004

Molecular basis for redox-Bohr and cooperative effects in cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies of oxidized and reduced high-resolution structures at pH 7.6.

The tetraheme cytochrome c3 is a small metalloprotein with ca. 13,000 Da found in sulfate-reducing bacteria, which is believed to act as a partner of hydrogenase. The three-dimensional structure of the oxidized and reduced forms of cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774 at pH 7.6 were determined using high-resolution X-ray crystallography and were compared with the previously determined oxidized form at pH 4.0. Theoretical calculations were performed with both structures, using continuum electrostatic calculations and Monte Carlo sampling of protonation and redox states, in order to understand the molecular basis of the redox-Bohr and cooperativity effects related to the coupled transfer of electrons and protons. We were able to identify groups that showed redox-linked conformational changes. In particular, Glu61, His76, and propionate D of heme II showed important contributions to the redox-cooperativity, whereas His76, propionate A of heme I, and propionate D of heme IV were the key residues for the redox-Bohr effect. Upon reduction, an important movement of the backbone region surrounding hemes I and II was also identified, that, together with a few redox-linked conformational changes in side-chain residues, results in a significant decrease in the solvent accessibility of hemes I and II.

Amino Acids↗

PubMed2004

Modeling electron transfer thermodynamics in protein complexes: interaction between two cytochromes c(3).

Redox protein complexes between type I and type II tetraheme cytochromes c(3) from Desulfovibrio vulgaris Hildenborough are here analyzed using theoretical methodologies. Various complexes were generated using rigid-body docking techniques, and the two lowest energy complexes (1 and 2) were relaxed using molecular dynamics simulations with explicit solvent and subjected to further characterization. Complex 1 corresponds to an interaction between hemes I from both cytochromes c(3). Complex 2 corresponds to an interaction between the heme IV from type I and the heme I from type II cytochrome c(3). Binding free energy calculations using molecular mechanics, Poisson-Boltzmann, and surface accessibility methods show that complex 2 is more stable than complex 1. Thermodynamic calculations on complex 2 show that complex formation induces changes in the reduction potential of both cytochromes c(3), but the changes are larger in the type I cytochrome c(3) (the largest one occurring on heme IV, of approximately 80 mV). These changes are sufficient to invert the global titration curves of both cytochromes, generating directionally in electron transfer from type I to type II cytochrome c(3), a phenomenon of obvious thermodynamic origin and consequences, but also with kinetic implications. The existence of processes like this occurring at complex formation may constitute a natural design of efficient redox chains.

Biophysics↗

PubMed2003

Investigation of protonatable residues in Rhodothermus marinus caa3 haem-copper oxygen reductase: comparison with Paracoccus denitrificans aa3 haem-copper oxygen reductase.

The Rhodothermus marinus caa(3 )haem-copper oxygen reductase contains all the residues of the so-called D- and K-proton channels, with the notable exception of the helix VI glutamate residue (Glu278(I) in Paracoccus denitrificans aa(3)), being nevertheless a true oxygen reductase reducing O(2) to water, and an efficient proton pump. Instead, in the same helix, but one turn below, it has a tyrosine residue (Tyr256(I), R. marinus caa(3) numbering), whose hydroxyl group occupies the same spatial position as the carboxylate group of Glu278(I), as deduced by comparative modelling techniques. Therefore, we proposed previously that this tyrosine residue could play an important role in the proton pathway. In this article we further study this hypothesis, by investigating the equilibrium thermodynamics of protonation in R. marinus caa(3), using theoretical methodologies based on the structural model previously obtained. Control calculations are also performed for the P. denitrificans aa(3) oxygen reductase. In both oxygen reductases we find several residues that are proton active (i.e., that display partial protonation) at physiological pH, some of them being redox sensitive (i.e., sensitive to the protein redox state). However, the caa(3 )Tyr256(I) is not proton active at physiological pH, in contrast to the aa(3) Glu278(I) which is both proton active at physiological pH and shows a high redox sensitivity. In R. marinus caa(3) we do not find any other residues in the same protein zone that can have this property. Therefore, there are no putative D-channel residues that are proton active in this oxidase. The protonatable residues of the K-channel are much more functionally conserved in both oxygen reductases than the same type of residues in the D-channel. Two (Tyr262(I) and Lys336(I), caa(3) numbering) out of three protonatable K-channel residues are proton active and redox sensitive in both proteins.

Amino Acids↗

PubMed2003

Redox-Bohr and other cooperativity effects in the nine-heme cytochrome C from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies.

The nine-heme cytochrome c is a monomeric multiheme cytochrome found in Desulfovibrio desulfuricans ATCC 27774. The polypeptide chain comprises 296 residues and wraps around nine hemes of type c. It is believed to take part in the periplasmic assembly of proteins involved in the mechanism of hydrogen cycling, receiving electrons from the tetraheme cytochrome c3. With the purpose of understanding the molecular basis of electron transfer processes in this cytochrome, we have determined the crystal structures of its oxidized and reduced forms at pH 7.5 and performed theoretical calculations of the binding equilibrium of protons and electrons in these structures. This integrated study allowed us to observe that the reduction process induced relevant conformational changes in several residues, as well as protonation changes in some protonatable residues. In particular, the surroundings of hemes I and IV constitute two areas of special interest. In addition, we were able to ascertain the groups involved in the redox-Bohr effect present in this cytochrome and the conformational changes that may underlie the redox-cooperativity effects on different hemes. Furthermore, the thermodynamic simulations provide evidence that the N- and C-terminal domains function in an independent manner, with the hemes belonging to the N-terminal domain showing, in general, a lower redox potential than those found in the C-terminal domain. In this way, electrons captured by the N-terminal domain could easily flow to the C-terminal domain, allowing the former to capture more electrons. A notable exception is heme IX, which has low redox potential and could serve as the exit path for electrons toward other proteins in the electron transfer pathway.

Binding Sites↗

PubMed2003

Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies.

Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.

Carboxylic Ester Hydrolases↗

PubMed2001

Studies of the reduction and protonation behavior of tetraheme cytochromes using atomic detail.

A comparative study of tetraheme cytochrome c3 molecules from several species was carried out using recently developed theoretical methods based on continuum electrostatics. The binding joint equilibrium of electrons and protons was simulated, revealing the complete thermodynamic aspects of electron-proton coupling in these molecules. The method yields excellent accuracy in terms of midpoint potentials, giving the correct reduction orders in all molecules examined, except for one heme site. The coupling between electrons and protons is shown to be present and significant at physiological pH in all cases. This phenomenon, known as the redox-Bohr effect, though of thermodynamic nature, is shown to have an intrinsic "dynamic" character at the molecular level (in the sense of the empty/occupied fluctuations at the microscopic level), with the binding states of redox and protonatable sites displaying both correlated averages and correlated fluctuations. The protonatable sites more directly involved in the redox-Bohr effect are identified using, among other properties, the statistical correlation between pairs of sites, which automatically reflects indirect effects mediated by other sites. Several sites are identified in this analysis. Propionate D of heme I seems to be the most interesting, generally showing a high correlation not only with its own heme, but also with heme II, corresponding to an indirect stabilization of the reduced forms of both hemes. Other interesting sites are the free histidines of two of the cytochromes and propionate D of heme IV, the latter being potentially associated with redox-induced structural changes. Among the set of cytochromes c3 analyzed in this study, significant differences are observed for several properties of the acidic cytochrome included in the set, from Desulfovibrio africanus, supporting the hypothesis of a different functional role.

Computer Simulation↗