Biomedical subjects

Victor Guallar

Publications and source records attributed to Victor Guallar.

14 recordsLinked to original sources

PubMed2006

Hydrogen bonding pathways in human dihydroorotate dehydrogenase.

Dihydroorotate dehydrogenase (DHOD) catalyzes the only redox reaction in the pathway for pyrimidine biosynthesis. In this reaction, a proton is transferred from a carbon atom of the substrate to a serine residue, and a hydride is transferred from another carbon atom of the substrate to a cofactor. The deprotonation of the substrate is postulated to involve a proton relay mechanism along a hydrogen bonding pathway in the active site. In this paper, molecular dynamics simulations are used to identify and characterize potential hydrogen bonding pathways that could facilitate the redox reaction catalyzed by human DHOD. The observed pathways involve hydrogen bonding of the active base serine to a water molecule, which is hydrogen bonded to the substrate carboxylate group or a threonine residue. The threonine residue is positioned to enable proton transfer to another water molecule leading to the bulk solvent. The impact of mutating the active base serine to cysteine is also investigated. This mutation is found to increase the average donor-acceptor distances for proton and hydride transfer and to disrupt the hydrogen bonding pathways observed for the wild-type enzyme. These effects could lead to a significant decrease in enzyme activity, as observed experimentally for the analogous mutant in Escherichia coli DHOD.

Carbon↗

PubMed2006

Modeling of ligation-induced helix/loop displacements in myoglobin: toward an understanding of hemoglobin allostery.

Combining quantum and molecular mechanics (QM/MM) methods and protein structure prediction algorithms, helix and loop movements are computed along the pathway of CO dissociation from myoglobin (Mb). The results are compared with high-resolution crystallographic data using sequence-displacement graphs. These graphs provide an unbiased method for evaluating main-chain segmental motions; they resolve an apparent disagreement between two sets of high-resolution crystal structures for MbCO and deoxyMb. The QM/MM modeling of the CO deligation reproduces the experimentally observed spin states and photodissociated crystal structure. The principal effect of CO dissociation is shown to be a concerted rotation of the E and F helices, which hold the heme like a clamshell. The rotation is a response to deligation forces, which impel the F helix away from the heme because of the Fe spin conversion, and which allow the E helix to collapse toward the heme as nonbonded contacts on the distal side are relieved. Additional helix and loop displacements stem from these primary events. In particular, the CD loop is found to be repositioned as a result of steric interactions with the water molecule that becomes H-bonded to the distal histidine in deoxyMb. A similar EF rotation and CD loop displacement are proposed to be the first steps along the allosteric pathway from the R to the T state in hemoglobin.

Allosteric Regulation↗

PubMed2006

The effect of heme environment on the hydrogen abstraction reaction of camphor in P450cam catalysis: a QM/MM study.

The discrepancies between the published QM/MM studies (Schöneboom, J. C.; Cohen, S.; Lin, H.; Shaik, S.; Thiel, W. J. Am. Chem. Soc. 2004, 126, 4017; Guallar, V.; Friesner, R. A. J. Am. Chem. Soc. 2004, 126, 8501) on H-abstraction of camphor in P450cam have largely been resolved. The crystallographic water molecule 903 situated near the oxo atom of Compound I acts as a catalyst for H-abstraction, lowering the barrier by about 4 kcal/mol. Spin density at the A-propionate side chain of heme can occur in the case of incomplete screening but has no major effect on the computed barrier.

Camphor↗

PubMed2006

The role of the heme propionates in heme biochemistry.

There are numerous studies, relying on both experimental and theoretical observations, illustrating the active role of the heme propionates in regulating electron delivery to the iron center as well as biochemical properties of the heme. Evidences for this come from a wide variety of heme containing systems: cytochromes, heme peroxidases, globins, etc. Here, we shortly summarize these studies and revisit previous theoretical calculations (V. Guallar, M.H. Baik, S.J. Lippard, R.A. Friesner, Proc. Natl. Acad. Sci. USA 100 (2003) 6998-7002) where the propionate groups induced the delocalization of the spin density in the cytochrome P450cam putative active species, Compound I. We introduce novel data, obtained by means of mixed quantum mechanics and molecular mechanics methods, indicating a larger electron delocalization into the protein. We also present novel results based on the recent migration of spin density observed by Barrows et al. (T.P. Barrows, T.L. Poulos, Biochemistry 44 (2005) 14062-68) on an ascorbate peroxidase mutant. All this data strongly supports the importance of the propionate groups in tuning the heme electronic properties.

Ascorbate Peroxidases↗

PubMed2005

Molecular dynamics simulations of salicylate effects on the micro- and mesoscopic properties of a dipalmitoylphosphatidylcholine bilayer.

