Allosteric regulation of glycogen synthase kinase 3ß: a theoretical study.

Glycogen synthase kinase 3ß (GSK-3ß) is a serine-threonine kinase belonging to the CMGC family that plays a key role in many biological processes, such as glucose metabolism, cell cycle regulation, and proliferation. Like most protein kinases, GSK-3ß is regulated via multiple pathways and sites. We performed all-atom molecular dynamics simulations on the unphosphorylated and phosphorylated unbound GSK-3ß and the phosphorylated GSK-3ß bound to a peptide substrate, its product, and a derived inhibitor. We found that GSK-3ß autophosphorylation at residue Tyr(216) results in widening of the catalytic groove, thereby facilitating substrate access. In addition, we studied the interactions of the phosphorylated GSK-3ß with a substrate and peptide inhibitor located at the active site and observed higher affinity of the inhibitor to the kinase. Furthermore, we detected a potential remote binding site which was previously identified in other kinases. In agreement with experiments we observed that binding of specific peptides at this remote site leads to stabilization of the activation loop located in the active site. We speculate that this stabilization could enhance the catalytic activity of the kinase. We point to this remote site as being structurally conserved and suggest that the allosteric phenomenon observed here may occur in the protein kinase superfamily.

How does the reductase help to regulate the catalytic cycle of cytochrome P450 3A4 using the conserved water channel?

Water molecules play a major role in the P450 catalytic cycle. Here, we locate the preferred water pathways and their gating mechanisms for the human cytochrome P450 3A4 (CYP3A4) and elucidate the role of the cytochrome P450 reductase (CPR) in turning on and activating these water channels. We perform explicit solvent molecular dynamic simulations of CYP3A4, unbound and bound to two substrates, and with and without the flavin mononucleotide (FMN)-binding domain of CPR. We observe in/out passage of water molecules via a water-specific and conserved channel (aqueduct) located between the active site and the heme proximal side. We find that the aqueduct gating mechanism is mediated by R375, the conserved arginine that salt bridges with the heme 7-propionate. When R375 rotates, it opens the aqueduct and establishes a connection between a cluster of active site water molecules network and the bulk solvent. The aqueduct region overlaps with the CPR binding-site to CYP3A4. Indeed, we find that when the FMN domain of CPR binds to CYP3A4, the aqueduct fully opens up, thereby allowing a flow of water molecules. The aqueduct's opening can permit proton transfer, shuttling the protons to the active site through ordered water molecules. In addition, the expulsion of water molecules via the aqueduct contributes to substrate binding. As such, the CPR binding has a function: it triggers the aqueduct's opening and thereby enables a proton shuttle pathway, which is needed for the dioxygen activation. This mechanism could be a general paradigm in P450s.

Theoretical characterization of substrate access/exit channels in the human cytochrome P450 3A4 enzyme: involvement of phenylalanine residues in the gating mechanism.

The human cytochrome P450 3A4 mono-oxygenates approximately 50% of all drugs. Its substrates/products enter/leave the active site by access/exit channels. Here, we perform steered molecular dynamics simulations, pulling the products temazepam and testosterone-6betaOH out of the P450 3A4 enzyme in order to identify the preferred substrate/product pathways and their gating mechanism. We locate six different egress pathways of products from the active site with different exit preferences for the two products and find that there is more than just one access/exit channel in CYP3A4. The so-called solvent channel manifests the largest opening for both tested products, thereby identifying this channel as a putative substrate channel. Most channels consist of one or two pi-stacked phenylalanine residues that serve as gate keepers. The oxidized drug breaks the hydrophobic interactions of the gating residues and forms mainly hydrophobic contacts with the gate. We argue that product exit preferences in P450s are regulated by protein-substrate specificity.

MolAxis: a server for identification of channels in macromolecules.

MolAxis is a freely available, easy-to-use web server for identification of channels that connect buried cavities to the outside of macromolecules and for transmembrane (TM) channels in proteins. Biological channels are essential for physiological processes such as electrolyte and metabolite transport across membranes and enzyme catalysis, and can play a role in substrate specificity. Motivated by the importance of channel identification in macromolecules, we developed the MolAxis server. MolAxis implements state-of-the-art, accurate computational-geometry techniques that reduce the dimensions of the channel finding problem, rendering the algorithm extremely efficient. Given a protein or nucleic acid structure in the PDB format, the server outputs all possible channels that connect buried cavities to the outside of the protein or points to the main channel in TM proteins. For each channel, the gating residues and the narrowest radius termed 'bottleneck' are also given along with a full list of the lining residues and the channel surface in a 3D graphical representation. The users can manipulate advanced parameters and direct the channel search according to their needs. MolAxis is available as a web server or as a stand-alone program at http://bioinfo3d.cs.tau.ac.il/MolAxis.

