Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes the coronavirus disease 2019 (COVID-19), spread rapidly around the globe. The main protease encoded by SARS-CoV-2 is essential for processing of the polyproteins translated from the viral RNA genome, making this protein a potential drug target. A recently reported mutation in the protease, P108S, may be responsible for milder symptoms observed in COVID-19 patients in Tokyo. Starting from a crystal structure of the SARS-CoV-2 main protease in the dimeric form, we performed triplicate 5.0-μs molecular dynamics simulations of the wild-type and P108S mutant. Our computational results suggest a link between the mutation P108S and dynamics of the catalytic sites in the main protease: The catalytic dyad become considerably inaccessible to substrates in the P108S mutant. Our results also demonstrate the potential of molecular dynamics simulations to complement experimental techniques and other computational methods, such as protein design calculations, which predict effects of mutations based on static crystal structures. Further studies are certainly necessary to quantitively understand the relationships between the P108S mutation and physical properties of the main protease, but the results of our study will immediately inform development of new protease inhibitors.
In the present study, we provide a reformulation of the theory originally proposed by Förster which allows for simple and convenient formulas useful to estimate the relative contributions of transition dipole moments of a donor and acceptor (chemical factors), their orientation factors (intermolecular structural factors), intermolecular center-to-center distances (intermolecular structural factors), spectral overlaps of absorption and emission spectra (photophysical factors), and refractive index (material factor) to the excitation energy transfer (EET) rate constant. To benchmark their validity, we focused on the EET occurring in C-phycocyanin (C-PC) chromophores. To this aim, we resorted to quantum chemistry calculations to get optimized molecular structures of the C-PC chromophores within the density functional theory (DFT) framework. The absorption and emission spectra, as well as transition dipole moments, were computed by using the time-dependent DFT (TDDFT). Our method was applied to several types of C-PCs showing that the EET rates are determined by an interplay of their specific physical, chemical, and geometrical features. These results show that our formulas can become a useful tool for a reliable estimation of the relative contributions of the factors regulating the EET transfer rate.
Cytochrome c oxidase (CcO), a terminal oxidase in the respiratory chain, catalyzes the reduction of O2 to water coupled with the proton pump across the membrane. Mitochondrial CcO exists in monomeric and dimeric forms, and as a monomer as part of the respiratory supercomplex, although the enzymatic reaction proceeds in the CcO monomer. Recent biochemical and crystallographic studies of monomeric and dimeric CcOs have revealed functional and structural differences among them. In solubilized mitochondrial membrane, the monomeric form is dominant, and a small amount of dimer is observed. The activity of the monomeric CcO is higher than that of the dimer, suggesting that the monomer is the active form. In the structure of monomeric CcO, a hydrogen bond network of water molecules is formed at the entrance of the proton transfer K-pathway, and in dimeric CcO, this network is altered by a cholate molecule binding between monomers. The specific binding of the cholate molecule at the dimer interface suggests that the binding of physiological ligands similar in size or shape to cholate could also trigger dimer formation as a physiological standby form. Because the dimer interface also contains weak interactions of nonspecifically bound lipid molecules, hydrophobic interactions between the transmembrane helices, and a Met-Met interaction between the extramembrane regions, these interactions could support the stabilization of the standby form. Structural analyses also suggest that hydrophobic interactions of cardiolipins bound to the trans-membrane surface of CcO are involved in forming the supercomplex.
Middle rhodopsin (MR) found from the archaeon Haloquadratum walsbyi is evolutionarily located between two different types of rhodopsins, bacteriorhodopsin (BR) and sensory rhodopsin II (SRII). Some isomers of the chromophore retinal and the photochemical reaction of MR are markedly different from those of BR and SRII. In this study, to obtain the structural information regarding its active center (i.e., retinal), we subjected MR embedded in lipid bilayers to solid-state magic-angle spinning nuclear magnetic resonance (NMR) spectroscopy. The analysis of the isotropic 13C chemical shifts of the retinal chromophore revealed the presence of three types of retinal configurations of dark-adapted MR: (13-trans, 15-anti (all-trans)), (13-cis, 15-syn), and 11-cis isomers. The higher field resonance of the 20-C methyl carbon in the all-trans retinal suggested that Trp182 in MR has an orientation that is different from that in other microbial rhodopsins, owing to the changes in steric hindrance associated with the 20-C methyl group in retinal. 13Cζ signals of Tyr185 in MR for all-trans and 13-cis, 15-syn isomers were discretely observed, representing the difference in the hydrogen bond strength of Tyr185. Further, 15N NMR analysis of the protonated Schiff base corresponding to the all-trans and 13-cis, 15-syn isomers in MR showed a strong electrostatic interaction with the counter ion. Therefore, the resulting structural information exhibited the property of stable retinal conformations of dark-adapted MR.
A cutinase-like enzyme from Saccharomonospora viridis AHK190, Cut190, can depolymerize polyethylene terephthalate (PET). As high activity at approximately 70°C is required for PET depolymerization, structure-based protein engineering of Cut190 was carried out. Crystal structure information of the Cut190 mutants was used for protein engineering and for evaluating the molecular basis of activity and thermal stability. A variety of biophysical methods were employed to unveil the mechanisms underlying the unique features of Cut190, which included the regulation of its activity and thermal stability by Ca2+. Ca2+ association and dissociation can change the enzyme conformation to regulate catalytic activity. Weak metal-ion binding would be required for the naïve conformational change of Cut190, while maintaining its fluctuation, to "switch" the enzyme on and off. The activity of Cut190 is regulated by the weak Ca2+ binding to the specific site, Site 1, while thermal stability is mainly regulated by binding to another Site 2, where a disulfide bond could be introduced to increase the stability. Recent results on the structure-activity relationship of engineered Cut190 are reviewed, including the application for PET depolymerization by enzymes.
ABEGO is a coarse-grained representation for poly-peptide backbone dihedral angles. The Ramachandran map is divided into four segments denoted as A, B, E, and G to represent the local conformation of poly-peptide chains in the character strings. Although the ABEGO representation is widely used in backbone building simulation for de novo protein design, it cannot capture minor differences in backbone dihedral angles, which potentially leads to ambiguity between two structurally distinct fragments. Here, I show a nontrivial example of two local motifs that could not be distinguished by their ABEGO representations. I found that two well-known local motifs αα-hairpins and αα-corners are both represented as α-GBB-α and thus indistinguishable in the ABEGO representation, although they show distinct arrangements of the flanking α-helices. I also found that α-GBB-α motifs caused a loss of efficiency in the ABEGO-based fragment-assembly simulations for de novo protein backbone design. Nevertheless, I was able to design amino-acid sequences that were predicted to fold into the target topologies that contained these α-GBB-α motifs, which suggests such topologies that are difficult to build by ABEGO-based simulations are designable once the backbone structures are modeled by some means. The finding that certain local motifs bottleneck the ABEGO-based fragment-assembly simulations for construction of backbone structures suggests that finer representations of backbone torsion angles are required for efficiently generating diverse topologies containing such indistinguishable local motifs.
The green fluorescent protein (GFP) derived from Pacific Ocean jellyfish is an essential tool in biology. GFP-solvent interactions can modulate the fluorescent property of GFP. We previously reported that glycine insertion is an effective mutation in the yellow variant of GFP, yellow fluorescent protein (YFP). Glycine insertion into one of the β-strands comprising the barrel structure distorts its structure, allowing water molecules to invade near the chromophore, enhancing hydrostatic pressure or solution hydrophobicity sensitivity. However, the underlying mechanism of how glycine insertion imparts environmental sensitivity to YFP has not been elucidated yet. To unveil the relationship between fluorescence and β-strand distortion, we investigated the effects of glycine insertion on the dependence of the optical properties of GFP variants named enhanced-GFP (eGFP) and its yellow (eYFP) and cyan (eCFP) variants with respect to pH, temperature, pressure, and hydrophobicity. Our results showed that the quantum yield decreased depending on the number of inserted glycines in all variants, and the dependence on pH, temperature, pressure, and hydrophobicity was altered, indicating the invasion of water molecules into the β-barrel. Peak shifts in the emission spectrum were observed in glycine-inserted eGFP, suggesting a change of the electric state in the excited chromophore. A comparative investigation of the spectral shift among variants under different conditions demonstrated that glycine insertion rearranged the hydrogen bond network between His148 and the chromophore. The present results provide important insights for further understanding the fluorescence mechanism in GFPs and suggest that glycine insertion could be a potent approach for investigating the relationship between water molecules and the intra-protein chromophore.
Synchronized movement of (both unicellular and multicellular) systems can be observed almost everywhere. Understanding of how organisms are regulated to synchronized behavior is one of the challenging issues in the field of collective motion. It is hypothesized that one or a few agents in a group regulate(s) the dynamics of the whole collective, known as leader(s). The identification of the leader (influential) agent(s) is very crucial. This article reviews different mathematical models that represent different types of leadership. We focus on the improvement of the leader-follower classification problem. It was found using a simulation model that the use of interaction domain information significantly improves the leader-follower classification ability using both linear schemes and information-theoretic schemes for quantifying influence. This article also reviews different schemes that can be used to identify the interaction domain using the motion data of agents.

