Search

Carboxypeptidase is a Zn-dependent protease that specifically recognises and hydrolyses peptides with a hydrophobic side chain at the C-terminal residue. According to hydrolysis mechanisms proposed in the literature, catalysis requires a water molecule to be close to the Zn ion so as to be activated as a nucleophile. Among small molecules that resemble the slowly hydrolysed Gly-Tyr peptide, which have been previously designed as inhibitors and characterised structurally, a variant with the terminal amino acid in a D-configuration has been the most effective. Our molecular dynamics simulations of carboxypeptidase complexed with different variants of those inhibitor ligands as well as variants of the Gly-Tyr peptide show that the strength of the inhibitory effect is not related to the binding strength of the ligand. Our data rather support an earlier notion that the inhibition is, at least partially, due to blocking a coordination site at the Zn ion by the ligand coordinating the metal ion in a bidentate fashion. [ABSTRACT FROM AUTHOR]

Supramolecular recognition of nucleobases and short sequences is an emerging research field focusing on possible applications to treat many diseases. Controlling the affinity and selectivity of synthetic receptors to target desired nucleotides or short sequences is a highly challenging task. Herein, we elucidate the effect of substituents in the phenyl ring of the anthracene‐benzene azacyclophane on the recognition of nucleoside triphosphates (NTPs) and double‐stranded DNA. We show that introducing phenyl rings increases the affinity for NTPs 10‐fold and implements groove and intercalation binding modes with double‐stranded DNA. NMR studies and molecular modeling calculations support the ability of cyclophanes to encapsulate nucleobases as part of nucleotides. [ABSTRACT FROM AUTHOR]

Pro-drugs, which ideally release their active compound only at the site of action, i.e., in a cancer cell, are a promising approach towards an increased specificity and hence reduced side effects in chemotherapy. A popular form of pro-drugs is esters, which are activated upon their hydrolysis. Since carboxylesterases that catalyse such a hydrolysis reaction are also abundant in normal tissue, it is of great interest whether a putative pro-drug is a probable substrate of such an enzyme and hence bears the danger of being activated not just in the target environment, i.e., in cancer cells. In this work, we study the binding mode of carboxylesters of the drug molecule camptothecin, which is an inhibitor of topoisomerase I, of varying size to human carboxylesterase 2 (HCE2) by molecular docking and molecular dynamics simulations. A comparison to irinotecan, known to be a substrate of HCE2, shows that all three pro-drugs analysed in this work can bind to the HCE2 protein, but not in a pose that is well suited for subsequent hydrolysis. Our data suggest, moreover, that for the irinotecan substrate, a reactant-competent pose is stabilised once the initial proton transfer from the putative nucleophile Ser202 to the His431 of the catalytic triad has already occurred. Our simulation work also shows that it is important to go beyond the static models obtained from molecular docking and include the flexibility of enzyme–ligand complexes in solvents and at a finite temperature. Under such conditions, the pro-drugs studied in this work are unlikely to be hydrolysed by the HCE2 enzyme, indicating a low risk of undesired drug release in normal tissue. [ABSTRACT FROM AUTHOR]

Cytochrome c Oxidase (CcO), a membrane protein of the respiratory chain, pumps protons against an electrochemical gradient by using the energy of oxygen reduction to water. The ("chemical") protons required for this reaction and those pumped are taken up via two distinct channels, named D-channel and K-channel, in a step-wise and highly regulated fashion. In the reductive phase of the catalytic cycle, both channels transport protons so that the pumped proton passes the D-channel before the "chemical" proton has crossed the K-channel. By performing molecular dynamics simulations of CcO in the O→E redox state (after the arrival of the first reducing electron) with various combinations of protonation states of the D- and K-channels, we analysed the effect of protonation on the two channels. In agreement with previous work, the amount of water observed in the D-channel was significantly higher when the terminal residue E286 was not (yet) protonated than when the proton arrived at this end of the D-channel and E286 was neutral. Since a sufficient number of water molecules in the channel is necessary for proton transport, this can be understood as E286 facilitating its own protonation. K-channel hydration shows an even higher dependence on the location of the excess proton in the K-channel. Also in agreement with previous work, the K-channel exhibits a very low hydration level that likely hinders proton transfer when the excess proton is located in the lower part of the K-channel, that is, on the N-side of S365. Once the proton has passed S365 (towards the reaction site, the bi-nuclear centre (BNC)), the amount of water in the K-channel provides hydrogen-bond connectivity that renders proton transfer up to Y288 at the BNC feasible. No significant direct effect of the protonation state of one channel on the hydration level, hydrogen-bond connectivity, or interactions between protein residues in the other channel could be observed, rendering proton conductivity in the two channels independent of each other. Regulation of the order of proton uptake and proton passage in the two channels such that the "chemical" proton leaves its channel last must, therefore, be achieved by other means of communication, such as the location of the reducing electron. [ABSTRACT FROM AUTHOR]