MolAxis: efficient and accurate identification of channels in macromolecules.

Channels and cavities play important roles in macromolecular functions, serving as access/exit routes for substrates/products, cofactor and drug binding, catalytic sites, and ligand/protein. In addition, channels formed by transmembrane (TM) proteins serve as transporters and ion channels. MolAxis is a new sensitive and fast tool for the identification and classification of channels and cavities of various sizes and shapes in macromolecules. MolAxis constructs corridors, which are pathways that represent probable routes taken by small molecules passing through channels. The outer medial axis of the molecule is the collection of points that have more than one closest atom. It is composed of two-dimensional surface patches and can be seen as a skeleton of the complement of the molecule. We have implemented in MolAxis a novel algorithm that uses state-of-the-art computational geometry techniques to approximate and scan a useful subset of the outer medial axis, thereby reducing the dimension of the problem and consequently rendering the algorithm extremely efficient. MolAxis is designed to identify channels that connect buried cavities to the outside of macromolecules and to identify TM channels in proteins. We apply MolAxis to enzyme cavities and TM proteins. We further utilize MolAxis to monitor channel dimensions along Molecular Dynamics trajectories of a human Cytochrome P450. MolAxis constructs high quality corridors for snapshots at picosecond time-scale intervals substantiating the gating mechanism in the 2e substrate access channel. We compare our results with previous tools in terms of accuracy, performance and underlying theoretical guarantees of finding the desired pathways. MolAxis is available on line as a web-server and as a stand alone easy-to-use program (http://bioinfo3d.cs.tau.ac.il/MolAxis/).

QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding.

Cytochrome P450 3A4 is involved in the metabolism of 50% of all swallowed drugs. The enzyme functions by means of a high-valent iron-oxo species, called compound I (Cpd I), which is formed after entrance of the substrate to the active site. We explored the features of Cpd I using hybrid quantum mechanical/molecular mechanical calculations on various models that are either substrate-free or containing one and two molecules of diazepam as a substrate. Mössbauer parameters of Cpd I were computed. Our major finding shows that without the substrate, Cpd I tends to elongate its Fe-S bond, localize the radical on the sulfur, and form hydrogen bonds with A305 and T309, which may hypothetically lead to Cpd I consumption by H-abstraction. However, the positioning of diazepam close to Cpd I, as enforced by the effector molecule, was found to strengthen the NH...S interactions of the conserved I443 and G444 residues with the proximal cysteinate ligand. These interactions are known to stabilize the Fe-S bond, and as such, the presence of the substrate leads to a shorter Fe-S bond and it prevents the localization of the radical on the sulfur. This diazepam-Cpd I stabilization was manifested in the 1W0E conformer. The effector substrate did not influence Cpd I directly but rather by positioning the active substrate close to Cpd I, thus displacing the hydrogen bonds with A305 and T309, and thereby giving preference to substrate oxidation. It is hypothesized that these effects on Cpd I, promoted by the restrained substrate, may be behind the special metabolic behavior observed in cases of multiple substrate binding (also called cooperative binding). This restraint constitutes a mechanism whereby substrates stabilize Cpd I sufficiently long to affect monooxygenation by P450s at the expense of Cpd I destruction by the protein residues.

Structural dynamics of the cooperative binding of organic molecules in the human cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is a key enzyme responsible for the metabolism of 50% of all orally administered drugs which exhibit an intriguing kinetic behavior typified by a sigmoidal dependence of the reaction velocity on the substrate concentration. There is evidence for the binding of two substrates in the active site of the enzyme, but the mechanism of this cooperative binding is unclear. Diazepam is such a drug that undergoes metabolism by CYP3A4 with sigmoidal dependence. Metabolism is initiated by hydrogen atom abstraction from the drug. To understand the factors that determine the cooperative binding and the juxtaposition of the C-H bond undergoing abstraction, we carried out molecular dynamics simulations for two enzymatic conformers and examined the differences between the substrate-free and the bound enzymes, with one and two diazepam molecules. Our results indicate that the effector substrate interacts both with the active substrate and with the enzyme, and that this interaction results in side chain reorientation with relatively minor long-range effects. In accord with experiment, we find that F304, in the interface between the active and effector binding sites, is a key residue in the mechanism of cooperative binding. The addition of the effector substrate stabilizes F304 and its environment, especially F213, and induces a favorable orientation of the active substrate, leading to a short distance between the targeted hydrogen for abstraction and the active species of the enzyme. In addition, in one conformer of the enzyme, residue R212 may strongly interact with F304 and counteract the effector's impact on the enzyme.