Christo Z. Christov

Conformational effects on the circular dichroism of Human Carbonic Anhydrase II: a multilevel computational study.

Circular Dichroism (CD) spectroscopy is a powerful method for investigating conformational changes in proteins and therefore has numerous applications in structural and molecular biology. Here a computational investigation of the CD spectrum of the Human Carbonic Anhydrase II (HCAII), with main focus on the near-UV CD spectra of the wild-type enzyme and it seven tryptophan mutant forms, is presented and compared to experimental studies. Multilevel computational methods (Molecular Dynamics, Semiempirical Quantum Mechanics, Time-Dependent Density Functional Theory) were applied in order to gain insight into the mechanisms of interaction between the aromatic chromophores within the protein environment and understand how the conformational flexibility of the protein influences these mechanisms. The analysis suggests that combining CD semi empirical calculations, crystal structures and molecular dynamics (MD) could help in achieving a better agreement between the computed and experimental protein spectra and provide some unique insight into the dynamic nature of the mechanisms of chromophore interactions.

Conformational Effects on the pro-S Hydrogen Abstraction Reaction in Cyclooxygenase-1: An Integrated QM/MM and MD Study

A key step in the cyclooxygenase reaction cycle of cyclooxygenase 1 (COX-1) is abstraction of the pro-S hydrogen atom of the arachidonic acid by a radical that is formed at the protein residue Tyr-385. Here we investigate this reaction step by a quantum-mechanics/molecular-mechanics approach in combination with molecular-dynamics simulations. The simulations identify the hydrogen abstraction angle as a crucial geometric determinant of the reaction, thus revealing the importance of the cyclooxygenase active site for calculating the potential energy surface of the reaction.

Integrated Synthetic, Pharmacological, and Computational Investigation of cis‐2‐(3,5‐Dichlorophenylcarbamoyl)cyclohexanecarboxylic Acid Enantiomers As Positive Allosteric Modulators of Metabotropic Glutamate Receptor Subtype 4

2‐(3,5‐Dichlorophenylcarbamoyl)cyclohexanecarboxylic acid (1) is a potent and selective positive allosteric modulator of metabotropic glutamate receptor subtype 4 (mGluR4). The activity of 1 was reported to reside in the cis diastereomer with equal potency between its enantiomeric forms (Niswender et al., Mol. Pharmacol. 2008, 74, 1345–1358). In the present study, the asymmetric synthesis of each of the cis enantiomers was performed, and their activities were compared with that of the racemic trans. In our assays, the cis enantiomers differ in potency, with one of them (1R,2S) higher and the other (1S,2R) lower than the racemic trans. High‐level quantum chemical calculations were carried out to characterize the structures of minimum energy in all‐isomer conformational space as well as particular intermediates between conformational transitions. Computational analysis identified structural features of 1 that can play a role in mGluR4 functionality and establish the basis for subsequent work, in which molecular chirality constructed on conformations derived from those found for the active (1R,2S) enantiomer can provide new ideas for drug discovery. Comparison between experimental and theoretical circular dichroism spectra confirmed both the absolute configuration of the (1R,2S) compound and its calculated most stable conformation, thereby supporting experimental and theoretical work.

Mechanisms of protein circular dichroism: insights from computational modeling.

Chirality is a fundamental property of molecular systems, and stereoselectivity underlines many fundamental biomolecular processes like biological recognition and catalysis. Circular dichroism (CD) which is a consequence of molecular chirality is an important method for the investigation of protein structure and structural changes during interactions with ligands, mutations, and folding. The development of computational methods allows powerful insight to be provided into the mechanisms of generation of CD spectra in complex systems as proteins and to explain experimental data, to validate predicted structures, and to explain fine details of biomolecular interactions. In this chapter, we provide a survey on several aspects of the current investigation on the CD phenomena: the emphasis is on its mechanisms and how they can be analyzed using computational methods in strong reference to the experimental data. We analyzed the mechanisms of interactions of the aromatic and disulfide chromophores mainly in the near-UV CD as the number is smaller and allows detailed analysis. We describe some of the investigations on model systems and improve the parameter sets for the approximated CD modeling methods and key calculations on several proteins. We also provide a comprehensive survey on the investigations done in our group on the chiropticity of class A β-lactamases.