Zinc ions are the second most abundant ions found in humans. Their role in proteins can be merely structural but also catalytic, owing to their transition metal character. Modelling their geometric–coordination versatility by empirical force fields is, thus, a challenging task. In this work, we evaluated three popular models, specifically designed to represent zinc ions with regard to their capability of preserving structural integrity. To this end, we performed molecular dynamics simulations of two zinc-containing protein–DNA complexes, which differed in their zinc coordination, i.e., four cysteines or two cysteines and two histidines. The most flexible non-bonded 12-6-4 Lennard–Jones-type model shows a preference for six-fold coordination of the Zn 2 + -ions in contradiction to the crystal structure. The cationic dummy atom model favours tetrahedral geometry, whereas the bonded extended zinc AMBER force field model, by construction, best preserves the initial geometry of a regular or slightly distorted tetrahedron. Our data renders the extended zinc AMBER force field the best model for structural zinc ions in a given geometry. In more complicated cases, though, more flexible models may be advantageous. [ABSTRACT FROM AUTHOR]

Cytochrome c oxidase is a membrane protein of the respiratory chain that consumes protons and molecular oxygen to produce water and uses the resulting energy to pump protons across the membrane. Our molecular dynamics simulations with an excess proton located at different positions in one of the proton-conducting channels, the K-channel, show a clear dependence of the number of water molecules inside the channel on the proton position. A higher hydration level facilitates the formation of hydrogen-bonded chains along which proton transfer can occur. However, a sufficiently high hydration level for such proton transport is observed only when the excess proton is located above S365, i.e., the lower third of the channel. From the channel entrance up to this point, proton transport is via water molecules as proton carriers. These hydronium ions move with their surrounding water molecules, up to K362, filling and widening the channel. The conformation of K362 depends on its own protonation state and on the hydration level, suggesting its role to be proton transport from a hydronium ion at the height of K362 to the upper part of the channel via a conformational change. The protonation-dependent conformational dynamics of E101 at the bottom of the channel renders proton transfer via E101 unlikely. Instead, its role is rather that of an amplifier of H96's proton affinity, suggesting H96 as the initial proton acceptor. [ABSTRACT FROM AUTHOR]

A mechanism involving two-metal ions and double Lewis-acid activation of the substrate, combined with leaving group stabilization is proposed, for a number of phosphoryltransfer enzymes, including the exonuclease subunit of DNA polymerase I. It is designed as a synthetic phosphoryl transfer catalyst the dicopper(II)macrocyclic complex LCu2 which is inspired by the active site structure of exonuclease subunit of DNA polymerase I enzyme.

The ability to locally tune solute–water interactions and thus control the hydrophilic/hydrophobic character of a solute is key to control molecular self-assembly and to develop new drugs and biocatalysts; it has been a holy grail in synthetic chemistry and biology. To date, the connection between (i) the hydrophobicity of a functional group; (ii) the local structure and thermodynamics of its hydration shell; and (iii) the relative influence of van der Waals (dispersion) and electrostatic interactions on hydration remains unclear. We investigate this connection using spectroscopic, classical simulation and ab initio methods by following the transition from hydrophile to hydrophobe induced by the step-wise fluorination of methyl groups. Along the transition, we find that water–solute hydrogen bonds are progressively transformed into dangling hydroxy groups. Each structure has a distinct thermodynamic, spectroscopic and quantum-mechanical signature connected to the associated local solute hydrophobicity and correlating with the relative contribution of electrostatics and dispersion to the solute–water interactions. [ABSTRACT FROM AUTHOR]

