Hyper-swivel head domain motions are required for complete mRNA-tRNA translocation and ribosome resetting.

Translocation of messenger RNA (mRNA) and transfer RNA (tRNA) substrates through the ribosome during protein synthesis, an exemplar of directional molecular movement in biology, entails a complex interplay of conformational, compositional, and chemical changes. The molecular determinants of early translocation steps have been investigated rigorously. However, the elements enabling the ribosome to complete translocation and reset for subsequent protein synthesis reactions remain poorly understood. Here, we have combined molecular simulations with single-molecule fluorescence resonance energy transfer imaging to gain insights into the rate-limiting events of the translocation mechanism. We find that diffusive motions of the ribosomal small subunit head domain to hyper-swivelled positions, governed by universally conserved rRNA, can maneuver the mRNA and tRNAs to their fully translocated positions. Subsequent engagement of peptidyl-tRNA and disengagement of deacyl-tRNA from mRNA, within their respective small subunit binding sites, facilitate the ribosome resetting mechanism after translocation has occurred to enable protein synthesis to resume.

Full-Length Computational Model of the SARS-CoV-2 Spike Protein and Its Implications for a Viral Membrane Fusion Mechanism.

The SARS-CoV-2 virus has now become one of the greatest causes of infectious death and morbidity since the 1918 flu pandemic. Substantial and unprecedented progress has been made in the elucidation of the viral infection process in a short time; however, our understanding of the structure-function dynamics of the spike protein during the membrane fusion process and viral uptake remains incomplete. Employing computational approaches, we use full-length structural models of the SARS-CoV-2 spike protein integrating Cryo-EM images and biophysical properties, which fill the gaps in our understanding. We propose a membrane fusion model incorporating structural transitions associated with the proteolytic processing of the spike protein, which initiates and regulates a series of events to facilitate membrane fusion and viral genome uptake. The membrane fusion mechanism highlights the notable role of the S1 subunit and eventual mature spike protein uptake through the host membrane. Our comprehensive view accounts for distinct neutralizing antibody binding effects targeting the spike protein and the enhanced infectivity of the SARS-CoV-2 variant.

Scaling molecular dynamics beyond 100,000 processor cores for large-scale biophysical simulations.

The growing interest in the complexity of biological interactions is continuously driving the need to increase system size in biophysical simulations, requiring not only powerful and advanced hardware but adaptable software that can accommodate a large number of atoms interacting through complex forcefields. To address this, we developed and implemented strategies in the GENESIS molecular dynamics package designed for large numbers of processors. Long-range electrostatic interactions were parallelized by minimizing the number of processes involved in communication. A novel algorithm was implemented for nonbonded interactions to increase single instruction multiple data (SIMD) performance, reducing memory usage for ultra large systems. Memory usage for neighbor searches in real-space nonbonded interactions was reduced by approximately 80%, leading to significant speedup. Using experimental data describing physical 3D chromatin interactions, we constructed the first atomistic model of an entire gene locus (GATA4). Taken together, these developments enabled the first billion-atom simulation of an intact biomolecular complex, achieving scaling to 65,000 processes (130,000 processor cores) with 1 ns/day performance. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Structural Transition and Antibody Binding of EBOV GP and ZIKV E Proteins from Pre-Fusion to Fusion-Initiation State.

Membrane fusion proteins are responsible for viral entry into host cells—a crucial first step in viral infection. These proteins undergo large conformational changes from pre-fusion to fusion-initiation structures, and, despite differences in viral genomes and disease etiology, many fusion proteins are arranged as trimers. Structural information for both pre-fusion and fusion-initiation states is critical for understanding virus neutralization by the host immune system. In the case of Ebola virus glycoprotein (EBOV GP) and Zika virus envelope protein (ZIKV E), pre-fusion state structures have been identified experimentally, but only partial structures of fusion-initiation states have been described. While the fusion-initiation structure is in an energetically unfavorable state that is difficult to solve experimentally, the existing structural information combined with computational approaches enabled the modeling of fusion-initiation state structures of both proteins. These structural models provide an improved understanding of four different neutralizing antibodies in the prevention of viral host entry.

Characterization of Conformational Ensembles of Protonated N-glycans in the Gas-Phase.

Ion mobility mass spectrometry (IM-MS) is a technique capable of investigating structural changes of biomolecules based on their collision cross section (CCS). Recent advances in IM-MS allow us to separate carbohydrate isomers with subtle conformational differences, but the relationship between CCS and atomic structure remains elusive. Here, we characterize conformational ensembles of gas-phase N-glycans under the electrospray ionization condition using molecular dynamics simulations with enhanced sampling. We show that the separation of CCSs between isomers reflects folding features of N-glycans, which are determined both by chemical compositions and protonation states. Providing a physicochemical basis of CCS for N-glycans helps not only to interpret IM-MS measurements but also to estimate CCSs of complex glycans.

Backbone Circularization Coupled with Optimization of Connecting Segment in Effectively Improving the Stability of Granulocyte-Colony Stimulating Factor.

