PubMed Health⌕ Search

Biomedical subjects

Jerome Baudry

Publications and source records attributed to Jerome Baudry.

13 recordsLinked to original sources

Van der waals interactions and decrease of the rotational barrier of methyl-sized rotators: a theoretical study.

Rotational barriers of methyl-sized molecular rotators are investigated theoretically using ab initio and empirical force field calculations in molecular models simulating various environmental conditions experienced by the molecular rotors. Calculations on neopentane surrounded by methyl groups suggest that the neopentane's methyl rotational potential energy barrier can be reduced by up to an order of magnitude by locating satellite functional groups around the rotator at a geometry that destabilizes the staggered conformation of the rotator through van der Waals repulsive interactions and reduces the staggered/eclipsed relative energy difference. Molecular mechanics and molecular dynamics calculations indicate that this barrier-reducing geometry can also be found in molecular rotators surface mounted on graphite surfaces or carbon nanotube models. In these models, molecular dynamics simulations show that the rotation of methyl-sized functional groups can be catalyzed by van der Waals interactions, thus making very rigid rotators become thermally activated at room temperature. These results are discussed in the context of design of nanostructures and use of methyl groups as markers for microenvironmental conditions.

Journal Article↗

Class-dependent sequence alignment strategy improves the structural and functional modeling of P450s.

Different procedures for obtaining homology models for P450s are investigated using various sequence alignments sharing various levels of sequence identity with available P450 crystal structures. In this analysis, we have investigated how well homology modeling can reproduce known crystal structures as well as how effectively these homology models can be used to reproduce known ligand-binding modes. Homology models obtained from sequence alignments that discriminate between Class I and Class II P450s are significantly closer to the experimental crystal structures and more closely reproduce known ligand's binding modes, than those obtained using sequence alignments that combine Class I and Class II P450s. The quality of the models is slightly improved by constructing hybrid-structure models that model three of the most variable regions of P450s independently from the rest of the protein: the B region that includes SRS1, the FG region that includes SRS2 and SRS3 and the beta4 region that includes SRS6.

Animals↗

Can proteins and crystals self-catalyze methyl rotations?

The chi (C(alpha)-C(beta)) torsional barrier in the dipeptide alanine (N-methyl-l-alanyl-N-methylamide) crystal was investigated using ab initio calculations at various levels of theory, molecular mechanics, and molecular dynamics. For one of the two molecules in the asymmetric unit the calculations suggest that rotation around the chi dihedral angle is catalyzed by the crystal environment, reducing by up to approximately 2kT the torsional barrier in the crystal with respect to that in the gas phase. This catalytic effect is present at both low and room temperature and originates from a van der Waals destabilization of the minima in the methyl dihedral potential coming from the nonbonded environment of the side chain. Screening of a subset of the Protein Data Bank with a pharmacophore model reproducing the crystal environment around this side chain methyl identified a protein containing an alanine residue with an environment similar to that in the crystal. Calculations indicate that this chi torsional barrier is also reduced in the protein at low temperature but not at room temperature. This suggests that environment-catalyzed rotation of methyl groups can occur both in the solid phase and in native biological structures, though this effect might be temperature-dependent. The relevance of this catalytic effect is discussed in terms of its natural occurrence and its possible contribution to the low-frequency vibrational modes of molecules.

Catalysis↗

Ile115Leu mutation in the SRS1 region of an insect cytochrome P450 (CYP6B1) compromises substrate turnover via changes in a predicted product release channel.

CYP6B1 represents the principal cytochrome P450 monooxygenase responsible for metabolizing furanocoumarins in Papilio polyxenes, an insect that specializes on host plants containing these toxins. Investigations of the amino acids responsible for the efficient metabolism of these plant toxins has identified Ile115 as one that modulates the rate of furanocoumarin metabolism even though it is predicted to be positioned at the edge of the heme plane and outside substrate contact regions. In contrast to previous expression studies conducted under conditions of limiting P450 reductase showing that the Ile115-to-Leu replacement enhances turnover of xanthotoxin and other furanocoumarins, studies conducted at high P450 reductase indicate that the Ile115-to-Leu replacement reduces turnover of these substrates. Further analysis of substrate binding affinities, heme spin state and NADPH consumption rates indicate that, whereas the I115L replacement mutant displays higher substrate affinity and heme spin state than the wild-type CYP6B1 protein, it utilizes NADPH more slowly than the wild-type CYP6B1 protein at high P450 reductase levels. Molecular models developed for the wild-type CYP6B1 and mutant protein suggest that more constricted channels extending from the catalytic site in the I115L mutant to the P450 surface limit the rate of product release from this mutant catalytic site under conditions not limited by the rate of electron transfer from NADPH.

Amino Acid Sequence↗

A retinoic acid binding cytochrome P450: CYP120A1 from Synechocystis sp. PCC 6803.

