CAVER Analyst 1.0: graphic tool for interactive visualization and analysis of tunnels and channels in protein structures

ABSTRACT Summary: The transport of ligands, ions or solvent molecules into proteins with buried binding sites or through the membrane is enabled by protein tunnels and channels. CAVER Analyst is a software tool for calculation, analysis and real-time visualization of access tunnels and channels in static and dynamic protein structures. It provides an intuitive graphic user interface for setting up the calculation and interactive exploration of identified tunnels/channels and their characteristics. Availability and Implementation: CAVER Analyst is a multi-platform software written in JAVA. Binaries and documentation are freely available for non-commercial use at http://www.caver.cz

A Five Residue Insertion Between Codons 28 And 29 Of The Hiv-1 Protease Gene Reduces The Replicative Capacity Of The Virus

HIV-1 protease (PR) is a 99 amino acid protein responsible for cleavage of the viral polyprotein. We have identified a novel clinical isolate, MDR/28, which contains a five residue insertion between codons 28 and 29 of a multi-drug resistant (MDR) PR. This clinical isolate displays reduced viral replicative capacity compared to the wild-type. As opposed to drug-resistance mutations, studies on insertions remain largely underrepresented in the literature, and the consequences of such insertions are largely unknown. To understand the mechanism leading to reduced replicative capacity, three PR models were created and subjected to 40ns molecular dynamics simulations: MDR/28, wild type, and MDR PR. In addition, PR inhibitors (PI) atazanavir (ATV), darunavir (DRV), lopinavir (LPV) and saquinavir (SQV), as well as cleavage peptide CA/p2 were docked to the three models. The MDR/28-PI complexes displayed decreased binding affinity when compared to WT complexes, likely due to an increased active site cavity volume and altered secondary structure at residues local to the insertion mutant. Additionally, in the active site of MDR/28 the predicted binding mode of the CA/p2 peptide did not include contact with the catalytic residues, and migrated from that position, a behavior not seen with any tested PIs or with either of the other PR models. These structural changes produced by the insertion suggest a mechanism for reduced replicative capacity of the mutant virus

Regulating effect of β-ketoacyl synthase domain of fatty acid synthase on fatty acyl chain length in de novo fatty acid synthesis

Fatty acid synthase (FAS) is a multifunctional homodimeric protein, and is the key enzyme required for the anabolic conversion of dietary carbohydrates to fatty acids. FAS synthesizes long-chain fatty acids from three substrates: acetyl-CoA as a primer, malonyl-CoA as a 2 carbon donor, and NADPH for reduction. The entire reaction is composed of numerous sequential steps, each catalyzed by a specific functional domain of the enzyme. FAS comprises seven different functional domains, among which the β-ketoacyl synthase (KS) domain carries out the key condensation reaction to elongate the length of fatty acid chain. Acyl tail length controlled fatty acid synthesis in eukaryotes is a classic example of how a chain building multienzyme works. Different hypotheses have been put forward to explain how those sub-units of FAS are orchestrated to produce fatty acids with proper molecular weight. In the present study, molecular dynamics simulation based binding free energy calculation and access tunnels analysis showed that the C16 acyl tail fatty acid, the major product of FAS, fits to the active site on KS domain better than any other substrates. These simulations supported a new hypothesis about the mechanism of fatty acid production ratio: the geometric shape of active site on KS domain might play a determinate role

Visual analysis of protein-ligand interactions

