The MoD-QM/MM methodology for structural refinement of photosystem II and other biological macromolecules

Molecular Structure Related to an S = 5/2 High-Spin S2 State Manganese Cluster of Photosystem II Investigated by Q-Band Pulse EPR Spectroscopy

The high-spin S₂ state of the photosynthetic oxygen-evolving cluster Mn₄CaO₅, corresponding to the g = 4.1 signal for X-band electron paramagnetic resonance (EPR), was investigated using Q-band pulsed EPR, which detected a main peak at g = 3.10 and satellite peaks at 5.25, 4.55, and 2.80. We evaluated the spin state as the zero-field splitting of D = 0.465 cm^-1 and E/D = 0.245 with S = 5/2. The temperature dependence of the T₁ relaxation time revealed that the excited-state energy was 28.7 cm^-1 higher than that of the S = 5/2 ground state. By comparing present quantum mechanical (QM) calculation models, a closed-cubane structure with the protonation state of two oxygens, W1 (= OH^-) and W2 (= H₂O), was the most probable structure for the S = 5/2 state. The three-pulse electron spin-echo envelope modulation (ESEEM) detected the nuclear signal, which was assigned to nitrogen as His332 ligated to the Mn1 ion. The obtained hyperfine constant for the nitrogen signal was significantly reduced from that in the S = 1/2 low-spin state. These results indicate that the S = 5/2 spin state arises from the closed-cubane structure.

Computational vaccinology guided design of multi-epitope subunit vaccine against a neglected arbovirus of the Americas

Mayaro virus (MAYV) is an arbovirus found in the Americas that can cause debilitating arthritogenic disease. Although it is an emerging virus, the only current approach is vector control, as there are no approved vaccines to prevent MAYV infection nor therapeutics to treat it. In search of an effective vaccine candidate against MAYV, we used immunoinformatics and molecular modeling to attempt to identify promiscuous T-cell epitopes of the nonstructural polyproteins (nsP1, nsP2, nsP3, and nsP4) from 127 MAYV genomes sequenced in the Americas (08 Bolivia, 72 Brazil, 04 French Guiana, 05 Haiti, 20 Peru, 04 Trinidad and Tobago, and 14 Venezuela). For this purpose, consensus sequences of 360 proteins were used to identify short protein sequences that can bind to MHC I class (MHC II). Our analysis revealed 56 potential MHC-I/TCD8+ (29 MHC-II/TCD4+) epitopes, but only 6 (16) TCD8+ (TCD4+) epitopes showed high antigenicity and conservation, non-allergenicity, non-toxicity, and excellent population coverage. Finally, classical and quantum mechanical calculations (QM:MM) were used to improve the quality of the docking calculations, with the QM part of the simulations performed using the density functional theory formalism (DFT). These results provide insights for the advancement of diagnostic platforms, vaccine development, and immunotherapeutic interventions.Communicated by Ramaswamy H. Sarma.

Flexible boundary layer using exchange for embedding theories. I. Theory and implementation

Embedding theory is a powerful computational chemistry approach to exploring the electronic structure and dynamics of complex systems, with Quantum Mechanical/Molecular Mechanics (QM/MM) being the prime example. A challenge arises when trying to apply embedding methodology to systems with diffusible particles, e.g., solvents, if some of them must be included in the QM region, for example, in the description of solvent-supported electronic states or reactions involving proton transfer or charge-transfer-to-solvent: without a special treatment, inter-diffusion of QM and MM particles will eventually lead to a loss of QM/MM separation. We have developed a new method called Flexible Boundary Layer using Exchange (FlexiBLE) that solves the problem by adding a biasing potential to the system that closely maintains QM/MM separation. The method rigorously preserves ensemble averages by leveraging their invariance to an exchange of identical particles. With a careful choice of the biasing potential and the use of a tree algorithm to include only important QM and MM exchanges, we find that the method has an MM-forcefield-like computational cost and thus adds negligible overhead to a QM/MM simulation. Furthermore, we show that molecular dynamics with the FlexiBLE bias conserves total energy, and remarkably, sub-diffusional dynamical quantities in the inner QM region are unaffected by the applied bias. FlexiBLE thus widens the range of chemistry that can be studied with embedding theory.

Structural Predictive Model of Presenilin-2 Protein and Analysis of Structural Effects of Familial Alzheimer's Disease Mutations

