Conformational analysis of retinoids and restriction of their dynamics by retinoid-binding proteins

Vitamin A and its natural derivatives

Vitamin A and derivatives, the natural retinoids, underpin signaling pathways of cellular differentiation, and are key chromophores in vision. These functions depend on transfer across membranes, and carrier proteins to shuttle retinoids to specific cell compartments. Natural retinoids, ultimately derived from plant carotenoids by metabolism to all-trans retinol, are lipophilic and consist of a cyclohexenyl (β-ionone) moiety linked to a polyene chain. This structure constrains the orientation of retinoids within lipid membranes. Cis-trans isomerization at double bonds of the polyene chain and s-cis/s-trans rotational isomerization at single bonds define the functional dichotomy of retinoids (signaling/vision) and specificities of interactions with specific carrier proteins and receptors. Metabolism of all-trans retinol to 11-cis retinal, transfer to photoreceptors, and removal and recycling of all-trans retinal generated by photoreceptor irradiation, is the key process underlying vision. All-trans retinol transferred into cells is metabolized to all-trans retinoic acid and shuttled to the cell nucleus to regulate gene expression controlling organ, tissue and cell differentiation, and cellular homeostasis. Research methods need to address the potential of photoisomerization in vitro to confound research results, and data should be interpreted in the context of membrane-association properties of retinoids and physiological concentrations in vivo. Despite a century of research, there are many fundamental questions of retinoid cellular biochemistry and molecular biology still to be answered. Computational modeling techniques will have an important role for understanding the nuances of vitamin A signaling and function.

Simulated (Un)Binding of Arachidonic Acid in the Cyclooxygenase site of Prostaglandin H2 Synthase-1

Molecular dynamics simulations with external forces are employed to study the unbinding and binding of arachidonic acid (AA) in the cyclooxygenase (COX) site of prostaglandin H2 synthase-1. Simulations with AA inside the COX binding channel reveal sequences of concerted bond rotations in the fatty acid alkyl chain which obviate the need for gross conformational changes in the protein and substrate during unbinding and binding. The all-cis structure of AA, with double bonds separated by two single bonds, facilitates easy access to the COX channel and correct positioning inside the active site for the COX chemistry to occur. Two derivatives of AA, one with a cis double bond changed to a trans configuration and the other with a double bond reduced to a single bond, are also studied. In both cases the concertedness of bond rotations in the fatty acid chain is diminished and larger forces are required to move the fatty acid inside the COX channel. Important motions of residues near the mouth of the COX channel are found and analyzed. In particular, a conformational “switch” involving Arg83, Glu524 and Arg120 is seen to mediate the movement of the substrate from the membrane to the channel.

The molecular basis of the interactions between synthetic retinoic acid analogues and the retinoic acid receptors

All-trans-retinoic acid (ATRA) and its synthetic analogues EC23 and EC19 direct cellular differentiation by interacting as ligands for the retinoic acid receptor (RARα, β and γ) family of nuclear receptor proteins. To date, a number of crystal structures of natural and synthetic ligands complexed to their target proteins have been solved, providing molecular level snap-shots of ligand binding. However, a deeper understanding of receptor and ligand flexibility and conformational freedom is required to develop stable and effective ATRA analogues for clinical use. Therefore, we have used molecular modelling techniques to define RAR interactions with ATRA and two synthetic analogues, EC19 and EC23, and compared their predicted biochemical activities to experimental measurements of relative ligand affinity and recruitment of coactivator proteins. A comprehensive molecular docking approach that explored the conformational space of the ligands indicated that ATRA is able to bind the three RAR proteins in a number of conformations with one extended structure being favoured. In contrast the biologically-distinct isomer, 9-cis-retinoic acid (; 9CRA), showed significantly less conformational flexibility in the RAR binding pockets. These findings were used to inform docking studies of the synthetic retinoids EC23 and EC19, and their respective methyl esters. EC23 was found to be an excellent mimic for ATRA, and occupied similar binding modes to ATRA in all three target RAR proteins. In comparison, EC19 exhibited an alternative binding mode which reduces the strength of key polar interactions in RARα/γ but is well-suited to the larger RARβ binding pocket. In contrast, docking of the corresponding esters revealed the loss of key polar interactions which may explain the much reduced biological activity. Our computational results were complemented using an in vitro binding assay based on FRET measurements, which showed that EC23 was a strongly binding, pan-agonist of the RARs, while EC19 exhibited specificity for RARβ, as predicted by the docking studies. These findings can account for the distinct behaviour of EC23 and EC19 in cellular differentiation assays, and additionally, the methods described herein can be further applied to the understanding of the molecular basis for the selectivity of different retinoids to RARα, β and γ.

