Qing-Chuan Zheng

Structural and Dynamic Basis of Human Cytochrome P450 7B1: A Survey of Substrate Selectivity and Major Active Site Access Channels

Cytochrome P450 (CYP) 7B1 is a steroid cytochrome P450 7α-hydroxylase that has been linked directly with bile salt synthesis and hereditary spastic paraplegia type 5 (SPG5). The enzyme provides the primary metabolic route for neurosteroids dehydroepiandrosterone (DHEA), cholesterol derivatives 25-hydroxycholesterol (25-HOChol), and other steroids such as 5α-androstane-3β,17β-diol (anediol), and 5α-androstene-3β,17β-diol (enediol). A series of investigations including homology modeling, molecular dynamics (MD), and automatic docking, combined with the results of previous experimental site-directed mutagenesis studies and access channels analysis, have identified the structural features relevant to the substrate selectivity of CYP7B1. The results clearly identify the dominant access channels and critical residues responsible for ligand binding. Both binding free energy analysis and total interaction energy analysis are consistent with the experimental conclusion that 25-HOChol is the best substrate. According to 20 ns MD simulations, the Phe cluster residues that lie above the active site, particularly Phe489, are proposed to merge the active site with the adjacent channel to the surface and accommodate substrate binding in a reasonable orientation. The investigation of CYP7B1-substrate binding modes provides detailed insights into the poorly understood structural features of human CYP7B1 at the atomic level, and will be valuable information for drug development and protein engineering.

Molecular Dynamics Simulations Suggest Ligand’s Binding to Nicotinamidase/Pyrazinamidase

The research on the binding process of ligand to pyrazinamidase (PncA) is crucial for elucidating the inherent relationship between resistance of Mycobacterium tuberculosis and PncA’s activity. In the present study, molecular dynamics (MD) simulation methods were performed to investigate the unbinding process of nicotinamide (NAM) from two PncA enzymes, which is the reverse of the corresponding binding process. The calculated potential of mean force (PMF) based on the steered molecular dynamics (SMD) simulations sheds light on an optimal binding/unbinding pathway of the ligand. The comparative analyses between two PncAs clearly exhibit the consistency of the binding/unbinding pathway in the two enzymes, implying the universality of the pathway in all kinds of PncAs. Several important residues dominating the pathway were also determined by the calculation of interaction energies. The structural change of the proteins induced by NAM’s unbinding or binding shows the great extent interior motion in some homologous region adjacent to the active sites of the two PncAs. The structure comparison substantiates that this region should be very important for the ligand’s binding in all PncAs. Additionally, MD simulations also show that the coordination position of the ligand is displaced by one water molecule in the unliganded enzymes. These results could provide the more penetrating understanding of drug resistance of M. tuberculosis and be helpful for the development of new antituberculosis drugs.

Exploring the molecular basis of dsRNA recognition by Mss116p using molecular dynamics simulations and free-energy calculations.

DEAD-box proteins are the largest family of helicase that are important in nearly all aspects of RNA metabolism. However, it is unclear how these proteins recognize and bind RNA. Here, we present a detailed analysis of the related DEAD-box protein Mss116p-RNA interaction, using molecular dynamics simulations with MM-GBSA calculations. The energetic analysis indicates that the two strands of double strands RNA (dsRNA) are recognized asymmetrically by Mss116p. The strand 1 of dsRNA provides the main binding affinity. Meanwhile, the nonpolar interaction provides the main driving force for the binding process. Although the contribution of polar interaction is small, it is vital in stabilizing the protein-RNA interaction. Compared with the wild type Mss116p, two studied mutants Q412A and D441A have obviously reduced binding free energies with dsRNA because of the decreasing of polar interaction. Three important residues Lys409, Arg415 and Arg438 lose their binding affinity significantly in mutants. In conclusion, these results complement previous experiments to advance comprehensive understanding of Mss116p-dsRNA interaction. The results also would provide support for the application of similar approaches to the understanding of other DEAD-box protein-RNA complexes.

