Dynamics of nitric oxide controlled by protein complex in bacterial system.

Nitric oxide (NO) plays diverse and significant roles in biological processes despite its cytotoxicity, raising the question of how biological systems control the action of NO to minimize its cytotoxicity in cells. As a great example of such a system, we found a possibility that NO-generating nitrite reductase (NiR) forms a complex with NO-decomposing membrane-integrated NO reductase (NOR) to efficiently capture NO immediately after its production by NiR in anaerobic nitrate respiration called denitrification. The 3.2-Å resolution structure of the complex of one NiR functional homodimer and two NOR molecules provides an idea of how these enzymes interact in cells, while the structure may not reflect the one in cells due to the membrane topology. Subsequent all-atom molecular dynamics (MD) simulations of the enzyme complex model in a membrane and structure-guided mutagenesis suggested that a few interenzyme salt bridges and coulombic interactions of NiR with the membrane could stabilize the complex of one NiR homodimer and one NOR molecule and contribute to rapid NO decomposition in cells. The MD trajectories of the NO diffusion in the NiR:NOR complex with the membrane showed that, as a plausible NO transfer mechanism, NO released from NiR rapidly migrates into the membrane, then binds to NOR. These results help us understand the mechanism of the cellular control of the action of cytotoxic NO.

Challenges and opportunities in connecting simulations with experiments via molecular dynamics of cellular environments.

Computer simulations are widely used to study molecular systems, especially in biology. As simulations have greatly increased in scale reaching cellular levels there are now significant challenges in managing, analyzing, and interpreting such data in comparison with experiments that are being discussed. Management challenges revolve around storing and sharing terabyte to petabyte scale data sets whereas the analysis of simulations of highly complex systems will increasingly require automated machine learning and artificial intelligence approaches. The comparison between simulations and experiments is furthermore complicated not just by the complexity of the data but also by difficulties in interpreting experiments for highly heterogeneous systems. As an example, the interpretation of NMR relaxation measurements and comparison with simulations for highly crowded systems is discussed.

Mechanistic insight into the blocking of CO diffusion in [NiFe]-hydrogenase mutants through multiscale simulation.

[NiFe]-hydrogenases are fascinating biological catalysts with potential application in biofuel cells. However, a severe problem in practical application is the strong sensitivity of hydrogenase to gaseous inhibitor molecules such as CO and O(2). Recently, a number of successful protein engineering studies have been reported that aimed at lowering the access of diatomic inhibitors to the active site pocket, but the molecular mechanism conferring increased resistance remained unclear. Here we use a multiscale simulation approach combining molecular dynamics with a master equation formalism to explain the steady drop in CO diffusion rate observed for the mutants V74M L122A, V74M L122M, and V74M of Desulfovibrio fructosovorans [NiFe]-hydrogenase. We find that diffusion in these variants is controlled by two gates, one between residues 74 and 476 and the other between residues 74 and 122. The existence of two control points in different locations explains why the reduction in the experimental diffusion rate does not simply correlate with the width of the main gas channel. We also find that in the more effective mutation (V74M) CO molecules are still able to reach the active site through transitions that are gated by the microsecond dihedral motions of the side chain of R476 and the thermal fluctuations of the width of the gas channel defined by M74 and L122. Reflecting on the molecular information gained from simulation, we discuss future mutation experiments that could further lower the diffusion rates of small ligands inhibiting [NiFe]-hydrogenase.

Uncovering a dynamically formed substrate access tunnel in carbon monoxide dehydrogenase/acetyl-CoA synthase.

The transport of small ligands to active sites of proteins is the basis of vital processes in biology such as enzymatic catalysis and cell signaling, but also of more destructive ones including enzyme inhibition and oxidative damage. Here, we show how a diffusion-reaction model solved by means of molecular dynamics and density functional theory calculations provides novel insight into the transport of small ligands in proteins. In particular, we unravel the existence of an elusive, dynamically formed gas channel, which CO2 takes to diffuse from the solvent to the active site (C-cluster) of the bifunctional multisubunit enzyme complex carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS). Two cavities forming this channel are temporarily created by protein fluctuations and are not apparent in the X-ray structures. The ligand transport is controlled by two residues at the end of this tunnel, His113 and His116, and occurs on the same time scale on which chemical binding to the active site takes place (0.1-1 ms), resulting in an overall binding rate on the second time scale. We find that upon reduction of CO2 to CO, the newly formed Fe-hydroxy ligand greatly strengthens the hydrogen-bond network, preventing CO from exiting the protein through the same way that CO2 takes to enter the protein. This is the basis for directional transport of CO from the production site (C-cluster of CODH subunit) to the utilization site (A-cluster of ACS subunit). In view of these results, a general picture emerges of how large proteins guide small ligands toward their active sites.

Multiscale simulation reveals multiple pathways for H2 and O2 transport in a [NiFe]-hydrogenase.

Hydrogenases are enzymes that catalyze the reversible conversion of hydrogen molecules to protons and electrons. The mechanism by which the gas molecules reach the active site is important for understanding the function of the enzyme and may play a role in the selectivity for hydrogen over inhibitor molecules. Here, we develop a general multiscale molecular simulation approach for the calculation of diffusion rates and determination of pathways by which substrate or inhibitor gases can reach the protein active site. Combining kinetic data from both equilibrium simulations and enhanced sampling, we construct a master equation describing the movement of gas molecules within the enzyme. We find that the time-dependent gas population of the active site can be fit to the same phenomenological rate law used to interpret experiments, with corresponding diffusion rates in very good agreement with experimental data. However, in contrast to the conventional picture, in which the gases follow a well-defined hydrophobic tunnel, we find that there is a diverse network of accessible pathways by which the gas molecules can reach the active site. The previously identified tunnel accounts for only about 60% of the total flux. Our results suggest that the dramatic decrease in the diffusion rate for mutations involving the residue Val74 could be in part due to the narrowing of the passage Val74-Arg476, immediately adjacent to the binding site, explaining why mutations of Leu122 had only a negligible effect in experiment. Our method is not specific to the [NiFe]-hydrogenase and should be generally applicable to the transport of small molecules in proteins.