Salicylate, an amphiphilic molecule and a popular member of the nonsteroidal anti-inflammatory drug family, is known to affect hearing through reduction of the electromechanical coupling in the outer hair cells of the ear. This reduction of electromotility by salicylate has been widely studied, but the molecular mechanism of the phenomenon is still unknown. In this study, we investigated one aspect of salicylate's action, namely the perturbation of electrical and mechanical membrane properties by salicylate in the absence of cytoskeletal or membrane-bound motor proteins such as prestin. In particular, we simulated the interaction of salicylate with a dipalmitoylphosphatidylcholine (DPPC) bilayer via atomically detailed molecular dynamics simulations to observe the effect of salicylate on the microscopic and mesoscopic properties of the bilayer. The results demonstrate that salicylate interacts with the bilayer by associating at the water-DPPC interface in a nearly perpendicular orientation and penetrating more deeply into the bilayer than either sodium or chloride. This association has several affects on the membrane properties. First, binding of salicylate to the membrane displaces chloride from the bilayer-water interface. Second, salicylate influences the electrostatic potential and dielectric properties of the bilayer, with significant changes at the water-lipid bilayer interface. Third, salicylate association results in structural changes, including decreased headgroup area per lipid and increased lipid tail order. However, salicylate does not significantly alter the mechanical properties of the DPPC bilayer; bulk compressibility, area compressibility, and bending modulus were only perturbed by small, statistically insignificant amounts by the presence of salicylate. The observations from these simulations are in qualitative agreement with experimental data and support the conclusion that salicylate influences the electrical but not the mechanical properties of DPPC membranes.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2005

Importance of accurate charges in molecular docking: quantum mechanical/molecular mechanical (QM/MM) approach.