A comparative analysis of binding sites between mouse CYP2C38 and CYP2C39 based on homology modeling, molecular dynamics simulation and docking studies.

Mouse CYP2C38 and CYP2C39 are two closely related enzymes with 91.8% sequence identity. But they exhibit different substrate binding features. In this study, three-dimensional models of CYP2C38 and CYP2C39 were constructed using X-ray crystal structure of human CYP2C8 as the template based on homology modeling methods and molecular dynamics simulations. Tolbutamide as the common substrate of CYP2C38 and CYP2C39 was docked into them and positioned in their active sites with different orientation. All-trans retinoic acid (atRA) is a specific substrate for CYP2C39 and not catalyzed by CYP2C38. By comparison of active site architectures between CYP2C38 and CYP2C39, the possible reasons affecting their substrate binding were proposed. In addition, Arg241, Glu300, Leu366 and Leu476 are identified as critical residue for substrates binding.

Improving GPX activity of selenium‐containing human single‐chain Fv antibody by site‐directed mutation based on the structural analysis

Glutathione peroxidase (GPX) is one of the important members of the antioxidant enzyme family. It can catalyze the reduction of hydroperoxides with glutathione to protect cells against oxidative damage. In previous studies, we have prepared the human catalytic antibody Se‐scFv‐B3 (selenium‐containing single‐chain Fv fragment of clone B3) with GPX activity by incorporating a catalytic group Sec (selenocysteine) into the binding site using chemical mutation; however, its activity was not very satisfying. In order to try to improve its GPX activity, structural analysis of the scFv‐B3 was carried out. A three‐dimensional (3D) structure of scFv‐B3 was constructed by means of homology modeling and binding site analysis was carried out. Computer‐aided docking and energy minimization (EM) calculations of the antibody‐GSH (glutathione) complex were also performed. From these simulations, Ala44 and Ala180 in the candidate binding sites were chosen to be mutated to serines respectively, which can be subsequently converted into the catalytic Sec group. The two mutated protein and wild type of the scFv were all expressed in soluble form in Escherichia coli Rosetta and purified by Ni2+‐immobilized metal affinity chromatography (IMAC), then transformed to selenium‐containing catalytic antibody with GPX activity by chemical modification of the reactive serine residues. The GPX activity of the mutated catalytic antibody Se‐scFv‐B3‐A180S was significantly increased compared to the original Se‐scFv‐B3. Copyright

Molecular basis of the recognition of arachidonic acid by cytochrome P450 2E1 along major access tunnel

Cytochrome P450 2E1 is widely known for its ability to oxidize both low molecular weight xenobiotics and endogenous fatty acids (e.g., arachidonic acid (AA)). In this study, we investigated the structural features of the AA‐bound CYP2E1 complex utilizing molecular dynamics (MD) and found that the distinct binding modes for both AA and fatty acid analog are conserved. Moreover, multiple random acceleration MD simulations and steered MD simulations uncovered the most possible tunnel for fatty acids. The main attractions are derived from three key residues, His107, Ala108, and His109, whose side chains reorient to keep ligands bound via hydrogen bonds during the initial unbinding process. More importantly, based on the calculated binding free energy results, we hypothesize that the hydrogen bonds between the receptor and the ligand are the most important contributors involved in the binding affinity. Thus, it is inferred that the hydrogen bonds between these three residues and the ligand may help offer insights into the structural basis of the different ligand egress mechanisms for fatty acids and small weight compounds. Our investigation provides detailed atomistic insights into the structural features of human CYP2E1–fatty acid complex structures. Furthermore, the ligand‐binding characteristics obtained in the present study are helpful for both experimental and computational studies of CYPs and may allow future researchers to achieve desirable changes in enzymatic activities.

Analysis of clinically relevant substrates of CYP2B6 enzyme by computational methods

