Exploring the Effects of Mutagenesis on FusionRed by Using Excited-State QM/MM Dynamics and Classical Force Field Simulations.

Fluorescent proteins (FPs) are a powerful tool for examining tissues, cells, and subcellular components in vivo and in vitro. FusionRed is a particular FP variant mutated from mKate2 that, in addition to lower cytotoxicity and aggregation rates, has shown potential for acting as a tunable photoswitch. This was posited to stem partially from the presence of a bulky side chain at position 158 and a further stabilizing residue at position 157. In this work, we apply computational techniques including classical molecular dynamics (MD) and combined quantum mechanics/molecular mechanics simulations (QM/MM) to explore the effect of mutagenesis at these locations in FusionRed on the chromophore structure, the excited-state surface, and relative positional stability of the chromophore in the protein pocket. We find specific connections between the statistical sampling of the underlying protein structure and the nonradiative decay mechanisms from excited-state dynamics. A single mutation (C158I) that restricts the motion of the chromophore through a favorable hydrophobic interaction corresponds to an increase in fluorescence quantum yield (FQY), while a second rescue mutation (C158I-A157N) partially restores the flexibility of the chromophore and photoswitchability with favorable water interactions on the surface of the protein that counteracts the original interaction. We suggest that applying this understanding of structural features that inhibit or favor rotation on the excited state can be applied for rational design of new, tunable and red photoswitches.

AutoParams: An Automated Web-Based Tool To Generate Force Field Parameters for Molecular Dynamics Simulations.

Many research questions benefit from molecular dynamics simulations to observe the motions and conformations of molecules over time, which rely on force fields that describe sets of common molecules by category. With the increase of importance for large data sets used in machine learning and growing computational efficiency, the ability to rapidly create large numbers of force field inputs is of high importance. Unusual molecules, such as nucleotide analogues, functionalized carbohydrates, and modified amino acids, are difficult to describe consistently using standard force fields, requiring the development of custom parameters for each unique molecule. While these parameters may be created by individual users, the process can become time-consuming or may introduce errors that may not be immediately apparent. We present an open-source automated parameter generation service, AutoParams, which requires minimal input from the user and creates useful Amber force field parameter sets for most molecules, particularly those that combine molecular types (e.g., a carbohydrate functionalized with a benzene). We include hierarchical atom-typing logic that makes it straightforward to expand with additional force fields and settings, and options for creating monomers in polymers, such as functionalized amino acids. It can be straightforwardly linked to any charge generation program and currently has interfaces to Psi4, PsiRESP, and TeraChem. It is open source and is available via GitHub. It includes error checking and testing protocols to ensure the parameters will be sufficient for subsequent molecular dynamics simulations and streamlines the creation of force field databases.

Characterizing Hydrogen-Bond Interactions in Pyrazinetetracarboxamide Complexes: Insights from Experimental and Quantum Topological Analyses.

Experimental and topological analyses of dipalladium(II) complexes with pyrazinetetracarboxamide ligands containing tetraethyl (1), tetrahexyl (2), and tetrakis(2-hydroxyethyl) ethyl ether (3) are described. The presence of two very short O---O distances between adjacent amide carbonyl groups in the pincer complexes revealed two protons, which necessitated two additional anions to satisfy charge requirements. The results of the crystal structures indicate carbonyl O---O separations approaching that of low barrier hydrogen bonds, ranging from 2.413(5) to 2.430(3) Å. Solution studies and quantum topological analyses, the latter including electron localization function, noncovalent interaction, and Bader's quantum theory of atoms in molecules, were carried out to probe the nature of the short hydrogen bonds and the influence of the ligand environment on their strength. Findings indicated that the ligand field, and, in particular, the counterion at the fourth coordination site, may play a subtle role in determining the degree of covalent association of the bridging protons with one or the other carbonyl groups.

Revealing Drivers for Carboxy-S-adenosyl-l-methionine Use by Neomorphic Variants of a DNA Methyltransferase.