Alzheimer's disease manifests itself in brain tissue by neuronal death, due to aggregation of β-amyloid, produced by senile plaques, and hyperphosphorylation of the tau protein, which produces neurofibrillary tangles. One of the genetic markers of the disease is the gene that translates the presenilin-2 protein, which has mutations that favor the appearance of the disease and has no reported crystallographic structure. In view of this, protein modeling is performed using prediction and structural refinement tools followed by an energetic and stereochemical characterization for its validation. For the simulation, four reported mutations are chosen, which are Met239Ile, Met239Val, Ser130Leu, and Thr122Arg, all associated with various functional responses. From a theoretical analysis, a preliminary bioinformatic study is made to find the phosphorylation patterns in the protein and the hydropathic index according to the polarity and chemical environment. Molecular visualization was carried out with the Chimera 1.14 software, and the theoretical calculation with the hybrid quantum mechanics/molecular mechanics system from the semi-empirical method, with Spartan18 software and an AustinModel1 basis. These relationships allow for studying the system from a structural approach with the determination of small distance changes, potential surfaces, electrostatic maps, and angle changes, which favor the comparison between wild-type and mutant systems. With the results obtained, it is expected to complement experimental data reported in the literature from models that would allow us to understand the effects of the selected mutations.

Quantum Chemical Calculation of Molecular and Periodic Peptide and Protein Structures

Special-purpose classical force fields (FFs) provide good accuracy at very low computational cost, but their application is limited to systems for which potential energy functions are available. This excludes most metal-containing proteins or those containing cofactors. In contrast, the GFN2-xTB semiempirical quantum chemical method is parametrized for almost the entire periodic table. The accuracy of GFN2-xTB is assessed for protein structures with respect to experimental X-ray data. Furthermore, the results are compared with those of two special-purpose FFs, HF-3c, PM6-D3H4X, and PM7. The test sets include proteins without any prosthetic groups as well as metalloproteins. Crystal packing effects are examined for a set of smaller proteins to validate the molecular approach. For the proteins without prosthetic groups, the special purpose FF OPLS-2005 yields the smallest overall RMSD to the X-ray data but GFN2-xTB provides similarly good structures with even better bond-length distributions. For the metalloproteins with up to 5000 atoms, a good overall structural agreement is obtained with GFN2-xTB. The full geometry optimizations of protein structures with on average 1000 atoms in wall-times below 1 day establishes the GFN2-xTB method as a versatile tool for the computational treatment of various biomolecules with a good accuracy/computational cost ratio.

Substrate water exchange in the S2 state of photosystem II is dependent on the conformation of the Mn4Ca cluster

The structural flexibility of the Mn₄Ca cluster in photosystem II supports the exchange of the central O5 bridge.

Quantum mechanics implementation in drug-design workflows: does it really help?

The pharmaceutical industry is progressively operating in an era where development costs are constantly under pressure, higher percentages of drugs are demanded, and the drug-discovery process is a trial-and-error run. The profit that flows in with the discovery of new drugs has always been the motivation for the industry to keep up the pace and keep abreast with the endless demand for medicines. The process of finding a molecule that binds to the target protein using in silico tools has made computational chemistry a valuable tool in drug discovery in both academic research and pharmaceutical industry. However, the complexity of many protein-ligand interactions challenges the accuracy and efficiency of the commonly used empirical methods. The usefulness of quantum mechanics (QM) in drug-protein interaction cannot be overemphasized; however, this approach has little significance in some empirical methods. In this review, we discuss recent developments in, and application of, QM to medically relevant biomolecules. We critically discuss the different types of QM-based methods and their proposed application to incorporating them into drug-design and -discovery workflows while trying to answer a critical question: are QM-based methods of real help in drug-design and -discovery research and industry?

Effects of aligned α-helix peptide dipoles on experimental electrostatic potentials

Aligned protein α-helix dipoles have been implicated in protein function and structure. The recent breakthroughs in high-resolution electron microscopy (EM) of macromolecules makes it possible to explore fundamental aspects of structural biology at the detailed molecular level. The electrostatic potential (ESP) generated by aligned protein α-helix dipole should be observable in high-resolution EM maps despite the fact that the effect may be partially screened by induced electric fields. Here, we show that aligned backbone dipoles in protein α-helices account for long-range features in the protein ESP functions. Our results are consistent with experimental EM maps and density functional theory calculations, including direct Fourier summation for proper calculation of the ESP due to the nonlocal nature of the ESP function from aligned dipoles and other partial atomic charges.

