A Hamiltonian electrostatic coupling scheme for hybrid Car–Parrinello molecular dynamics simulations

The molecular mechanism of the ligand exchange reaction of an antibody against a glutathione-coated gold cluster

The labeling of proteins with heavy atom clusters is of paramount importance in biomedical research, but its detailed molecular mechanism remains unknown. Here we uncover it for the particular case of the anti-influenza N9 neuraminidase NC10 antibody against a glutathione-coated gold cluster by means of ab initio QM/MM calculations. We show that the labeling reaction follows an associative double S_N2-like reaction mechanism, involving a proton transfer, with low activation barriers only if one of the two distinct peptide/peptidic ligands (the one that occupies the side position) is substituted. Positively charged residues in the vicinity of the incoming thiol result in strong interactions between the antibody and the AuMPC, favoring the ligand exchange reaction for suitable protein mutants. These results pave the way for future investigations aimed at engineering biomolecules to increase their reactivity towards a desired gold atom cluster.

Copper(II) and the pathological H50Q α-synuclein mutant: Environment meets genetics

Copper is one of the metals described to bind the Parkinson disease-related protein α-synuclein (aSyn), and to promote its aggregation. Although histidine at position 50 in the aSyn sequence is one of the most studied copper-anchoring sites, its precise role in copper binding and aSyn aggregation is still unclear. Previous studies suggested that this residue does not significantly affect copper-mediated aSyn aggregation. However, our findings showed that the aggregation of the pathological H50Q aSyn mutant is enhanced by copper hints otherwise. Despite the inexistence of a model for aSyn H50Q-copper complexation, we discuss possible mechanisms by which this metal contributes to the misfolding and self-assembly of this particular aSyn mutant. Considering the genetic association of the H50Q mutation with familial forms of Parkinson disease, and the fact that copper homeostasis is deregulated in this disorder, understanding the interplay between both factors will shed light into the molecular and cellular mechanisms triggering the development and spreading of the aSyn pathology.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Photophysics of a copper phenanthroline elucidated by trajectory and wavepacket-based quantum dynamics: a synergetic approach

On-the-fly excited state molecular dynamics is a valuable method for studying non-equilibrium processes in excited states and is beginning to emerge as a mature approach much like its ground state counterparts.

Different QM/MM Approaches To Elucidate Enzymatic Reactions: Case Study on ppGalNAcT2

Hybrid QM/MM computational studies can provide invaluable insight into the mechanisms of enzymatic reactions that can be exploited for rational drug design. Various approaches are available for such studies. However, their strengths and weaknesses may not be immediately apparent. Using the retaining glycosyltransferase ppGalNAcT2 as a case study, we compare different methodologies used to obtain reaction paths and transition state information. Ab Initio MD using CPMD coupled with the String Method is used to derive the minimum free energy reaction path. The geometrical features of the free energy path, especially around the transition state, agree with the minimum potential energy path obtained by the much less computationally expensive Nudged Elastic Band method. The barrier energy, however, differs by 8 kcal/mol. The free energy surface generated by metadynamics provides a rough overview of the reaction and can confirm the physical relevance of optimized paths or provide an initial guess for path optimization methods. Calculations of enzymatic reactions are usually performed at best at the DFT level of theory. A comparison of widely used functionals with high-level DLPNO-CCSD(T)/CBS data on the potential energy profile serves as a validation of the usability of DFT for this type of enzymatic reaction. The M06-2X meta-hybrid functional in particular matches the DLPNO-CCSD(T)/CBS reference extremely well with errors within 1 kcal/mol. However, even pure-GGA functional OPBE provides sufficient accuracy for this type of study.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

A quantum-mechanics molecular-mechanics scheme for extended systems