Conformational Dynamics, Ligand Binding and Effects of Mutations in NirE an S-Adenosyl-L-Methionine Dependent Methyltransferase

Heme d1, a vital tetrapyrrol involved in the denitrification processes is synthesized from its precursor molecule precorrin-2 in a chemical reaction catalysed by an S-adenosyl-L-methionine (SAM) dependent Methyltransferase (NirE). The NirE enzyme catalyses the transfer of a methyl group from the SAM to uroporphyrinogen III and serves as a novel potential drug target for the pharmaceutical industry. An important insight into the structure-activity relationships of NirE has been revealed by elucidating its crystal structure, but there is still no understanding about how conformational flexibility influences structure, cofactor and substrate binding by the enzyme as well as the structural effects of mutations of residues involved in binding and catalysis. In order to provide this missing but very important information we performed a comprehensive atomistic molecular dynamics study which revealed that i) the binding of the substrate contributes to the stabilization of the structure of the full complex; ii) conformational changes influence the orientation of the pyrrole rings in the substrate, iii) more open conformation of enzyme active site to accommodate the substrate as an outcome of conformational motions; and iv) the mutations of binding and active site residues lead to sensitive structural changes which influence binding and catalysis.

Mechanistic Insights into the Reaction of Chlorination of Tryptophan Catalyzed by Tryptophan 7-Halogenase

Tryptophan 7-halogenase catalyzes chlorination of free tryptophan to 7-chlorotryptophan, which is the first step in the antibiotic pyrrolnitrin biosynthesis. Many biologically and pharmaceutically active natural products contain chlorine and thus, an understanding of the mechanism of its introduction into organic molecules is important. Whilst enzyme-catalyzed chlorination is accomplished with ease, it remains a difficult task for the chemists. Therefore, utilizing enzymes in the synthesis of chlorinated organic compounds is important, and providing atomistic mechanistic insights about the reaction mechanism of tryptophan 7-halogenase is vital and timely. In this work, we examined a mechanism for the reaction of tryptophan chlorination, performed by tryptophan 7-halogenase, by calculating potential energy and free energy surfaces using two different Combined Quantum Mechanical/Molecular Mechanical (QM/MM) methods both employing Density Functional Theory (DFT) for the QM region. Both computational strategies agree on the nature of the rate-limiting step and provided close results for the reaction barriers of the two reaction steps. The calculations for both the potential energy and the free energy profiles showed very similar geometric features and hydrogen bonding interactions for the characterized stationary points.

Effects of Mutations on Structure–Function Relationships of Matrix Metalloproteinase-1

Matrix metalloproteinase-1 (MMP-1) is one of the most widely studied enzymes involved in collagen degradation. Mutations of specific residues in the MMP-1 hemopexin-like (HPX) domain have been shown to modulate activity of the MMP-1 catalytic (CAT) domain. In order to reveal the structural and conformational effects of such mutations, a molecular dynamics (MD) study was performed of in silico mutated residues in the X-ray crystallographic structure of MMP-1 complexed with a collagen-model triple-helical peptide (THP). The results indicate an important role of the mutated residues in MMP-1 interactions with the THP and communication between the CAT and the HPX domains. Each mutation has a distinct impact on the correlated motions in the MMP-1•THP. An increased collagenase activity corresponded to the appearance of a unique anti-correlated motion and decreased correlated motions, while decreased collagenase activity corresponded both to increased and decreased anti-correlated motions.

Structural Insights from Molecular Dynamics Simulations of Tryptophan 7-Halogenase and Tryptophan 5-Halogenase

