The role of protein-bound water molecules in microbial rhodopsins

Protein-bound internal water molecules are essential features of the structure and function of microbial rhodopsins. Besides structural stabilization, they act as proton conductors and even proton storage sites. Currently, the most understood model system exhibiting such features is bacteriorhodopsin (bR). During the last 20 years, the importance of water molecules for proton transport has been revealed through this protein. It has been shown that water molecules are as essential as amino acids for proton transport and biological function. In this review, we present an overview of the historical development of this research on bR. We furthermore summarize the recently discovered protein-bound water features associated with proton transport. Specifically, we discuss a pentameric water/amino acid arrangement close to the protonated Schiff base as central proton-binding site, a protonated water cluster as proton storage site at the proton-release site, and a transient linear water chain at the proton uptake site. We highlight how protein conformational changes reposition or reorient internal water molecules, thereby guiding proton transport. Last, we compare the water positions in bR with those in other microbial rhodopsins to elucidate how protein-bound water molecules guide the function of microbial rhodopsins. This article is part of a Special Issue entitled: Retinal Proteins — You can teach an old dog new tricks.     

His74 conservation in the bilin reductase PcyA family reflects an important role in protein-substrate structure and dynamics

Phycocyanobilin:ferredoxin oxidoreductase (PcyA) catalyzes the proton-coupled four-electron reduction of biliverdin IXα’s two vinyl groups to produce phycocyanobilin, an essential chromophore for phytochromes, cyanobacteriochromes and phycobiliproteins. Previous site directed mutagenesis studies indicated that the fully conserved residue His74 plays a critical role in the H-bonding network that permits proton transfer. Here, we exploit X-ray crystallography, enzymology and molecular dynamics simulations to understand the functional role of this invariant histidine. The structures of the H74A, H74E and H74Q variants of PcyA reveal that a “conserved” buried water molecule that bridges His74 and catalytically essential His88 is not required for activity. Despite distinct conformations of Glu74 and Gln74 in the H74E and H74Q variants, both retain reasonable activity while the H74A variant is inactive, suggesting smaller residues may generate cavities that increase flexibility, thereby reducing enzymatic activity. Molecular dynamic simulations further reveal that the crucial active site residue Asp105 is more dynamic in H74A compared to wild-type PcyA and the two other His74 variants, supporting the conclusion that the Ala74 mutation has increased the flexibility of the active site.     

Determination of Lennard-Jones interaction parameters using a new procedure

For theoretical and chemical engineering applications, accurate and, if possible, simple models of molecular interactions are needed. We have recently proposed a new procedure for determining Lennard-Jones interaction parameters for fluids, forcing agreement between the values of the pressure obtained from empirical equations of state and those obtained from computer simulations. In this work we obtain new intermolecular Lennard-Jones parameters for non-polar molecules, taking into account their deviation from the spherical shape by means of an acentric factor. Our procedure could help to connect the microscopic and macroscopic worlds and it will be progressively implemented in order to obtain a better representation of other substances and mixtures of chemical interest.     

Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors

Since the emergence of the severe acute respiratory syndrome (SARS) to date, neither an effective antiviral drug nor a vaccine against SARS is available. However, it was found that a mixture of two HIV-1 proteinase inhibitors, lopinavir and ritonavir, exhibited some signs of effectiveness against the SARS virus. To understand the fine details of the molecular interactions between these proteinase inhibitors and the SARS virus via complexation, molecular dynamics simulations were carried out for the SARS-CoV 3CL^pro free enzyme (free SARS) and its complexes with lopinavir (SARS-LPV) and ritonavir (SARS-RTV). The results show that flap closing was clearly observed when the inhibitors bind to the active site of SARS-CoV 3CL^pro. The binding affinities of LPV and RTV to SARS-CoV 3CL^pro do not show any significant difference. In addition, six hydrogen bonds were detected in the SARS-LPV system, while seven hydrogen bonds were found in SARS-RTV complex.     

Computational analysis on conformational dynamics of bone morphogenetic protein-2 (BMP-2)

BMP-2 is widely used for bone regeneration because of its ability to induce osteoblast differentiation and proliferation. The pharmaceutical application of BMP-2 as bone implant makes the studies on stability and conformational dynamics very relevant as proteins are functional only in their native three-dimensional state. Knowing the factors affecting BMP-2 structure becomes essential for designing bone implants activated by BMP-2. In order to explore the influence of temperature and hydration on protein conformation, we have performed the molecular dynamics (MD) simulations at the time scale of 100 ns with two different force fields. We have examined the dynamic behaviour of BMP-2 monomer and dimer in aqueous medium as well as in vacuum at four different temperatures (300, 350, 400 and 450 K). MD simulation of BMP-2 monomer and dimer in water and vacuum environments shows the major contribution of water in structure stabilization. Temperature of the system affects the secondary structure differently in case of monomer and dimer simulation and the dynamics also depends on the environment viz. vacuum and aqueous. Vacuum simulations show very early loss of the major secondary structure content. On the other hand, BMP-2 monomer and dimer in aqueous environment show the unfolding of α-helix with increasing temperature. This unfolded α-helix is converted into β-sheet at 400 K in monomer of BMP-2. Contrary to this, we did not observe β-sheet formation in dimer BMP-2 even at 450 K indicating that monomers are more aggregation prone entity as compared to dimers of BMP-2.     

Molecular mechanism of the association and dissociation of Deltarasin from the heterodimeric KRas4B-PDEδ complex

The formation of the KRas4B-PDEδ complex activates different signaling pathways required for the development and maintenance of cancer. Previous experimental and theoretical studies have allowed researchers to design an inhibitor of the KRas4B-PDEδ complex, “Deltarasin.” This inhibitor binds to the prenyl-binding pocket of PDEδ and subsequently inhibits the proliferation of human pancreatic ductal adenocarcinoma cells that depend on oncogenic KRas4B. Nevertheless, structural and energetic information about the inhibitory effects of Deltarasin on the KRas4B-PDEδ complex are not available. In this study, we explore the properties of Deltarasin in inhibiting the formation of wild-type and mutant KRas4B-PDEδ complexes present in different cell lines expressing mutant RAS genes (G12D, G12C, G12V, G13D, Q61L, and Q61R) using 1.7 μs molecular dynamics (MD) simulations in combination with the MMGBSA approach. Our results revealed the energetic and structural mechanisms that suggest a higher affinity of Deltarasin for PDEδ than the farnesylated HVR. Moreover, Deltarasin exerts another dissociative effect by binding to the protein-protein dimeric interface of wild-type KRas4B-PDEδ, whereas associative and dissociative effects were observed for mutant KRas4B-PDEδ, providing a mechanistic explanation for the inhibitory effects of Deltarasin on different cancer cell lines.