Alternative mechanisms for O2release and O–O bond formation in the oxygen evolving complex of photosystem II

A new transition state for O₂release has been found. An alternative, nearly competitive, mechanism for O–O bond formation is described, which is very similar to the one previously suggested.

Convergence of QM/MM and Cluster Models for the Spectroscopic Properties of the Oxygen-Evolving Complex in Photosystem II

The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

Water oxidation mechanism in photosystem II, including oxidations, proton release pathways, O-O bond formation and O2 release

The present status of DFT studies on water oxidation in photosystem II is described. It is argued that a full understanding of all steps is close. In each S-transition, the manganese that is oxidized and the proton released are strongly implicated, and structures of all intermediates have been determined. For the S2-state, recent important experimental findings support key elements of the structure and the mechanism. In this mechanism, the O-O bond is formed between an oxyl radical in the center of the cluster and an Mn-bridging μ-oxo ligand, which was suggested already in 2006. The DFT structure of the oxygen evolving complex, suggested in 2008, is very similar to the recent high-resolution X-ray structure. Some new aspects of the interaction between P680 and the OEC are suggested. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

3D-RISM-Dock: A New Fragment-Based Drug Design Protocol

We explore a new approach in the rational design of specificity in molecular recognition of small molecules based on statistical-mechanical integral equation theory of molecular liquids in the form of the three-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH). The numerically stable iterative solution of conventional 3D-RISM equations includes the fragmental decomposition of flexible ligands, which are treated as distinct species in solvent mixtures of arbitrary complexity. The computed density functions for solution (including ligand) molecules are obtained as a set of discrete spatial grids that uniquely describe the continuous solvent-site distribution around the protein solute. Potentials of mean force derived from these distributions define the scoring function interfaced with the AutoDock program for an automated ranking of docked conformations. As a case study in terms of solvent composition, we analyze cooperative interactions encountered in the binding of a flexible thiamine molecule to the prion protein at near-physiological conditions. The predicted location and residency times of computed binding modes are in excellent agreement with the available experimental data.

QM/MM Prediction of the Stark Shift in the Active Site of a Protein

Recent developments in the biophysical characterization of proteins have provided a means of directly measuring electrostatic fields by introducing a probe molecule to the system of interest and interpreting photon absorption in the context of the Stark effect. To fully account for this effect, the development of accurate atomistic models is of paramount importance. However, suitable computational protocols for evaluating Stark shifts in proteins are yet to be established. In this work, we present a comprehensive computational method to predict the change in absorption frequency of a probe functional group as a direct result of a perturbation in its surrounding electrostatic field created by a protein environment, i.e., the Stark shift. We apply the method to human aldose reductase, a key protein enzyme that catalyzes the reduction of monosaccharides. We develop a protocol based on a combination of molecular dynamics and moving-domain QM/MM methods, which achieves quantitative agreement with experiment. We outline the difficulties in predicting localized electrostatic field changes within a protein environment, and by extension the Stark shift, due to a protein site mutation. Furthermore, the combined use of Stark effect spectroscopy and computational modeling is used to predict the protonation state of ionizable residues in the vicinity of the electrostatic probe.

Theoretical evaluation of structural models of the S2 state in the oxygen evolving complex of Photosystem II: protonation states and magnetic interactions