The extent to which accuracy of electric charges plays a role in protein-ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed "Survival of the Fittest" (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein-ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications.

Algorithms↗

PubMed2005

A binding mechanism in protein-nucleotide interactions: implication for U1A RNA binding.

We present a close electronic view of the protein-base interface for the N-terminal domain of the human protein U1A. Combining accurate mixed quantum mechanics/molecular mechanics techniques and protein structure prediction methods, we provide a detailed electronic structure description of the protein-RNA stacking interactions. Our analysis indicates the evolution of the protein structure optimizing the interaction between Asp-92 and the RNA bases. The results show a direct coupling of the C-terminal tail and Asp-92, providing a direct rationalization of the experimentally determined role of the C-terminal domain in RNA binding. Here, we propose a mechanism where a protein side chain, with a delocalized electronic pi system, assists in the nucleotide binding. The binding mechanism involves a short-range interaction of the side chain with the nucleotide base and an electronic long-range interaction through a sandwich-stacking motif. The structural motif of the binding mechanism is observed in similar protein-RNA interactions and in various protein-ATP-binding sites.

Algorithms↗

PubMed2005

Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis.

We describe large scale ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic reactions. First, technical aspects of the methodology are reviewed, including the hybrid density functional theory (DFT) methods that are typically employed for the QM aspect of the calculations, and various approaches to defining the interface between the QM and MM regions in QM/MM approaches. The modeling of the enzymatic catalytic cycle for three examples--methane monooxygenase, cytochrome P450, and triose phosphate isomerase--are discussed in some depth, followed by a brief summary of other systems that have been investigated by ab initio methods over the past several years. Finally, a discussion of the qualitative and quantitative conclusions concerning enzymatic catalysis that are available from modern ab initio approaches is presented, followed by a conclusion briefly summarizing future prospects.

Antibodies↗

PubMed2004

Cytochrome P450CAM enzymatic catalysis cycle: a quantum mechanics/molecular mechanics study.

The catalytic pathway of cytochrome P450cam is studied by means of a hybrid quantum mechanics/molecular mechanics method. Our results reveal an active role of the enzyme in the different catalytic steps. The protein initially controls the energy gap between the high- and low-spin states in the substrate binding process, allowing thermodynamic reduction by putidaredoxin reductase and molecular oxygen addition. A second electron reduction activates the delivery of protons to the active site through a selective interaction of Thr252 and the distal oxygen causing the O--O cleavage. Finally, the protein environment catalyzes the substrate hydrogen atom abstraction step with a remarkably low free energy barrier ( approximately 8 kcal/mol). Our results are consistent with the effect of mutations on the enzymatic efficacy and provide a satisfactory explanation for the experimental failure to trap the proposed catalytically competent species, a ferryl Fe(IV) heme.

Camphor↗

PubMed2004

Computational modeling of the catalytic reaction in triosephosphate isomerase.

We present a comprehensive analysis of the catalytic cycle of the enzyme triosephosphate isomerase (TIM), including both the reactive chemistry and the catalytic loop and side-chain motions. Combining accurate mixed quantum mechanics/molecular mechanics (QM/MM) and protein structure prediction methods, we have modeled both the structural and chemical aspects of the reversible isomerization of dihydroxyacetone phosphate (DHAP) to d-glyceraldehyde 3-phosphate (GAP), for which there is a wealth of experimental data. The conjunction of this novel computational approach with the use of the recent near-atomic resolution TIM-DHAP Michaelis complex PDB structure, 1NEY.pdb, has enabled us to obtain robust qualitative and, where available, quantitative agreement with a wide range of experimental data. Among the principal conclusions that we are able to draw are the importance of the monoanionic (as opposed to dianioic) form of the substrate phosphate group in the catalytic cycle, detailed positioning and energetics of the key catalytic residues in the active-site, the flexible nature of Glu165, which favors its direct involvement in the formation of the enediol intermediate, energetics of the open and closed form of the catalytic loop region in the presence and absence of substrate, and quantitative reproduction of various experimentally measured reaction rates, typically to within approximately 1 kcal/mol. Our results are consistent with the available experimental data, and provide an initial picture as to why loop opening when GAP is the product has a higher barrier than when DHAP is the product.

Algorithms↗

PubMed2003

Peripheral heme substituents control the hydrogen-atom abstraction chemistry in cytochromes P450.

We elucidate the hydroxylation of camphor by cytochrome P450 with the use of density functional and mixed quantum mechanics/molecular mechanics methods. Our results reveal that the enzyme catalyzes the hydrogen-atom abstraction step with a remarkably low free-energy barrier. This result provides a satisfactory explanation for the experimental failure to trap the proposed catalytically competent high-valent heme Fe(IV) oxo (oxyferryl) species responsible for this hydroxylation chemistry. The primary and previously unappreciated contribution to stabilization of the transition state is the interaction of positively charged residues in the active-site cavity with carboxylate groups on the heme periphery. A similar stabilization found in dioxygen binding to hemerythrin, albeit with reversed polarity, suggests that this mechanism for controlling the relative energetics of redox-active intermediates and transition states in metalloproteins may be widespread in nature.

Camphor↗

PubMed2002

Dynamics of alkane hydroxylation at the non-heme diiron center in methane monooxygenase.

Semiclassical molecular dynamics simulations have been combined with quantum chemistry calculations to provide detailed modeling of the methane and ethane hydroxylation reactions catalyzed by the hydroxylase enzymes of the soluble methane monooxygenase system. The experimental distribution of enantiomeric alcohols in the reaction of ethanes made chiral by the use of hydrogen isotopes is quantitatively reproduced and explained. The reaction dynamics involve a mixture of concerted and bound radical trajectories, and we characterize each of these reactive channels in detail. Diffusion of the bound radical intermediate at the active site core determines the global rate constant. The results also provide a qualitative rationale for the lack of ring-opened products derived from certain radical clock substrate probes and for the relative rate constants and kinetic isotope effects exhibited by a variety of substrates.

Binding Sites↗

PubMed2002

Proton-transfer dynamics in the activation of cytochrome P450eryF.

Molecular dynamics simulations are combined with quantum chemistry calculations of instantaneous proton-transfer energy profiles to investigate proton-transfer events in the transient pathway of cytochrome P450eryF (6-deoxyerythronolide B hydroxylase; CYP107A1), from the oxyferrous species to the catalytically active ferryl oxygen species (compound I). This reaction is one of the most fundamental unresolved aspects in the mechanism of oxidation that is common to all cytochrome P450s. We find that this process involves an ultrafast proton transfer from the crystallographic water molecule W519 to the distal oxygen bound to the heme group, and a subsequent proton-transfer event from W564 to W519. Both proton-transfer events are found to be endothermic in the oxyferrous state, suggesting that the oxyferrous reduction is mechanistically linked to the proton-transfer dynamics. These findings indicate that the hydrogen bond network, proximate to the O(2)-binding cleft, plays a crucial functional role in the enzymatic activation of P450s. Our results are consistent with the effect of mutations on the enzymatic efficacy.

Bacterial Proteins↗

PubMed2002

Quantum chemical studies of methane monooxygenase: comparision with P450.

The catalytic pathways of soluble methane monooxygenase (sMMO) and cytochrome P450CAM, iron-containing enzymes, are described and compared. Recent extensive density functional ab initio electronic structure calculations have revealed many similarities in a number of the key catalytic steps, as well as some important differences. A particularly interesting and significant contrast is the role played by the protein in each system. For sMMO, the protein stabilizes various species in the catalytic cycle through a series of carboxylate shifts. This process is adequately described by a relatively compact model of the active site ( approximately 100 atoms), providing a reasonable description of the energetics of hydrogen atom abstraction. For P450CAM, in contrast, the inclusion of the full protein is necessary for an accurate description of the hydrogen atom abstraction.

Binding Sites↗