Discharge and photo-luminance properties of a parallel plates electron emission lighting device.

The gas discharge and photo-luminance properties of a planar lighting source featuring highly uniform light emission and mercury-free design were studied. The current density-voltage characteristics and the associated gas discharge of the devices operating with the values of the ratio of electric field to gas pressure (E/p) between 4.3 kV/Torr-cm and 35.7 kV/Torr-cm indicate that the width of the cathode fall extends over the entire gap between the two electrodes and the device is mostly in the obstructed discharge regime. The optical emission analysis confirmed the electron collision-induced gas emissions and strong effect of gas pressure on the phosphor emission when operated at constant current density, both are indicative of the primary roles played by the electron energy.

A microscopic model for gas diffusion dynamics in a [NiFe]-hydrogenase.

We describe and apply a microscopic model for the calculation of gas diffusion rates in a [NiFe]-hydrogenase. This enzyme has attracted much interest for use as a H(2) oxidising catalyst in biofuel cells, but a major problem is their inhibition by CO and O(2). In our model, the diffusive hopping of gas molecules in the protein interior is coarse grained using a master equation approach with transition rates estimated from equilibrium and non-equilibrium pulling simulations. Propagating the rate matrix in time, we find that the probability for a gas molecule to reach the enzyme active site follows a mono-exponential increase. Fits to a phenomenological rate law give an effective diffusion rate constant for CO that is in very good agreement with experimental measurements. We find that CO prefers to move along the canonical 'hydrophobic' main channel towards the active site, in contrast to O(2) and H(2), which were previously shown to explore larger fractions of the protein. Differences in the diffusion of the three gases are discussed in light of recent efforts to engineer a gas selectivity filter in the enzyme.

Glycan engineering of the SARS-CoV-2 receptor-binding domain elicits cross-neutralizing antibodies for SARS-related viruses.

Broadly protective vaccines against SARS-related coronaviruses that may cause future outbreaks are urgently needed. The SARS-CoV-2 spike receptor-binding domain (RBD) comprises two regions, the core-RBD and the receptor-binding motif (RBM); the former is structurally conserved between SARS-CoV-2 and SARS-CoV. Here, in order to elicit humoral responses to the more conserved core-RBD, we introduced N-linked glycans onto RBM surfaces of the SARS-CoV-2 RBD and used them as immunogens in a mouse model. We found that glycan addition elicited higher proportions of the core-RBD-specific germinal center (GC) B cells and antibody responses, thereby manifesting significant neutralizing activity for SARS-CoV, SARS-CoV-2, and the bat WIV1-CoV. These results have implications for the design of SARS-like virus vaccines.

Author Correction: The oxidative inactivation of FeFe hydrogenase reveals the flexibility of the H-cluster.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Crowding in Cellular Environments at an Atomistic Level from Computer Simulations.

The effects of crowding in biological environments on biomolecular structure, dynamics, and function remain not well understood. Computer simulations of atomistic models of concentrated peptide and protein systems at different levels of complexity are beginning to provide new insights. Crowding, weak interactions with other macromolecules and metabolites, and altered solvent properties within cellular environments appear to remodel the energy landscape of peptides and proteins in significant ways including the possibility of native state destabilization. Crowding is also seen to affect dynamic properties, both conformational dynamics and diffusional properties of macromolecules. Recent simulations that address these questions are reviewed here and discussed in the context of relevant experiments.

Slow-Down in Diffusion in Crowded Protein Solutions Correlates with Transient Cluster Formation.

For a long time, the effect of a crowded cellular environment on protein dynamics has been largely ignored. Recent experiments indicate that proteins diffuse more slowly in a living cell than in a diluted solution, and further studies suggest that the diffusion depends on the local surroundings. Here, detailed insight into how diffusion depends on protein-protein contacts is presented based on extensive all-atom molecular dynamics simulations of concentrated villin headpiece solutions. After force field adjustments in the form of increased protein-water interactions to reproduce experimental data, translational and rotational diffusion was analyzed in detail. Although internal protein dynamics remained largely unaltered, rotational diffusion was found to slow down more significantly than translational diffusion as the protein concentration increased. The decrease in diffusion is interpreted in terms of a transient formation of protein clusters. These clusters persist on sub-microsecond time scales and follow distributions that increasingly shift toward larger cluster size with increasing protein concentrations. Weighting diffusion coefficients estimated for different clusters extracted from the simulations with the distribution of clusters largely reproduces the overall observed diffusion rates, suggesting that transient cluster formation is a primary cause for a slow-down in diffusion upon crowding with other proteins.

Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin.

The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.

The oxidative inactivation of FeFe hydrogenase reveals the flexibility of the H-cluster.

Nature is a valuable source of inspiration in the design of catalysts, and various approaches are used to elucidate the mechanism of hydrogenases, the enzymes that oxidize or produce H2. In FeFe hydrogenases, H2 oxidation occurs at the H-cluster, and catalysis involves H2 binding on the vacant coordination site of an iron centre. Here, we show that the reversible oxidative inactivation of this enzyme results from the binding of H2 to coordination positions that are normally blocked by intrinsic CO ligands. This flexibility of the coordination sphere around the reactive iron centre confers on the enzyme the ability to avoid harmful reactions under oxidizing conditions, including exposure to O2. The versatile chemistry of the diiron cluster in the natural system might inspire the design of novel synthetic catalysts for H2 oxidation.