We introduce and discuss a hybrid quantum-mechanics molecular-mechanics (QM-MM) approach for Car-Parrinello DFT simulations with pseudopotentials and planewaves basis, designed for the treatment of periodic systems. In this implementation the MM atoms are considered as additional QM ions having fractional charges of either sign, which provides conceptual and computational simplicity by exploiting the machinery already existing in planewave codes to deal with electrostatics in periodic boundary conditions. With this strategy, both the QM and MM regions are contained in the same supercell, which determines the periodicity for the whole system. Thus, while this method is not meant to compete with non-periodic QM-MM schemes able to handle extremely large but finite MM regions, it is shown that for periodic systems of a few hundred atoms, our approach provides substantial savings in computational times by treating classically a fraction of the particles. The performance and accuracy of the method is assessed through the study of energetic, structural, and dynamical aspects of the water dimer and of the aqueous bulk phase. Finally, the QM-MM scheme is applied to the computation of the vibrational spectra of water layers adsorbed at the TiO2 anatase (1 0 1) solid-liquid interface. This investigation suggests that the inclusion of a second monolayer of H2O molecules is sufficient to induce on the first adsorbed layer, a vibrational dynamics similar to that taking place in the presence of an aqueous environment. The present QM-MM scheme appears as a very interesting tool to efficiently perform molecular dynamics simulations of complex condensed matter systems, from solutions to nanoconfined fluids to different kind of interfaces.

A polarizable QM/MM approach to the molecular dynamics of amide groups solvated in water

The infrared (IR) spectra of polypeptides are dominated by the so-called amide bands. Because they originate from the strongly polar and polarizable amide groups (AGs) making up the backbone, their spectral positions sensitively depend on the local electric fields. Aiming at accurate computations of these IR spectra by molecular dynamics (MD) simulations, which derive atomic forces from a hybrid quantum and molecular mechanics (QM/MM) Hamiltonian, here we consider the effects of solvation in bulk liquid water on the amide bands of the AG model compound N-methyl-acetamide (NMA). As QM approach to NMA we choose grid-based density functional theory (DFT). For the surrounding MM water, we develop, largely based on computations, a polarizable molecular mechanics (PMM) model potential called GP6P, which features six Gaussian electrostatic sources (one induced dipole, five static partial charge distributions) and, therefore, avoids spurious distortions of the DFT electron density in hybrid DFT/PMM simulations. Bulk liquid GP6P is shown to have favorable properties at the thermodynamic conditions of the parameterization and beyond. Lennard-Jones (LJ) parameters of the DFT fragment NMA are optimized by comparing radial distribution functions in the surrounding GP6P liquid with reference data obtained from a "first-principles" DFT-MD simulation. Finally, IR spectra of NMA in GP6P water are calculated from extended DFT/PMM-MD trajectories, in which the NMA is treated by three different DFT functionals (BP, BLYP, B3LYP). Method-specific frequency scaling factors are derived from DFT-MD simulations of isolated NMA. The DFT/PMM-MD simulations with GP6P and with the optimized LJ parameters then excellently predict the effects of aqueous solvation and deuteration observed in the IR spectra of NMA. As a result, the methods required to accurately compute such spectra by DFT/PMM-MD also for larger peptides in aqueous solution are now at hand.

Enzymatic Cleavage of Glycosidic Bonds: Strategies on How to Set Up and Control a QM/MM Metadynamics Simulation

Carbohydrates play crucial roles in many biological processes, from cell-cell adhesion to chemical signaling. Their complexity and diversity, related to α/β anomeric configuration, ring substituents, and conformational variations, require a diverse set of enzymes for their processing. Among them, glycoside hydrolases (GHs) are responsible for the hydrolysis of one of the strongest bonds in nature: the glycosidic bond. These highly specialized biological catalysts select particular conformations their carbohydrate substrates to enhance catalysis. The evolution of this conformation during the reaction of glycosidic bond cleavage, known as the conformational catalytic itinerary, is of fundamental interest in glycobiology, with impact on inhibitor and drug design. Here we review some of the aspects and the main strategies one needs to take into account when simulating a reaction in a GH enzyme using QM/MM metadynamics. Several specific aspects are highlighted, from the importance of the distortion of the substrate at the Michaelis complex to the variable control during the metadynamics simulation or the analysis of the reaction mechanism and conformational itinerary. The increasing speed of computer power and methodological advances have added a vital tool to the study of GH mechanisms, as shown here and recent reviews. It is hoped that this chapter will serve as a first guide for those attempting to perform a metadynamics simulation of these relevant and fascinating enzymes.