At least 35 cytochrome P450 (P450, CYP) or cytochrome P450-like genes have been identified in 10 cyanobacterial genomes yet none have been functionally characterized. CYP110 and CYP120 represent the two largest cyanobacterial P450 families with 16 and four members, respectively, identified to date. The Synechocystis sp. PCC 6803 CYP120A1 protein sequence shares high degrees of conservation with CYP120A2 from Trichodesmium erythraeum IMS101 and CYP120B1 and CYP120C1 from Nostoc punctiforme PCC 73102. In this communication, we report the cloning, expression, purification, and characterization of CYP120A1 from Synechocystis. Homology modeling predictions of the three-dimensional structure of CYP120A1 coupled with in silico screening for potential substrates and experimental spectroscopic analyses have identified retinoic acid as a compound binding with high affinity to this P450's catalytic site. These characterizations of Synechocystis CYP120A1 lay the initial foundations for understanding the basic role of cytochrome P450s in cyanobacteria and related organisms.

Bacterial Proteins↗

Structural and functional divergence of insect CYP6B proteins: From specialist to generalist cytochrome P450.

How polyphagous herbivores cope with the diversity and unpredictability of plant defenses remains largely unknown at both the genetic and molecular levels. To examine whether generalist counterdefense enzymes are structurally more flexible and functionally more diverse, two counterdefensive allelochemical-metabolizing cytochrome P450 proteins, CYP6B1 from the specialist Papilio polyxenes, feeding on furanocoumarin-containing plants, and CYP6B8 from the generalist Helicoverpa zea, feeding on hundreds of host plant species, are compared structurally and functionally. Molecular modeling indicates that CYP6B8 has more flexible overall folding, a more elastic catalytic pocket, and one more substrate access channel than CYP6B1. Baculovirus-mediated expression of the CYP6B8 and CYP6B1 proteins demonstrates that CYP6B8 metabolizes six biosynthetically diverse plant allelochemicals (xanthotoxin, quercetin, flavone, chlorogenic acid, indole-3-carbinol, and rutin) and three insecticides (diazinon, cypermethrin, and aldrin), whereas CYP6B1 metabolizes only two allelochemicals (xanthotoxin and flavone) and one insecticide (diazinon) of those tested. These results indicate that generalist counterdefense proteins are capable of accepting a more structurally diverse array of compounds compared with specialist counterdefense proteins.

Amino Acid Sequence↗

Identification of variable amino acids in the SRS1 region of CYP6B1 modulating furanocoumarin metabolism.

The homology model of Papilio polyxenes CYP6B1 places Ile115, one of two variable amino acids, in the SRS1 of various CYP6B subfamily proteins in close proximity to the heme and Ala113, another variable amino acid, in a more distal position. We have constructed mutant CYP6B1 proteins altered at either of these positions and homology models of each based on multiple alignments with crystallized P450 proteins. The homology models suggest the existence of significant structural diversity in the hydrogen bond network surrounding the heme as a result of single point mutations in SRS1. Mutagenesis of Ile115 or Ala113 to other residues present in the insect CYP6B subfamily indicates that these amino acids control the spin state of the heme and, as a result, the catalytic activity of this monooxygenase. In particular, the I115L mutation significantly increases the spin state of the heme coordinately with 2- to 4-fold increases in its turnover of linear furanocoumarins. Other A113V, A113L, A113Q, and A113E mutations display more variation in their effects but, in each case, strong correlations exist between furanocoumarin turnover and heme spin state. These data demonstrate that variable amino acids in SRS1 of the insect CYP6B subfamily exert dramatic effects on the range of furanocoumarins metabolized, even when they occur in positions potentially distal from the substrate. These effects are possibly mediated through rearrangement of the local hydrogen bond network.

Amino Acid Sequence↗

Classical force field parameters for the heme prosthetic group of cytochrome c.

Accurate force fields are essential for describing biological systems in a molecular dynamics simulation. To analyze the docking of the small redox protein cytochrome c (cyt c) requires simulation parameters for the heme in both the reduced and oxidized states. This work presents parameters for the partial charges and geometries for the heme in both redox states with ligands appropriate to cyt c. The parameters are based on both protein X-ray structures and ab initio density functional theory (DFT) geometry optimizations at the B3LYP/6-31G* level. The simulations with the new parameter set reproduce the geometries of the X-ray structures and the interaction energies between water and heme prosthetic group obtained from B3LYP/6-31G* calculations. The parameter set developed here will provide new insights into docking processes of heme containing redox proteins.

Algorithms↗

Molecular docking of substrates and inhibitors in the catalytic site of CYP6B1, an insect cytochrome p450 monooxygenase.