Conformational Analysis of Retinoic Acids: Effects of Steric Interactions on Nonplanar Conjugated Polyenes

Retinoic acids and other vitamin A analogs contain a trimethylcyclohexenyl ring in conjugation with a polyene chain joined at carbon-6 (C6) and carbon-7 (C7). A MP2-SCS/cc-pVDZ// B3LYP/6-31G(d) 2-D potential energy surface was computed for all-trans retinoic acid, which had 6 minima (3 enantiomeric pairs). The global minima were distorted s-gauche enantiomers (_6-7 = 53°) with half-chair conformations of the ring. Distorted s-gauche enantiomers (_6-7 = 55°) with inverted half-chair ring conformations were 1.7 kJ/mol above the global minima. The s-trans enantiomers (_6-7 = 164°) were 11.3 kJ/mol above the global minima. Steric energies were computed by the method of Guo and Karplus to identify key structural elements in retinoic acids which determines their conformation. Small molecule crystal structures in the CCDC database with trimethylcyclohexenyl ring and exocyclic double bonds have ring-chain geometries near to one of the 6 energy minima of retinoic acids, except for retinaldehyde iminium cations.

Conformational Analysis of Free and Bound Retinoic Acid

The conformational profiles of unbound all-trans and 9-cis retinoic acid (RA) have been determined using classical and quantum mechanical calculations. Sixty-six all-trans-RA (ATRA) and forty-eight 9-cis-RA energy minimum conformers were identified via HF/6-31G* geometry optimizations in vacuo. Their relative conformational energies were estimated utilizing the M06, M06-2x and MP2 methods combined with the 6-311+G(d,p), aug-cc-pVDZ and aug-cc-pVTZ basis sets, as well as complete basis set MP2 extrapolations using the latter two basis sets. Single-point energy calculations performed with the M06-2x density functional were found to yield similar results to MP2/CBS for the low-energy retinoic acid conformations. Not unexpectedly, the conformational propensities of retinoic acid were governed by the orientation and arrangement of the torsion angles associated with the polyene tail. We also used previously reported QM/MM X-ray refinement results on four ATRA-protein crystal structures plus one newly refined 9-cis-RA complex (PDB ID 1XDK) in order to investigate the conformational preferences of bound retinoic acid. In the re-refined RA conformers the conjugated double bonds are nearly coplanar, which is consistent with the global minimum identified by the Omega/QM method rather than the corresponding crystallographically determined conformations given in the PDB. Consequently, a 91.3% average reduction of the local strain energy in the gas phase, as well as 92.1% in PCM solvent, was observed using the QM/MM refined structures versus the PDB deposited RA conformations. These results thus demonstrate that our QM/MM X-ray refinement approach can significantly enhance the quality of X-ray crystal structures refined by conventional refinement protocols, thereby providing reliable drug-target structural information for use in structure-based drug discovery applications.

Insights into the catalytic mechanism of cofactor-independent phosphoglycerate mutase from X-ray crystallography, simulated dynamics and molecular modeling