Protonation states of water ligands and oxo bridges are intimately involved in tuning the electronic structures and oxidation potentials of the oxygen evolving complex (OEC) in Photosystem II, steering the mechanistic pathway, which involves at least five redox state intermediates S(n) (n = 0-4) resulting in the oxidation of water to molecular oxygen. Although protons are practically invisible in protein crystallography, their effects on the electronic structure and magnetic properties of metal active sites can be probed using spectroscopy. With the twin purpose of aiding the interpretation of the complex electron paramagnetic resonance (EPR) spectroscopic data of the OEC and of improving the view of the cluster at the atomic level, a complete set of protonation configurations for the S(2) state of the OEC were investigated, and their distinctive effects on magnetic properties of the cluster were evaluated. The most recent X-ray structure of Photosystem II at 1.9 Å resolution was used and refined to obtain the optimum structure for the Mn(4)O(5)Ca core within the protein pocket. Employing this model, a set of 26 structures was constructed that tested various protonation scenarios of the water ligands and oxo bridges. Our results suggest that one of the two water molecules that are proposed to coordinate the outer Mn ion (Mn(A)) of the cluster is deprotonated in the S(2) state, as this leads to optimal experimental agreement, reproducing the correct ground state spin multiplicity (S = 1/2), spin expectation values, and EXAFS-derived metal-metal distances. Deprotonation of Ca(2+)-bound water molecules is strongly disfavored in the S(2) state, but dissociation of one of the two water ligands appears to be facile. The computed isotropic hyperfine couplings presented here allow distinctions between models to be made and call into question the assumption that the largest coupling is always attributable to Mn(III). The present results impose limits for the total charge and the proton configuration of the OEC in the S(2) state, with implications for the cascade of events in the Kok cycle and for the water splitting mechanism.

Quantum chemical description of absorption properties and excited-state processes in photosynthetic systems

The theoretical description of the initial steps in photosynthesis has gained increasing importance over the past few years. This is caused by more and more structural data becoming available for light-harvesting complexes and reaction centers which form the basis for atomistic calculations and by the progress made in the development of first-principles methods for excited electronic states of large molecules. In this Review, we discuss the advantages and pitfalls of theoretical methods applicable to photosynthetic pigments. Besides methodological aspects of excited-state electronic-structure methods, studies on chlorophyll-type and carotenoid-like molecules are discussed. We also address the concepts of exciton coupling and excitation-energy transfer (EET) and compare the different theoretical methods for the calculation of EET coupling constants. Applications to photosynthetic light-harvesting complexes and reaction centers based on such models are also analyzed.

Structural models of the manganese complex of photosystem II and mechanistic implications

Photosynthetic water oxidation and O₂ formation are catalyzed by a Mn₄Ca complex bound to the proteins of photosystem II (PSII). The catalytic site, including the inorganic Mn₄CaO(n)H(x) core and its protein environment, is denoted as oxygen-evolving complex (OEC). Earlier and recent progress in the endeavor to elucidate the structure of the OEC is reviewed, with focus on recent results obtained by (i) X−ray spectroscopy (specifically by EXAFS analyses), and (ii) X-ray diffraction (XRD, protein crystallography). Very recently, an impressive resolution of 1.9Å has been achieved by XRD. Most likely however, all XRD data on the Mn₄CaO(n)H(x) core of the OEC are affected by X-ray induced modifications (radiation damage). Therefore and to address (important) details of the geometric and electronic structure of the OEC, a combined analysis of XRD and XAS data has been approached by several research groups. These efforts are reviewed and extended using an especially comprehensive approach. Taking into account XRD results on the protein environment of the inorganic core of the Mn complex, 12 alternative OEC models are considered and evaluated by quantitative comparison to (i) extended-range EXAFS data, (ii) polarized EXAFS of partially oriented PSII membrane particles, and (iii) polarized EXAFS of PSII crystals. We conclude that there is a class of OEC models that is in good agreement with both the recent crystallographic models and the XAS data. On these grounds, mechanistic implications for the O−O bond formation chemistry are discussed. This article is part of a Special Issue entitled: Photosystem II.

QM/MM methods: looking inside heme proteins biochemistry

Mixed quantum mechanics/molecular mechanics methods offer a valuable computational tool for understanding biochemical events. When combined with conformational sampling techniques, they allow for an exhaustive exploration of the enzymatic mechanism. Heme proteins are ubiquitous and essential for every organism. In this review we summarize our efforts towards the understanding of heme biochemistry. We present: 1) results on ligand migration on globins coupled to the ligand binding event, 2) results on the localization of the spin density in compound I of cytochromes and peroxidases, 3) novel methodologies for mapping the electron transfer pathways and 4) novel data on Tryptophan 2,3-dioxygenase. For this enzyme our results strongly indicate that the distal oxygen will end up on the C3 indole carbon, whereas the proximal oxygen will end up in the C2 position. Interestingly, the process involves the formation of an epoxide and a heme ferryl intermediate. The overall energy profile indicates an energy barrier of approximately 18 kcal/mol and an exothermic driving force of almost 80 kcal/mol.