Furanocoumarins represent plant toxins that are used in the treatment of a variety of skin diseases and are metabolized by cytochrome p450 monooxygenases (p450s) existing in insects such as Papilio polyxenes (the black swallowtail). To elucidate the active site in the CYP6B1 protein that is the principal p450 existing in this species, we have constructed a homology model of it based on sequence and structure alignments with the bacterial CYP102 protein whose crystal structure has been defined and with the insect CYP6B4 protein that also metabolizes furanocoumarins. In the derived CYP6B1 model, Phe116 and His117 in SRS1, Phe371 in SRS5 and Phe484 in SRS6 contribute to the formation of a resonant network that stabilizes the p450's catalytic site and allows for interactions with its furanocoumarin substrates. The first two of these residues are absolutely conserved in all members of the insect CYP6B subfamily and the last two are variable in different members of the CYP6B subfamily. A combination of theoretical and experimental docking analyses of two substrates (xanthotoxin and bergapten) and two inhibitors (coumarin and pilocarpine) of this p450 provide significant information on the positioning of furanocoumarins within this catalytic pocket. Molecular replacement models based on the results of variations at two of these critical amino acids provide support for our furanocoumarin-docked model and begin to rationalize the altered substrate reactivities observed in experimental analyses.

5-Methoxypsoralen↗

Common active site architecture and binding strategy of four phenylpropanoid P450s from Arabidopsis thaliana as revealed by molecular modeling.

Despite extensive primary sequence diversity, crystal structures of several bacterial cytochrome P450 monooxygenases (P450s) and a single eukaryotic P450 indicate that these enzymes share a structural core of alpha-helices and beta-sheets and vary in the loop regions contacting individual substrates. To determine the extent to which individual structural features are conserved among divergent P450s existing in a single biosynthetic pathway, we have modeled the structures of four highly divergent P450s (CYP73A5, CYP84A1, CYP75B1, CYP98A3) in the Arabidopsis phenylpropanoid pathway synthesizing lignins, flavonoids and anthocyanins. Analysis of these models has indicated that, despite primary sequence identities as low as 13%, the structural cores and several loop regions of these P450s are highly conserved. Substrate docking indicated that all four enzymes employ a common strategy to identify their substrates in that their cinnamate-derived substrates align along helix I with their aromatic ring positioned towards the C-terminus of this helix and their aliphatic tails positioned towards the N-terminus. Further similarity was observed in the way the substrates contact the consensus P450 substrate recognition sites (SRS). Residues predicted to contact the aromatic ring region exist in SRS5, SRS6 and the C-terminal portion of SRS4 and residues contacting the distal end of each substrate exist in SRS1, SRS2 and the N-terminal portion of SRS4. Alignments of the regions contacting the aromatic ring region indicate that SRS4, SRS5 and SRS6 share higher degrees of sequence conservation than found in SRS1, SRS2 or the full-length protein.

Amino Acid Sequence↗

A beta, aging, and Alzheimer's disease: a tale, models, and hypotheses.

In this paper we explore the potential functional role of the A beta peptides in the context of Alzheimer's disease (AD). We begin by defining the morphology of the amyloid deposits in relation to surrounding glial cells and, more importantly, in relation to the brain vasculature. Amyloid accumulation in the brain's microvasculature causes disturbances in the blood-brain barrier (BBB), and in larger arteries, impairment in control of regional cerebral blood flow due to myocyte degeneration. We postulate that the deposition of vascular amyloid may represent a hydrophobic protein plaster to seal leaks in the BBB, occasionally observed in aging and catastrophically common in AD. The vasoconstrictive activity of A beta may also be related to leaky vessels whereby decreasing the arterial diameter may also help to control breaches in the BBB. The admission of plasma neurotoxic proteins into the brain may be controlled by activation of microglia elicited by soluble A beta peptides creating a subtle, but permanent brain inflammatory reaction. We also delve into the influence that cholesterol metabolism may have in membrane topology and A beta production, and the close correlations that exist between cardiovascular disease and AD. Finally, we speculate about the possibility of a peripheral source of A beta that may, by crossing the BBB, contribute to the vascular and parenchymal deposits of A beta in the AD brain.

Aging↗

Molecular basis for the subtype discrimination of the estrogen receptor-beta-selective ligand, diarylpropionitrile.

Although the two subtypes of the human estrogen receptor (ER), ERalpha and ERbeta, share only 56% amino acid sequence identity in their ligand binding domain (LBD), the residues that surround the ligand are nearly identical; nevertheless, subtype-selective ligands are known. To understand the molecular basis by which diarylpropionitrile (DPN), an ERbeta-selective ligand, is able to discriminate between the two ERs, we examined its activity on ER mutants and chimeric constructs generated by DNA shuffling. The N-terminal region of the ERbeta LBD (through helix 6) appears to be fully responsible for the ERbeta selectivity of DPN. In fact, a single ERalpha point mutation (L384M) was largely sufficient to switch the DPN response of this ER to that of the ERbeta type, but residues in helix 3 are also important in achieving the full ERbeta selectivity of DPN. Using molecular modeling, we found an energetically favorable fit for the S-DPN enantiomer in ERbeta, in which the proximal phenol mimics the A ring of estradiol, and the nitrile engages in stabilizing interactions with residues in the ligand-binding pocket of ERbeta. Our findings highlight that a limited number of critical interactions of DPN with the ERbeta ligand-binding pocket underlie its ER subtype-selective character.

Amino Acid Sequence↗