Carnosine and Homocarnosine Degradation Mechanisms by the Human Carnosinase Enzyme CN1: Insights from Multiscale Simulations

The endogenous dipeptide l-carnosine, and its derivative homocarnosine, prevent and reduce several pathologies like amytrophic lateral sclerosis (ALS), Alzheimer's disease, and Parkinson's disease. Their beneficial action is severely hampered because of the hydrolysis by carnosinase enzymes, in particular the human carnosinase, hCN1. This belongs to the metallopeptidase M20 family, where a cocatalytic active site is formed by two Zn(2+) ions, bridged by a hydroxide anion. The protein may exist as a monomer and as a dimer in vivo. Here we used hybrid quantum mechanics/molecular mechanics simulations based on the dimeric apoenzyme's structural information to predict the Michaelis complexes with l-carnosine and its derivative homocarnosine. On the basis of our calculations, we suggest that (i) l-carnosine degradation occurs through a nucleophilic attack of a Zn(2+)-coordinated bridging moiety for both monomer and dimer. This mechanistic hypothesis for hCN1 catalysis differs from previous proposals, while it is in agreement with available experimental data. (ii) The experimentally measured higher affinity of homocarnosine for the enzyme relative to l-carnosine might be explained, at least in part, by more extensive interactions inside the monomeric and dimeric hCN1's active site. (iii) Hydrogen bonds at the binding site, present in the dimer but absent in the monomer, might play a role in the experimentally observed higher activity of the dimeric form. Investigations of the enzymatic reaction are required to establish or disprove this hypothesis. Our results may serve as a basis for the design of potent hCN1 inhibitors.

CPMD/GULP QM/MM interface for modeling periodic solids: Implementation and its application in the study of Y-zeolite supported Rhn clusters

We report here the development of hybrid quantum mechanics/molecular mechanics (QM/MM) interface between the plane-wave density functional theory based CPMD code and the empirical force-field based GULP code for modeling periodic solids and surfaces. The hybrid QM/MM interface is based on the electrostatic coupling between QM and MM regions. The interface is designed for carrying out full relaxation of all the QM and MM atoms during geometry optimizations and molecular dynamics simulations, including the boundary atoms. Both Born-Oppenheimer and Car-Parrinello molecular dynamics schemes are enabled for the QM part during the QM/MM calculations. This interface has the advantage of parallelization of both the programs such that the QM and MM force evaluations can be carried out in parallel to model large systems. The interface program is first validated for total energy conservation and parallel scaling performance is benchmarked. Oxygen vacancy in α-cristobalite is then studied in detail and the results are compared with a fully QM calculation and experimental data. Subsequently, we use our implementation to investigate the structure of rhodium cluster (Rhn ; n = 2 to 6) formed from Rh(C2 H4 )2 complex adsorbed within a cavity of Y-zeolite in a reducible atmosphere of H2 gas. © 2016 Wiley Periodicals, Inc.

Deacylation Mechanism and Kinetics of Acyl-Enzyme Complex of Class C β-Lactamase and Cephalothin

Understanding the molecular details of antibiotic resistance by the bacterial enzymes β-lactamases is vital for the development of novel antibiotics and inhibitors. In this spirit, the detailed mechanism of deacylation of the acyl-enzyme complex formed by cephalothin and class C β-lactamase is investigated here using hybrid quantum-mechanical/molecular-mechanical molecular dynamics methods. The roles of various active-site residues and substrate in the deacylation reaction are elucidated. We identify the base that activates the hydrolyzing water molecule and the residue that protonates the catalytic serine (Ser64). Conformational changes in the active sites and proton transfers that potentiate the efficiency of the deacylation reaction are presented. We have also characterized the oxyanion holes and other H-bonding interactions that stabilize the reaction intermediates. Together with the kinetic and mechanistic details of the acylation reaction, we analyze the complete mechanism and the overall kinetics of the drug hydrolysis. Finally, the apparent rate-determining step in the drug hydrolysis is scrutinized.

