Parametrized quantum-mechanical approaches combined with the fragment molecular orbital method

Analysis of Site Energies and Excitonic Couplings: The Role of Symmetry and Polarization

An excitonic coupling model is developed based on an equation-of-motion coupled cluster combined with the fragment molecular orbital method. The effects of polarization and excitonic coupling on the splitting of quasi-degenerate levels in systems containing multiple chromophores are elucidated on dimers of formaldehyde, water, formic acid, hydrogen fluoride, and carbon monoxide. It is shown that the level structure is mainly determined by the mutual polarization of chromophores and to a lesser extent by the excitonic coupling. The role of symmetry in excitonic coupling in dimers is discussed. The excitonic coupling between all residues in the photoactive yellow protein (PDB: 2PHY) is analyzed.

Site-Specific Ionization Potentials and Electron Affinities in Large Molecular Systems at Coupled Cluster Level

A many-body expansion of ionization potentials and electron affinities is developed based on a combination of the fragment molecular orbital method and equation-of-motion coupled-cluster (EOM-CC). In addition to site-specific values, obtained as one-body properties, pair and triple corrections are added to account for nonlocal EOM-CC contributions of the molecular environment of a chromophore. The developed method is applied to carboxylic acids, alkyl cations, a protein ubiquitin (Protein Data Bank ID 1UBQ), and a nano ribbon of white graphene elucidating the effect of environment on ionization potential and electron affinity.

How E-, L-, and P-Selectins Bind to sLex and PSGL-1: A Quantification of Critical Residue Interactions

Selectins and their ability to interact with specific ligands are a cornerstone in cell communication. Over the last three decades, a considerable wealth of experimental and molecular modeling insights into their structure and modus operandi were gathered. Nonetheless, explaining the role of individual selectin residues on a quantitative level remained elusive, despite its importance in understanding the structure-function relationship in these molecules and designing their inhibitors. This work explores essential interactions of selectin-ligand binding, employing a multiscale approach that combines molecular dynamics, quantum-chemical calculations, and residue interaction network models. Such an approach successfully reproduces most of the experimental findings. It proves to be helpful, with the potential for becoming an established tool for quantitative predictions of residue contribution to the binding of biomolecular complexes. The results empower us to quantify the importance of particular residues and functional groups in the protein-ligand interface and to pinpoint differences in molecular recognition by the three selectins. We show that mutations in the E-, L-, and P-selectins, e.g., different residues in positions 46, 85, 97, and 107, present a crucial difference in how the ligand is engaged. We assess the role of sulfation of tyrosine residues in PSGL-1 and suggest that TyrSO3- in position 51 interacting with Arg85 in P-selectin is a significant factor in the increased affinity of P-selectin to PSGL-1 compared to E- and L-selectins. We propose an original pharmacophore targeting five essential PSGL-binding sites based on the analysis of the selectin···PSGL-1 interactions.

Analytic Gradient for Time-Dependent Density Functional Theory Combined with the Fragment Molecular Orbital Method

The analytic energy gradient of energy with respect to nuclear coordinates is derived for the fragment molecular orbital (FMO) method combined with time-dependent density functional theory (TDDFT). The response terms arising from the use of a polarizable embedding are derived. The obtained analytic FMO-TDDFT gradient is shown to be accurate in comparison to both numerical FMO-TDDFT and unfragmented TDDFT gradients, at the level of two- and three-body expansions. The gradients are used for geometry optimizations, molecular dynamics, vibrational calculations, and simulations of IR and Raman spectra of excited states. The developed method is used to optimize the geometry of the ground and excited electronic states of the photoactive yellow protein (PDB: 2PHY).