Methylation of DNA plays a key role in diverse biological processes spanning from bacteria to mammals. DNA methyltransferases (MTases) typically employ S-adenosyl-l-methionine (SAM) as a critical cosubstrate and the relevant methyl donor for modification of the C5 position of cytosine. Recently, work on the CpG-specific bacterial MTase, M.MpeI, has shown that a single N374K point mutation can confer the enzyme with the neomorphic ability to use the sparse, naturally occurring metabolite carboxy-S-adenosyl-l-methionine (CxSAM) in order to generate the unnatural DNA modification, 5-carboxymethylcytosine (5cxmC). Here, we aimed to investigate the mechanistic basis for this DNA carboxymethyltransferase (CxMTase) activity by employing a combination of computational modeling and in vitro characterization. Modeling of substrate interactions with the enzyme variant allowed us to identify a favorable salt bridge between CxSAM and N374K that helps to rationalize selectivity of the CxMTase. Unexpectedly, we also discovered a potential role for a key active site E45 residue that makes a bidentate interaction with the ribosyl sugar of CxSAM, located on the opposite face of the CxMTase active site. Prompted by these modeling results, we further explored the space-opening E45D mutation and found that the E45D/N374K double mutant in fact inverts selectivity, preferring CxSAM over SAM in biochemical assays. These findings provide new insight into CxMTase active site architecture and may offer broader utility given the numerous opportunities offered by using SAM analogs for selective molecular labeling in concert with nucleic acid or even protein-modifying MTases.

Combining Evolutionary Conservation and Quantum Topological Analyses To Determine Quantum Mechanics Subsystems for Biomolecular Quantum Mechanics/Molecular Mechanics Simulations.

Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.

Divergence in Dimerization and Activity of Primate APOBEC3C.

The APOBEC3 (A3) family of single-stranded DNA cytidine deaminases are host restriction factors that inhibit lentiviruses, such as HIV-1, in the absence of the Vif protein that causes their degradation. Deamination of cytidine in HIV-1 (-)DNA forms uracil that causes inactivating mutations when uracil is used as a template for (+)DNA synthesis. For APOBEC3C (A3C), the chimpanzee and gorilla orthologues are more active than human A3C, and we determined that Old World Monkey A3C from rhesus macaque (rh) is not active against HIV-1. Biochemical, virological, and coevolutionary analyses combined with molecular dynamics simulations showed that the key amino acids needed to promote rhA3C antiviral activity, 44, 45, and 144, also promoted dimerization and changes to the dynamics of loop 1, near the enzyme active site. Although forced evolution of rhA3C resulted in a similar dimer interface with hominid A3C, the key amino acid contacts were different. Overall, our results determine the basis for why rhA3C is less active than human A3C and establish the amino acid network for dimerization and increased activity. Based on identification of the key amino acids determining Old World Monkey antiviral activity we predict that other Old World Monkey A3Cs did not impart anti-lentiviral activity, despite fixation of a key residue needed for hominid A3C activity. Overall, the coevolutionary analysis of the A3C dimerization interface presented also provides a basis from which to analyze dimerization interfaces of other A3 family members.

Computational Investigation of APOBEC3H Substrate Orientation and Selectivity.

APOBEC3H is a cytidine deaminase protein most well-known for its involvement in antiretroviral activity in humans. It acts upon a single stranded DNA (ssDNA) substrate with preferential targeting of a 5'-TCA-3' motif. Currently available crystal structures do not include the ssDNA substrate in the A3H system, nor is the mechanism of recognition for the preferred sequence known. To determine the position and orientation of the substrate in the active site, we used high-performance computing to perform molecular dynamics simulations on several systems of APOBEC3H. We examined different DNA sequences in the active site to determine the structural and chemical mechanism by which the preferred sequence is recognized. We found residues N49, K50, K51, and K52 to be relevant to the recognition of 3'-adenine and residues S86 and S87 to be relevant to the recognition of 5'-thymine, with both recognitions primarily driven by electrostatic nonbonded interactions.

Single-nucleotide polymorphism of the DNA cytosine deaminase APOBEC3H haplotype I leads to enzyme destabilization and correlates with lung cancer.

A number of APOBEC family DNA cytosine deaminases can induce mutations in tumor cells. APOBEC3H haplotype I is one of the deaminases that has been proposed to cause mutations in lung cancer. Here, we confirmed that APOBEC3H haplotype I can cause uracil-induced DNA damage in lung cancer cells that results in γH2AX foci. Interestingly, the database of cancer biomarkers in DNA repair genes (DNArCdb) identified a single-nucleotide polymorphism (rs139298) of APOBEC3H haplotype I that is involved in lung cancer. While we thought this may increase the activity of APOBEC3H haplotype I, instead we found through computational modeling and cell-based experiments that this single-nucleotide polymorphism causes the destabilization of APOBEC3H Haplotype I. Computational analysis suggests that the resulting K121E change affects the structure of APOBEC3H leading to active site disruption and destabilization of the RNA-mediated dimer interface. A K117E mutation in a K121E background stabilized the APOBEC3H haplotype I, thus enabling biochemical study. Subsequent analysis showed that K121E affected catalytic activity, single-stranded DNA binding and oligomerization on single-stranded DNA. The destabilization of a DNA mutator associated with lung cancer supports the model that too much APOBEC3-induced mutation could result in immune recognition or death of tumor cells.