Interplay of histidine residues of the Alzheimer's disease Aβ peptide governs its Zn-induced oligomerization

Conformational changes of Aβ peptide result in its transformation from native monomeric state to the toxic soluble dimers, oligomers and insoluble aggregates that are hallmarks of Alzheimer's disease (AD). Interactions of zinc ions with Aβ are mediated by the N-terminal Aβ(1-16) domain and appear to play a key role in AD progression. There is a range of results indicating that these interactions trigger the Aβ plaque formation. We have determined structure and functional characteristics of the metal binding domains derived from several Aβ variants and found that their zinc-induced oligomerization is governed by conformational changes in the minimal zinc binding site 6HDSGYEVHH14. The residue H6 and segment 11EVHH14, which are part of this site are crucial for formation of the two zinc-mediated interaction interfaces in Aβ. These structural determinants can be considered as promising targets for rational design of the AD-modifying drugs aimed at blocking pathological Aβ aggregation.

Combined Electrostatically Embedded Multiconfiguration Molecular Mechanics and Molecular Mechanical Method: Application to Molecular Dynamics Simulation of a Chemical Reaction in Aqueous Solution with Hybrid Density Functional Theory

We here combine the electrostatically embedded multiconfiguration molecular mechanics (EE-MCMM) method for generating global potential energy surfaces in the presence of an electrostatic potential with molecular mechanics (MM). The resulting EE-MCMM/MM method is illustrated by applying it to carry out a molecular dynamics simulation for the symmetric bimolecular reaction Cl(-) + CH3Cl' → ClCH3 + Cl'(-) in aqueous solution with hybrid density functional theory as the quantum mechanical level. The potential of mean force is calculated, and the free energy barrier is found to be 25.3 kcal/mol, which is in good agreement with previous work. The advantage of the combined EE-MCMM and MM method is that the number of quantum mechanical calculations required for the active subsystem is very small compared to straight direct dynamics.

An Efficient Real Space Multigrid QM/MM Electrostatic Coupling

A popular strategy for simulating large systems where quantum chemical effects are important is the use of mixed quantum mechanical/molecular mechanics methods (QM/MM). While the cost of solving the Schrödinger equation in the QM part is the bottleneck of these calculations, evaluating the Coulomb interaction between the QM and the MM part is surprisingly expensive. In fact it can be just as time-consuming as solving the QM part. We present here a novel real space multigrid approach that handles Coulomb interactions very effectively and implement it in the CP2K code. This novel scheme cuts the cost of this part of the calculation by 2 orders of magnitude. The method does not need very fine-tuning or adjustable parameters, and it is quite accurate, leading to a dynamics with very good energy conservation. We exemplify the validity of our algorithms with simulations of water and of a zwitterionic dipeptide solvated in water.

Fighting Cancer with Transition Metal Complexes: From Naked DNA to Protein and Chromatin Targeting Strategies

Many transition metal complexes have unique physicochemical properties that can be efficiently exploited in medicinal chemistry for cancer treatment. Traditionally, double-stranded DNA has been assumed to be the main binding target; however, recent studies have shown that nucleosomal DNA as well as proteins can act as dominant molecular binding partners. This has raised new questions about the molecular determinants that govern DNA versus protein binding selectivity, and has offered new ways to rationalize their biological activity and possible side effects. To address these questions, molecular simulations at an atomistic level of detail have been used to complement, support, and rationalize experimental data. Herein we review some relevant studies-focused on platinum and ruthenium compounds-to illustrate the power of state-of-the-art molecular simulation techniques and to demonstrate how the interplay between molecular simulations and experiments can make important contributions to elucidating the target preferences of some promising transition metal anticancer agents. This contribution aims at providing relevant information that may help in the rational design of novel drug-discovery strategies.

Peptide Hydrolysis in Thermolysin: Ab Initio QM/MM Investigation of the Glu143-Assisted Water Addition Mechanism