The analysis of protein-ligand interactions is complex because of the many factors at play. Most current methods for visual analysis provide this information in the form of simple 2D plots, which, besides being quite space hungry, often encode a low number of different properties. In this paper we present a system for compact 2D visualization of molecular simulations. It purposely omits most spatial information and presents physical information associated to single molecular components and their pairwise interactions through a set of 2D InfoVis tools with coordinated views, suitable interaction, and focus+context techniques to analyze large amounts of data. The system provides a wide range of motifs for elements such as protein secondary structures or hydrogen bond networks, and a set of tools for their interactive inspection, both for a single simulation and for comparing two different simulations. As a result, the analysis of protein-ligand interactions of Molecular Simulation trajectories is greatly facilitated.Peer ReviewedPostprint (author's final draft

Modeling Kinetics and Thermodynamics of Guest Encapsulation into the [ML] 12- Supramolecular Organometallic Cage

Altres ajuts: Acord transformatiu CRUE-CSICThe encapsulation of molecular guests into supramolecular hosts is a complex molecular recognition process in which the guest displaces the solvent from the host cavity, while the host deforms to let the guest in. An atomistic description of the association would provide valuable insights on the physicochemical properties that guide it. This understanding may be used to design novel host assemblies with improved properties (i.e., affinities) toward a given class of guests. Molecular simulations may be conveniently used to model the association processes. It is thus of interest to establish efficient protocols to trace the encapsulation process and to predict the associated magnitudes Δ G and Δ G ⧧. Here, we report the calculation of the Gibbs energy barrier and Gibbs binding energy by means of explicit solvent molecular simulations for the [GaL] 12- metallocage encapsulating a series of cationic molecules. The Δ G ⧧ for encapsulation was estimated by means of umbrella sampling simulations. The steps involved were identified, including ion-pair formation and naphthalene rotation (from L ligands of the metallocage) during the guest's entrance. The Δ G values were computed using the attach-pull-release method. The results reveal the sensitivity of the estimates on the force field parameters, in particular on atomic charges, showing that higher accuracy is obtained when charges are derived from implicit solvent quantum chemical calculations. Correlation analysis identified some indicators for the binding affinity trends. All computed magnitudes are in very good agreement with experimental observations. This work provides, on one side, a benchmarked way to computationally model a highly charged metallocage encapsulation process. This includes a nonstandard parameterization and charge derivation procedure. On the other hand, it gives specific mechanistic information on the binding processes of [GaL] 12- at the molecular level where key motions are depicted. Taken together, the study provides an interesting option for the future design of metal-organic cages

Gas channel rerouting in a primordial enzyme: Structural insights of the carbon-monoxide dehydrogenase/acetyl-CoA synthase complex from the acetogen Clostridium autoethanogenum

Clostridium autoethanogenum, the bacterial model for biological conversion of waste gases into biofuels, grows under extreme carbon-monoxide (CO) concentrations. The strictly anaerobic bacterium derives its entire cellular energy and carbon from this poisonous gas, therefore requiring efficient molecular machineries for CO-conversion. Here, we structurally and biochemically characterized the key enzyme of the CO-converting metabolism: the CO-dehydrogenase/Acetyl-CoA synthase (CODH/ACS). We obtained crystal structures of natively isolated complexes from fructose-grown and CO-grown C. autoethanogenum cultures. Both contain the same isoforms and if the overall structure adopts the classic alpha(2)beta(2) architecture, comparable to the model enzyme from Moorella thermoacetica, the ACS binds a different position on the CODH core. The structural characterization of a proteolyzed complex and the conservation of the binding interface in close homologs rejected the possibility of a crystallization artefact. Therefore, the internal CO-channeling system, critical to transfer CO generated at the C-cluster to the ACS active site, drastically differs in the complex from C. autoethanogenum. The 1.9-angstrom structure of the CODH alone provides an accurate picture of the new CO-routes, leading to the ACS core and reaching the surface. Increased gas accessibility would allow the simultaneous CO-oxidation and acetyl-CoA production. Biochemical experiments showed higher flexibility of the ACS subunit from C. autoethanogenum compared to M. thermoacetica, albeit monitoring similar CO-oxidation and formation rates. These results show a reshuffling of internal CO-tunnels during evolution of these Firmicutes, putatively leading to a bidirectional complex that ensure a high flux of CO-conversion toward energy conservation, acting as the main cellular powerplant

Single active-site mutants are sufficient to enhance serine:pyruvate α-transaminase activity in an ω-transaminase

We have analysed the natural evolution of transaminase structure and sequence between an α-transaminase serine-pyruvate aminotransferase, and an ω-transaminase from Chromobacterium violaceum with <20% sequence identity, and identified the active-site regions which are least conserved structurally. We also show that these structural changes correlate strongly with transaminase substrate specificity during evolution and therefore might normally be presumed to be essential determinants of substrate specificity. However, key residues are often conserved spatially during evolution and yet come from within a different region of the sequence via structural reorganisations. Here we also show that α-transaminase-type serine-pyruvate aminotransferase activity, can be engineered into the CV2025 ω-transaminase scaffold with any one of many possible single point mutations at three key positions, without the requirement for significant backbone remodeling, or repositioning of the residue from a different region of sequence. This finding has significant implications for enzyme redesign in which solutions to substrate specificity changes may be found that are significantly more efficient than by engineering in all sequence and structure determinants identified by correlation to substrate specificity. This article is protected by copyright. All rights reserved

New Thermophilic α/β Class Epoxide Hydrolases Found in Metagenomes From Hot Environments

This is the final version. Available from Frontiers Media via the DOI in this record.Two novel epoxide hydrolases (EHs), Sibe-EH and CH65-EH, were identified in the metagenomes of samples collected in hot springs in Russia and China, respectively. The two α/β hydrolase superfamily fold enzymes were cloned, over-expressed in Escherichia coli, purified and characterized. The new EHs were active toward a broad range of substrates, and in particular, Sibe-EH was excellent in the desymmetrization of cis-2,3-epoxybutane producing the (2R,3R)-diol product with ee exceeding 99%. Interestingly these enzymes also hydrolyse (4R)-limonene-1,2-epoxide with Sibe-EH being specific for the trans isomer. The Sibe-EH is a monomer in solution whereas the CH65-EH is a dimer. Both enzymes showed high melting temperatures with the CH65-EH being the highest at 85°C retaining 80% of its initial activity after 3 h thermal treatment at 70°C making it the most thermal tolerant wild type epoxide hydrolase described. The Sibe-EH and CH65-EH have been crystallized and their structures determined to high resolution, 1.6 and 1.4 Å, respectively. The CH65-EH enzyme forms a dimer via its cap domains with different relative orientation of the monomers compared to previously described EHs. The entrance to the active site cavity is located in a different position in CH65-EH and Sibe-EH in relation to other known bacterial and mammalian EHs