Phosphoglycerate mutases catalyze the isomerization of 2 and 3-phosphoglycerates, and are essential for glucose metabolism in most organisms. Here, we further characterize the 2,3-bisphosphoglycerate-independent phosphoglycerate mutase (iPGM) from Bacillus stearothermophilus by determination of a high-resolution (1.4A) crystal structure of the wild-type enzyme and the crystal structure of its S62A mutant. The mutant structure surprisingly showed the replacement of one of the two catalytically essential manganese ions with a water molecule, offering an additional possible explanation for its lack of catalytic activity. Crystal structures invariably show substrate phosphoglycerate to be entirely buried in a deep cleft between the two iPGM domains. Flexibility analyses were therefore employed to reveal the likely route of substrate access to the catalytic site through an aperture created in the enzyme's surface during certain stages of the catalytic process. Several conserved residues lining this aperture may contribute to orientation of the substrate as it enters. Factors responsible for the retention of glycerate within the phosphoenzyme structure in the proposed mechanism are identified by molecular modeling of the glycerate complex of the phosphoenzyme. Taken together, these results allow for a better understanding of the mechanism of action of iPGMs. Many of the results are relevant to a series of evolutionarily related enzymes. These studies will facilitate the development of iPGM inhibitors which, due to the demonstrated importance of this enzyme in many bacteria, would be of great potential clinical significance.

Functional concerted motions in the bovine serum retinol-binding protein

The large concerted motions in the apo/holo bovine serum retinol-binding protein were studied using molecular dynamics simulation and 'essential dynamics' analysis. Initially, concerted motions were calculated from conformational differences between various crystal structures. The dynamic behaviour of the protein in the configurational space directions, described by these concerted motions, is analysed. This reveals that the large backbone dynamics of the protein is not influenced by the presence of retinol. Study of free retinol dynamics and retinol in the retinol binding site reveals that the protein binds retinol in a favourable conformation, as opposed to what has been previously described for the bovine cellular retinol-binding protein.

Unbinding of retinoic acid from its receptor studied by steered molecular dynamics

Retinoic acid receptor (RAR) is a ligand-dependent transcription factor that regulates the expression of genes involved in cell growth, differentiation, and development. Binding of the retinoic acid hormone to RAR is accompanied by conformational changes in the protein which induce transactivation or transrepression of the target genes. In this paper we present a study of the hormone binding/unbinding process in order to clarify the role of some of the amino acid contacts and identify possible pathways of the all-trans retinoic acid binding/unbinding to/from human retinoic acid receptor (hRAR)-gamma. Three possible pathways were explored using steered molecular dynamics simulations. Unbinding was induced on a time scale of 1 ns by applying external forces to the hormone. The simulations suggest that the hormone may employ one pathway for binding and an alternative "back door" pathway for unbinding.

Essential dynamics of DNA containing a cis.syn cyclobutane thymine dimer lesion

Conformational properties of a UV-damaged DNA decamer containing a cis.syn cyclobutane thymine dimer (PD) have been investigated by molecular dynamics (MD) simulations. Results from MD simulations of the damaged decamer DNA show a kink of approximately 21.7 degrees at the PD damaged site and a disruption of H bonding between the 5'-thymine of the PD and its complementary adenine. However, no extra-helical flipping of the 3'-adenine complementary to the PD was observed. Comparison to two undamaged DNA decamers, one with the same sequence and the other with an AT replacing the TT sequence, indicates that these properties are specific to the damaged DNA. Essential dynamics (ED) derived from the MD trajectories of the three DNAs show that the backbone phosphate between the two adenines complementary to the PD of the damaged DNA has considerably larger mobility than the rest of the molecule and occurs only in the damaged DNA. As observed in the crystal structure of T4 endonuclease V in a complex with the damaged DNA, the interaction of the enzyme with the damaged DNA can lead to bending along the flexible joint and to induction of adenine flipping into an extra-helical position. Such motions may play an important role in damage recognition by repair enzymes.