Thermolysin (TLN) is one of the best-studied zinc metalloproteases. Yet the mechanism of action is still under debate. In order to investigate the energetic feasibility of the currently most favored mechanism, we have docked a tripeptide to the active site of TLN and computed the free energy profile at the quantum mechanics/molecular mechanics level of theory. The mechanism consists of three distinct steps:  (i) a Zn-bound water molecule is deprotonated by Glu143 and attacks the carbonyl bond of the substrate; (ii) Glu143 transfers the proton to the amide nitrogen atom; (iii) the nitrogen atom is protonated and the peptide bond is irreversibly broken. The free energy barriers for steps i and iii have almost equal heights, 14.8 and 14.7 kcal/mol, respectively, and are in good agreement with the effective experimental activation barrier obtained for similar substrates, 12.1-13.6 kcal/mol. Transition state stabilization for nucleophilic attack is achieved by formation of a weak coordination bond between the substrate carbonyl oxygen atom and the Zn ion and of three strong hydrogen bonds between the substrate and protonated His231 and two solvent molecules. The transition state for the nucleophilic attack (step i) is more tightly bonded than the enzyme-substrate complex, implying that TLN complies with Pauling's hypothesis regarding transition-state stabilization. Glu143, at first unfavorably oriented for protonation of the amide nitrogen atom, displayed large structural fluctuations that facilitated reorganization of the local hydrogen-bond network and transport of the proton to the leaving group on the nanosecond time scale. The present simulations give further evidence that Glu143 is a highly effective proton shuttle which should be assigned a key role in any reaction mechanism proposed for TLN.

Magnetostructural Dynamics from Hubbard-U Corrected Spin-Projection: [2Fe-2S] Complex in Ferredoxin

A Hubbard-corrected spin-projected two-determinant approach, EBS+Uscf, is introduced to treat low-spin ground states of antiferromagnetically coupled transition metal complexes. In addition to providing access to total energies, forces, and ab initio simulations, it allows one to readily compute Heisenberg's exchange coupling J(t) on the fly. By studying the binuclear [2Fe-2S] cofactor in a metalloprotein, Anabaena Fd, within this consistent nonempirical procedure in combination with a QM/MM framework, it is illustrated that spin-projection, self-interaction corrections, thermal fluctuations, and protein matrix shifts are crucial in obtaining ⟨J⟩ close to the experiment.

Quantitative Assessment of Electrostatic Embedding in Density Functional Theory Calculations of Biomolecular Systems

We evaluate the accuracy of density functional theory quantum calculations of biomolecular subsystems using a simple electrostatic embedding scheme. Our scheme is based on dividing the system of interest into a primary and secondary subsystem. A finite difference discretization of the Kohn-Sham equations is used for the primary subsystem, while its electrostatic environment is modeled with a simple one-electron potential. Force-field atomic partial charges are used to generate smeared Gaussian charge densities and to model the secondary subsystem. We illustrate the utility of this approach with calculations of truncated dipeptide chains. We analyze quantitatively the accuracy of this approach by calculating atomic forces and comparing results with full QM calculations. The impact of the choice made in terminating dangling bonds at the frontier of the QM region is also investigated.

Unraveling the Catalytic Pathway of Metalloenzyme Farnesyltransferase through QM/MM Computation

The protein farnesyltransferase (FTase) is a Zn(2+)-metalloenzyme that catalyzes the farnesylation reaction, i.e., the transfer of the 15-carbon atom farnesyl group from farnesyl diphosphate (FPP) to a specific cysteine of protein substrates. Oncogenic Ras proteins, which are among the FTase substrates, are observed in about 20-30% of human cancer cells. Thus, FTase represents a target for anticancer drug design. Herein, we present a classical force-field-based and quantum mechanics/molecular mechanics (QM/MM) computational study of the FTase reaction mechanism. Our findings offer a detailed picture of the FTase catalytic pathway, describing structural features and the energetics of its saddle points. A moderate dissociation of the diphosphate group from the FPP is observed during the nucleophilic attack of the zinc-bound thiolate. At the transition state, a resonance structure is observed, which indicates the formation of a metastable carbocation. However, no stable intermediate is found along the reaction pathway. Thus, the reaction occurs via an associative mechanism with dissociative character, in agreement with the mechanism proposed by Fierke et al. ( Biochemistry 2000, 39, 2593-2602 and Biochemistry 2003, 42, 9741-9748 ). Moreover, a fluorine-substituted FPP analogue (CF3-FPP) is used to investigate the inhibitory effect of fluorine, which in turn provides additional agreement with experimental data.