Many natural organic compounds with pharmaceutical applications, including antibiotics (chlortetracycline and vancomycin), antifungal compounds (pyrrolnitrin), and chemotherapeutics (salinosporamide A and rebeccamycin) are chlorinated. Halogenating enzymes like tryptophan 7-halogenase (PrnA) and tryptophan 5-halogenase (PyrH) perform regioselective halogenation of tryptophan. In this study, the conformational dynamics of two flavin-dependent tryptophan halogenases—PrnA and PyrH—was investigated through molecular dynamics simulations, which are in agreement with the crystallographic and kinetic experimental studies of both enzymes and provide further explanation of the experimental data at an atomistic level of accuracy. They show that the binding sites of the cofactor-flavin adenine dinucleotide and the substrate do not come into close proximity during the simulations, thus supporting an enzymatic mechanism without a direct contact between them. Two catalytically important active site residues, glutamate (E346/E354) and lysine (K79/K75) in PrnA and PyrH, respectively, were found to play a key role in positioning the proposed chlorinating agent, hypochlorous acid. The changes in the regioselectivity between PrnA and PyrH arise as a consequence of differences in the orientation of substrate in its binding site.

Targeting cancer-specific glycans by cyclic peptide lectinomimics

The transformation from normal to malignant phenotype in human cancers is associated with aberrant cell-surface glycosylation. Thus, targeting glycosylation changes in cancer is likely to provide not only better insight into the roles of carbohydrates in biological systems, but also facilitate the development of new molecular probes for bioanalytical and biomedical applications. In the reported study, we have synthesized lectinomimics based on odorranalectin 1; the smallest lectin-like cyclic peptide isolated from the frog Odorrana grahami skin, and assessed the ability of these peptides to bind specific carbohydrates on molecular and cellular levels. In addition, we have shown that the disulfide bond found in 1 can be replaced with a lactam bridge. However, the orientation of the lactam bridge, peptides 2 and 3, influenced cyclic peptide‘s conformation and thus these peptides’ ability to bind carbohydrates. Naturally occurring 1 and its analog 3 that adopt similar conformation in water bind preferentially l-fucose, and to a lesser degree d-galactose and N-acetyl-d-galactosamine, typically found within the mucin O-glycan core structures. In cell-based assays, peptides 1 and 3 showed a similar binding profile to Aleuria aurantia lectin and these two peptides inhibited the migration of metastatic breast cancer cell lines in a Transwell assay. Altogether, the reported data demonstrate the feasibility of designing lectinomimics based on cyclic peptides.

Dimerization and ligand binding in tyrosylprotein sulfotransferase-2 are influenced by molecular motions

Tyrosylprotein sulfotransferase-2 catalyses important, but a less explored posttranslational modification of proteins. The crystallographic structure of the full complex, containing the enzyme, PAPS cofactor and the peptide substrate (WTFC) revealed important structural and atomistic details about the overall structure, binding sites and interactions in TPST-2, but is not informative about the conformational flexibility, which is a fundamental protein property. In a recent computational study we analyzed the impact of conformational flexibility on key geometric determinants and interactions in WTFC. However, there are still unexplained effects about how the conformational dynamics influences the formation of the catalytically active dimer form of the enzyme and also how the binding of the cofactor and the substrate to the apoenzyme influences the structure and the pattern of correlated motions in the enzyme. In order to provide the missing knowledge we performed a molecular dynamics study on the dimerization mutant W113A, the apoenzyme (APO), enzyme–cofactor complex (ES), and enzyme–substrate complex (ES) and compared the results to the WTFC. We identified new hydrophobic interactions important for the stabilization of the enzyme dimer and provided structural insight about the role of the key mutation W113A for the dimerization in the absence of experimentally-derived structure. In addition we showed that the binding of the substrate and cofactor to the apoenzyme contributes to the stability of the whole active complex, influences the local interactions in the binding site and importantly, affects the pattern of the correlated motions in the entire molecule.