The S1 state geometries of benzonitrile, p-cyanophenol, o-cyanophenol, chlorobenzene, and p-chlorophenol were determined by Franck-Condon simulations and a fit of the geometry to the vibronic intensities and effective rotational constants in the harmonic limit based on ab initio force constants. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Molecular dynamics (MD) simulations in both crystalline and aqueous solution environments is utilized to examine the Diels-Alderase ribozyme, an in vitro-evolved ribonucleic acid enzyme which accelerates the formation of carbon-carbon bonds between an anthracene diene and a maleimide dienophile in a [4 + 2] cycloaddition. The simulation results obtained for the reaction with broad application in organic chemistry demonstrated that full catalytic activity of the reaction depends on the [Mg.sup.2+] ion concentration.

The catalytic mechanism of backbone phosphodiester hydrolysis by the restriction enzyme EcoRV is examined by using hybrid quantum mechanical/molecular mechanical calculations. The protein environment is tuned to electrostatically stabilize the transition state structures and the direct catalytic impact of essential residues is analyzed.

The transfer of protons through proton translocating channels is a complex process, for which direct samplings of different protonation states and side chain conformations in a transition network calculation provide an efficient, bias-free description. In principle, a new transition network calculation is required for every unsampled change in the system of interest, e.g. an unsampled protonation state change, which is associated with significant computational costs. Transition networks void of or including an unsampled change are termed unperturbed or perturbed, respectively. Here, we present a prediction method, which is based on an extensive coarse-graining of the underlying transition networks to speed up the calculations. It uses the minimum spanning tree and a corresponding sensitivity analysis of an unperturbed transition network as initial guess and refinement parameter for the determination of an unknown, perturbed transition network. Thereby, the minimum spanning tree defines a sub-network connecting all nodes without cycles and minimal edge weight sum, while the sensitivity analysis analyzes the stability of the minimum spanning tree towards individual edge weight reductions. Using the prediction method, we are able to reduce the calculation costs in a model system by up to 80%, while important network properties are maintained in most predictions. [ABSTRACT FROM AUTHOR]

Abstract: Apurinic/apyrimidinic endonuclease 1 (Ape1) is an important metal‐dependent enzyme in the base excision repair mechanism, responsible for the backbone cleavage of abasic DNA through a phosphate hydrolysis reaction. Molecular dynamics simulations of Ape1 complexed to its substrate DNA performed for models containing 1 or 2 Mg2+‐ions as cofactor located at different positions show a complex with 1 metal ion bound on the leaving group site of the scissile phosphate to be the most likely reaction‐competent conformation. Active‐site residue His309 is found to be protonated based on pKa calculations and the higher conformational stability of the Ape1‐DNA substrate complex compared to scenarios with neutral His309. Simulations of the D210N mutant further support the prevalence of protonated His309 and strongly suggest Asp210 as the general base for proton acceptance by a nucleophilic water molecule. [ABSTRACT FROM AUTHOR]

Distortions in the DNA sequence, such as damage or mispairs, are specifically recognized and processed by DNA repair enzymes. Many repair proteins and, in particular, glycosylases flip the target base out of the DNA helix into the enzyme's active site. Our molecular dynamics simulations of DNA with intact and damaged (oxidized) methyl-cytosine show that the probability of being flipped is similar for damaged and intact methyl-cytosine. However, the accessibility of the different 5-methyl groups allows direct discrimination of the oxidized forms. Hydrogen-bonded patterns that vary between methyl-cytosine forms carrying a carbonyl oxygen atom are likely to be detected by the repair enzymes and may thus help target site recognition. [ABSTRACT FROM AUTHOR]