The Role of Hydrogen Bonds in the Stabilization of Silver-Mediated Cytosine Tetramers

DNA oligomers can form silver-mediated duplexes, stable in gas phase and solution, with potential for novel biomedical and technological applications. The nucleobase-metal bond primarily drives duplex formation, but hydrogen (H-) bonds may also be important for structure selection and stability. To elucidate the role of H-bonding, we conducted theoretical and experimental studies of a duplex formed by silver-mediated cytosine homopobase DNA strands, two bases long. This silver-mediated cytosine tetramer is small enough to permit accurate, realistic modeling by DFT-based quantum mechanics/molecular mechanics methods. In gas phase, our calculations found two energetically favorable configurations distinguished by H-bonding, one with a novel interplane H-bond, and the other with planar H-bonding of silver-bridged bases. Adding solvent favored silver-mediated tetramers with interplane H-bonding. Overall agreement of electronic circular dichroism spectra for the final calculated structure and experiment validates these findings. Our results can guide use of these stabilization mechanisms for devising novel metal-mediated DNA structures.

Computational metallomics of the anticancer drug cisplatin

Cisplatin, cis-diamminedichlorido-platinum(II), is an important therapeutic tool in the struggle against different tumors, yet it is plagued with the emergence of resistance mechanisms after repeated administrations. This hampers greatly its efficacy. Overcoming resistance problems requires first and foremost an integrated and systematic understanding of the structural determinants and molecular recognition processes involving the drug and its cellular targets. Here we review a strategy that we have followed for the last few years, based on the combination of modern tools from computational chemistry with experimental biophysical methods. Using hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) simulations, validated by spectroscopic experiments (including NMR, and CD), we have worked out for the first time at atomic level the structural determinants in solution of platinated cellular substrates. These include the copper homeostasis proteins Ctr1, Atox1, and ATP7A. All of these proteins have been suggested to influence the pre-target resistance mechanisms. Furthermore, coupling hybrid QM/MM simulations with classical Molecular Dynamics (MD) and free energy calculations, based on force field parameters refined by the so-called "Force Matching" procedure, we have characterized the structural modifications and the free energy landscape associated with the recognition between platinated DNA and the protein HMGB1, belonging to the chromosomal high-mobility group proteins HMGB that inhibit the repair of platinated DNA. This may alleviate issues relative to on-target resistance process. The elucidation of the mechanisms by which tumors are sensitive or refractory to cisplatin may lead to the discovery of prognostic biomarkers. The approach reviewed here could be straightforwardly extended to other metal-based drugs.

How do Water Solvent and Glutathione Ligands Affect the Structure and Electronic Properties of Au25(SR)18(-)?

The effects of aqueous solvent and biological ligands on the structural and electronic properties of thiolate-protected Au25(SR)18(-) clusters have been studied by performing quantum mechanics/molecular mechanics (QM/MM) simulations. Analysis of bond distances and angles show that the solvated nanocluster experiences modest structural changes, which are reflected as flexibility of the Au core. The hydrophilic glutathione ligands shield the metallic core effectively and distort its symmetry via sterical hindrance effects. We show that the previously reported agreement between the calculated HOMO-LUMO gap of the cluster and the optical measurement is due to cancellation of errors, where the typical underestimation of the theoretical band gap compensates the effect of the missing solvent. The use of a hybrid functional results in a HOMO-LUMO gap value of 1.5 eV for the solvated nanocluster with glutathione ligands, in good agreement with optical measurements. Our results demonstrate that ligand/solvent effects should be considered for a proper comparison between theory and experiment.

Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation

The hemoprotein myoglobin is a model system for the study of protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 femtoseconds, with the C, F, and H helices moving away from the heme cofactor and the E and A helices moving toward it. These collective movements are predicted by hybrid quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, our calculations support the prediction that an immediate collective response of the protein occurs upon ligand dissociation, as a result of heme vibrational modes coupling to global modes of the protein.

Excited States and Photochemistry of Chromophores in the Photoactive Proteins Explored by the Combined Quantum Mechanical and Molecular Mechanical Calculations

A photoactive protein usually contains a unique chromophore that is responsible for the initial photoresponse and functions of the photoactive protein are determined by the interaction between the chromophore and its protein surroundings. The combined quantum mechanical and molecular mechanical (QM/MM) approach is demonstrated to be a very useful tool for exploring structures and functions of a photoactive protein with the chromophore and its protein surroundings treated by the QM and MM methods, respectively. In this review, we summarize the basic formulas of the QM/MM approach and emphasize its applications to excited states and photoreactions of chromophores in rhodopsin protein, photoactive yellow protein, and green fluorescent protein.

Active Site Dynamical Effects in the Hydrogen Transfer Rate-limiting Step in the Catalysis of Linoleic Acid by Soybean Lipoxygenase-1 (SLO-1): Primary and Secondary Isotope Contributions

Using ab initio molecular dynamics (AIMD) simulations that facilitate the treatment of rare events, we probe the active site participation in the rate-determining hydrogen transfer step in the catalytic oxidation of linoleic acid by soybean lipoxygenase-1 (SLO-1). The role of two different active site components is probed. (a) On the hydrogen atom acceptor side of the active site, the hydrogen bonding propensity between the acceptor side hydroxyl group, which is bound to the iron cofactor, and the backbone carboxyl group of isoleucine (residue number 839) is studied toward its role in promoting the hydrogen transfer event. Primary and secondary (H/D) isotope effects are also probed and a definite correlation with subtle secondary H/D isotope effects is found. With increasing average nuclear kinetic energy, the increase in transfer probability is enhanced due to the presence of the hydrogen bond between the backbone carbonyl of I839 and the acceptor oxygen. Further increase in average nuclear kinetic energy reduces the strength of this secondary hydrogen bond which leads to a deterioration in hydrogen transfer rates and finally embrances an Arrhenius-like behavior. (b) On the hydrogen atom donor side, the coupling between vibrational modes predominantly localized on the donor-side linoleic acid group and the reactive mode is probed. There appears to be a qualitative difference in the coupling between modes that belong to linoleic acid and the hydrogen transfer mode, for hydrogen and deuterium transfer. For example, the donor side secondary hydrogen atom is much more labile (by nearly a factor of 5) during deuterium transfer as compared to the case for hydrogen transfer. This appears to indicate a greater coupling between the modes belonging to the linoleic acid scaffold and the deuterium transfer mode and also provides a new rationalization for the abnormal (nonclassical) secondary isotope effect results obtained by Knapp, Rickert, and Klinman in J. Am. Chem. Soc. , 2002 , 124 , 3865 . To substantiate our findings noted in point a above, we have suggested an I839 → A839 or I839 → V839 mutation. This will modify the bulkiness of hydrogen the bonding residue, allowing greater flexibility in the secondary hydrogen bond formation highlighted above and adversely affecting the reaction rate.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

A potential with low point charges for pure siliceous zeolites

A modified CHARMM force-field (ZHB potential) with low point charges for silica was previously proposed by Zimmerman et al. (J. Chem. Theory Comput. 2011, 7, 1695). The ZHB potential is advantageous for quantum mechanics/molecular mechanics simulations as it minimizes the electron spill-out problems. In the same spirit, here we propose a modified ZHB potential (MZHB) by reformulating its bonding potential, while retaining the nonbonding potential as in the ZHB force-field. We show that several structural and dynamic properties of silica, like the IR spectrum, distribution functions, mechanical properties, and negative thermal expansion computed using the MZHB potential agree well with experimental data. Further, transferability of MZHB is also tested for reproducing the crystallographic structures of several polymorphs of silica. © 2015 Wiley Periodicals, Inc.