Mounting evidence thus far indicates that human cytochrome P450 2B6 (CYP2B6), an enzyme expressed at a relatively low level functionally, is primarily responsible for the metabolism of several clinically relevant drugs, including propofol, efavirenz, bupropion, mephobarbital, and the propofol analog 2,6-di-sec-butyl phenol. We used molecular dynamics and molecular docking methods to predict such interactions and to compare with experimentally measured metabolisms. Insight II and Discover Studio 2.5 were used to carry out the docking of these substrates into CYP2B6 to explore the critical residues and interaction energies of the complexes. Phe297, Glu301, Thr302 and Val367 were identified as major drug-binding residues, which is consistent with previous data on site-directed mutagenesis, crystallography structure, and from modeling and docking studies. In addition, our docking results suggest that nonpolar amino acid clusters and heme also participate in binding to mediate drug oxidative metabolism. The binding modes of the five clinically relevant substrates mentioned above for metabolism on CYP2B6 are presented.

Homology modeling and molecular dynamics study on N-acetylneuraminate lyase

With homology modeling techniques, molecular mechanics and molecular dynamics methods, a 3D structure model of N-acetylneuraminate lyase from human (hNAL, EC 4.1.3.3) was created and refined. This model was further assessed by Profile-3D and PROCHECK, which confirms that the refined model is reliable. Furthermore, the docking results of the substrates (sialic acid and KDO) into the active site of hNAL indicate that hNAL can cleave the sialic acid and KDO. Thr51 and Tyr143 may be the key amino acids residues as they have strong hydrogen bonding interactions with the substrates, which is in good agreement with the experimental results by Izard et al. (Structure 2:361–369. doi:10.1016/S0969-2126(00)00038-1 (1994)). From the docking studies, we also suggest that Asp176 and Ser218 only form hydrogen bonds with sialic acid, therefore, they may help sialic acid interact with hNAL steadly.

How Does (E)-2-(Acetamidomethylene)succinate Bind to Its Hydrolase? From the Binding Process to the Final Result

The binding of (E)-2-(acetamidomethylene)succinate (E-2AMS) to E-2AMS hydrolase is crucial for biological function of the enzyme and the last step reaction of vitamin B6 biological degradation. In the present study, several molecular simulation methods, including molecular docking, conventional molecular dynamics (MD), steered MD (SMD), and free energy calculation methods, were properly integrated to investigate the detailed binding process of E-2AMS to its hydrolase and to assign the optimal enzyme-substrate complex conformation. It was demonstrated that the substrate binding conformation with trans-form amide bond is energetically preferred conformation, in which E-2AMSs pose not only ensures hydrogen bond formation of its amide oxygen atom with the vicinal oxyanion hole but also provides probability of the hydrophobic interaction between its methyl moiety and the related enzymes hydrophobic cavity. Several key residues, Arg146, Arg167, Tyr168, Arg179, and Tyr259, orientate the E-2AMSs pose and stabilize its conformation in the active site via the hydrogen bond interaction with E-2AMS. Sequentially, the binding process of E-2AMS to E-2AMS hydrolase was studied by SMD simulation, which shows the surprising conformational reversal of E-2AMS. Several important intermediate structures and some significant residues were identified in the simulation. It is stressed that Arg146 and Arg167 are two pivotal residues responsible for the conformational reversal of E-2AMS in the binding or unbinding. Our research has shed light onto the full binding process of the substrate to E-2AMS hydrolase, which could provide more penetrating insight into the interaction of E-2AMS with the enzyme and would help in the further exploration on the catalysis mechanism.

Exploring the mechanism how Marburg virus VP35 recognizes and binds dsRNA by molecular dynamics simulations and free energy calculations.

Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled‐stranded RNA to deactivate interferon. The interface of VP35–dsRNA would be a feasible target for structure‐based antiviral agent design. To explore the essence of VP35–dsRNA interaction, molecular dynamics simulation combined with MM‐GBSA calculations were performed on Marburg virus VP35–dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35–dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35–dsRNA recognition mechanism will provide some insights for development of antiviral drug.