Backbone circularization of protein is a powerful method to improve its structural stability. In this paper, we presumed that a tight connection leads to much higher stability. Therefore, we designed circularized variants of a granulocyte-colony stimulating factor (G-CSF) with a structurally optimized terminal connection. To estimate the appropriate length of the connection, we surveyed the Protein Data Bank to find local structures as a model for the connecting segment. We set the library of local structures composed of "helix-loop-helix," subsequently selected entries similar to the G-CSF terminus, and finally sorted the hit structures according to the loop length. Two, five, or nine loop residues were frequently observed; thus, three circularized variants (C163, C166, and C170) were constructed, prepared, and evaluated. All circularized variants demonstrated a higher thermal stability than linear G-CSF (L175). In particular, C166 that retained five connecting residues demonstrated apparent Tm values of 69.4 °C, which is 8.7 °C higher than that of the circularized variant with no truncation (C177), indicating that the optimization of the connecting segment is effective for enhancing the overall structural stability. C166 also showed higher proteolytic stability against both endoprotease and exopeptidase than L175. We anticipate that the present study will contribute to the improvement in the general design of circularized protein and development of G-CSF biobetters.

Mechanisms for Two-Step Proton Transfer Reactions in the Outward-Facing Form of MATE Transporter.

Bacterial pathogens or cancer cells can acquire multidrug resistance, which causes serious clinical problems. In cells with multidrug resistance, various drugs or antibiotics are extruded across the cell membrane by multidrug transporters. The multidrug and toxic compound extrusion (MATE) transporter is one of the five families of multidrug transporters. MATE from Pyrococcus furiosus uses H(+) to transport a substrate from the cytoplasm to the outside of a cell. Crystal structures of MATE from P. furiosus provide essential information on the relevant H(+)-binding sites (D41 and D184). Hybrid quantum mechanical/molecular mechanical simulations and continuum electrostatic calculations on the crystal structures predict that D41 is protonated in one structure (Straight) and, both D41 and D184 protonated in another (Bent). All-atom molecular dynamics simulations suggest a dynamic equilibrium between the protonation states of the two aspartic acids and that the protonation state affects hydration in the substrate binding cavity and lipid intrusion in the cleft between the N- and C-lobes. This hypothesis is examined in more detail by quantum mechanical/molecular mechanical calculations on snapshots taken from the molecular dynamics trajectories. We find the possibility of two proton transfer (PT) reactions in Straight: the 1st PT takes place between side-chains D41 and D184 through a transient formation of low-barrier hydrogen bonds and the 2nd through another H(+) from the headgroup of a lipid that intrudes into the cleft resulting in a doubly protonated (both D41 and D184) state. The 1st PT affects the local hydrogen bond network and hydration in the N-lobe cavity, which would impinge on the substrate-binding affinity. The 2nd PT would drive the conformational change from Straight to Bent. This model may be applicable to several prokaryotic H(+)-coupled MATE multidrug transporters with the relevant aspartic acids.

Mosaic of Water Orientation Structures at a Neutral Zwitterionic Lipid/Water Interface Revealed by Molecular Dynamics Simulations.

Ordering of water structures near the surface of biological membranes has been recently extensively studied using interface-selective techniques like vibrational sum frequency generation (VSFG) spectroscopy. The detailed structures of interface water have emerged for charged lipids, but those for neutral zwitterionic lipids remain obscure. We analyze an all-atom molecular dynamics (MD) trajectory of a hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer to characterize the orientation of interface waters in different chemical environments. The structure and dynamics of interfacial waters strongly depend on both their vertical position along the bilayer normal as well as vicinal lipid charged groups. Water orientation in the vicinity of phosphate groups is opposite to that around choline groups. The results are consistent with observed VSFG spectra and demonstrate that a mosaic of water orientation structures exists on the surface of a neutral zwitterionic phospholipid bilayer, reflecting rapid water exchange and the influence of local chemical environments.

Confident identification of isomeric N-glycan structures by combined ion mobility mass spectrometry and hydrophilic interaction liquid chromatography.

RATIONALE: A central issue in glycan mass analysis is the ambiguity of structural assignments due to the heterogeneity and complexity of glycan structures. Ion mobility mass spectrometry (IM-MS) has the potential to separate isomeric glycans depending on their unique collisional cross section especially when coupled with hydrophilic interaction liquid chromatography (HILIC). METHODS: Ten pyridylaminated biantennary N-glycans including isomeric structures were measured by electrospray ionization quadrupole-time-of-flight mass spectrometry with an ion mobility phase. We investigated which adduct ions would be suitable for good separation in the ion mobility phase. The differences in observed drift time of isomeric glycans were assessed by molecular dynamics (MD) simulations in vacuum. Connecting an HILIC system with IM-MS provided another, augmented separation mode. RESULTS: By selecting doubly protonated precursor ion species, we succeeded in separating a pair of isomeric glycans in the ion mobility phase with reasonable resolution. MD simulations of monogalactosylated glycan isomers indicate that the galactosylated Man α1-3 branch preferentially folds back to the core chitobiose portion to form a compact structure. IM-MS combined with HILIC resulted in even clearer separation of isomeric glycans within 15 min. CONCLUSIONS: A combination of IM-MS with an HILIC system is eminently suitable for the confident and rapid distinction of glycan structures within a defined mixture.