Distortions in the DNA sequence such as damages or mispairs are specifically recognized and processed by DNA repair enzymes. A particular challenge for the enzymatic specificity is the recognition of a wrongly-placed native nucleotide such as thymine in T:G mispairs. An important step of substrate binding which is observed in many repair proteins is the flipping of the target base out of the DNA helix into the enzyme's active site. In this work we investigate how much the intrinsic dynamics of mispaired DNA is changed compared to canonical DNA. Our molecular dynamics simulations of DNA with and without T:G mispairs show significant differences in the conformation of paired and mispaired DNA. The wobble pair T:G shows local distortions such as twist, shear and stretch which deviate from canonical B form values. Moreover, the T:G mispair is found to be kinetically less stable, exhibiting two states with respect to base opening: a closed state comparable to the canonical base pairs, and a more open state, indicating a proneness for base flip. In addition, we observe that the thymine base in a T:G mispair is significantly more probable to be flipped than thymine in a T:A pair or cytosine in a C:G pair. Such local deformations and in particular the existence of a second, more-open state can be speculated to help the target-site recognition by repair enzymes. [ABSTRACT FROM AUTHOR]

Understanding how the chemical environment modulatesthe predominantconformations and kinetics of flexible molecules is a core interestof biochemistry and a prerequisite for the rational design of syntheticcatalysts. This study combines molecular dynamics simulation and Markovstate models (MSMs) to a systematic computational strategy for investigatingthe effect of the chemical environment of a molecule on its conformationsand kinetics. MSMs allow quantities to be computed that are otherwisedifficult to access, such as the metastable sets, their free energies,and the relaxation time scales related to the rare transitions betweenmetastable states. Additionally, MSMs are useful to identify observablesthat may act as sensors for the conformational or binding state ofthe molecule, thus guiding the design of experiments. In the presentstudy, the conformation dynamics of UDP-GlcNAc are studied in vacuum,water, water + Mg2+, and in the protein UDP-GlcNAc 2-epimerase.It is found that addition of Mg2+significantly affectsthe conformational stability, thermodynamics, and kinetics of UDP-GlcNAc.In particular, the slowest structural process, puckering of the GlcNAcsugar, depends on the overall conformation of UDP-GlcNAc and may thusact as a sensor of whether Mg2+is bound or not. Interestingly,transferring the molecule from vacuum to water makes the protein-bindingconformations UDP-GlcNAc first accessible, while adding Mg2+further stabilizes them by specifically associating to binding-competentconformations. While Mg2+is not cocrystallized in theUDP-GlcNAc 2-epimerase complex, the selectively stabilized Mg2+/UDP-GlcNAc complex may be a template for the bound state,and Mg2+may accompany the binding-competent ligand conformationto the binding pocket. This serves as a possible explanation of theenhanced epimerization rate in the presence of Mg2+. Thisrole of Mg2+has previously not been described and opensthe question whether “binding co-factors” may be a conceptof general relevance for protein–ligand binding. [ABSTRACT FROM AUTHOR]

The Diels-Alderase ribozyme is an in vitro-evolved ribonucleic acid enzyme that catalyzes a [4 + 2] cycloaddition reaction between an anthracene diene and a maleimide dienophile. The ribozyme can in principle be used to selectively synthesize only one product enantiomer, depending on which of the two entrances to the catalytic pocket, 'front' or 'back', the substrate is permitted to use. Here, we investigate stereoselection and substrate recognition in the ribozyme by means of multiple molecular dynamics simulations, performed on each of the two substrates individually in the pocket, on the reactant state, and on the product state. The results are consistent with a binding mechanism in which the maleimide likely binds first followed by the anthracene, which enters preferentially through the front door. The free energy profiles for anthracene binding indicate that the pre-( R, R)-enantiomer conformation is slightly preferred, in agreement with the experimentally observed small enantiomeric excess of the ( R, R)-enantiomer of the product. The reactant state is stabilized by the simultaneous presence of both substrates bound to their binding sites in the hydrophobic pocket as well as by stacking interactions between them. © 2012 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]