Effect of bisecting GlcNAc and core fucosylation on conformational properties of biantennary complex-type N-glycans in solution.

The introduction of bisecting GlcNAc and core fucosylation in N-glycans is essential for fine functional regulation of glycoproteins. In this paper, the effect of these modifications on the conformational properties of N-glycans is examined at the atomic level by performing replica-exchange molecular dynamics (REMD) simulations. We simulate four biantennary complex-type N-glycans, namely, unmodified, two single-substituted with either bisecting GlcNAc or core fucose, and disubstituted forms. By using REMD as an enhanced sampling technique, five distinct conformers in solution, each of which is characterized by its local orientation of the Manα1-6Man glycosidic linkage, are observed for all four N-glycans. The chemical modifications significantly change their conformational equilibria. The number of major conformers is reduced from five to two and from five to four upon the introduction of bisecting GlcNAc and core fucosylation, respectively. The population change is attributed to specific inter-residue hydrogen bonds, including water-mediated ones. The experimental NMR data, including nuclear Overhauser enhancement and scalar J-coupling constants, are well reproduced taking the multiple conformers into account. Our structural model supports the concept of "conformer selection", which emphasizes the conformational flexibility of N-glycans in protein-glycan interactions.

Conformational flexibility of N-glycans in solution studied by REMD simulations.

Protein-glycan recognition regulates a wide range of biological and pathogenic processes. Conformational diversity of glycans in solution is apparently incompatible with specific binding to their receptor proteins. One possibility is that among the different conformational states of a glycan, only one conformer is utilized for specific binding to a protein. However, the labile nature of glycans makes characterizing their conformational states a challenging issue. All-atom molecular dynamics (MD) simulations provide the atomic details of glycan structures in solution, but fairly extensive sampling is required for simulating the transitions between rotameric states. This difficulty limits application of conventional MD simulations to small fragments like di- and tri-saccharides. Replica-exchange molecular dynamics (REMD) simulation, with extensive sampling of structures in solution, provides a valuable way to identify a family of glycan conformers. This article reviews recent REMD simulations of glycans carried out by us or other research groups and provides new insights into the conformational equilibria of N-glycans and their alteration by chemical modification. We also emphasize the importance of statistical averaging over the multiple conformers of glycans for comparing simulation results with experimental observables. The results support the concept of "conformer selection" in protein-glycan recognition.

Screw motion regulates multiple functions of T4 phage protein gene product 5 during cell puncturing.

Bacteriophage T4 penetrates the outer membrane of Escherichia coli using a multifunctional device composed of a gene product 5 (gp5) protein trimer. We report that gp5 sequentially exerts distinct functions along the course of penetration stages induced by screw motion. A triple-stranded ß-helix of gp5 acts as a cell-puncturing drill bit to make a hole on the membrane and then send the lipids upward efficiently by strong charge interactions. The gp5 lysozyme domains, which degrade the peptidoglycan layer later, are shown to play novel roles to enlarge the hole and control the release of the ß-helix. The lysozyme active site is protected from lipid binding during the penetration and is exposed after the ß-helix release. Intrinsic multiple functions of gp5 are shown to be served in turn regulated by gradual change of interdomain interactions, which enables the initial infection process with single protein trimer by continuous screw motion. The results of lysozyme domain should be understood as the case where a single-function protein acquired multiple chemical functions through interplay with other domains in a multidomain protein.

DTA: dihedral transition analysis for characterization of the effects of large main-chain dihedral changes in proteins.

MOTIVATION: The biological function of proteins is associated with a variety of motions, ranging from global domain motion to local motion of side chain. We propose a method, dihedral transition analysis (DTA), to identify significant dihedral angle changes between two distinct protein conformations and for characterization of the effect of these transitions on both local and global conformation. RESULTS: Applying DTA to a comprehensive and non-redundant dataset of 459 high-resolution pairs of protein structures, we found that a dihedral transition occurs in 82% of proteins. Multiple dihedral transitions are shown to occur cooperatively along the sequence, which allows us to separate a polypeptide chain into fragments with and without transitions, namely transition fragments (TFs) and stable fragments (SFs), respectively. By characterizing the magnitude of TF conformational change and the effect of the transition on the neighboring fragments, flap and hinge motions are identified as typical motions. DTA is also useful to detect protein motions, subtle in RMSD but significant in terms of dihedral angle changes, such as the peptide-plane flip, the side-chain flip and path-preserving motions. We conclude that DTA is a useful tool to extract potential functional motions, some of which might have been missed using conventional methods for protein motion analysis. AVAILABILITY: http://dynamics.iam.u-tokyo.ac.jp/DTA/