Accelerated Molecular Dynamics and AlphaFold Uncover a Missing Conformational State of Transporter Protein OxlT.

Transporter proteins change their conformations to carry their substrate across the cell membrane. The conformational dynamics is vital to understanding the transport function. We have studied the oxalate transporter (OxlT), an oxalate:formate antiporter from Oxalobacter formigenes, significant in avoiding kidney stone formation. The atomic structure of OxlT has been recently solved in the outward-open and occluded states. However, the inward-open conformation is still missing, hindering a complete understanding of the transporter. Here, we performed a Gaussian accelerated molecular dynamics simulation to sample the extensive conformational space of OxlT and successfully predicted the inward-open conformation where cytoplasmic substrate formate binding was preferred over oxalate binding. We also identified critical interactions for the inward-open conformation. The results were complemented by an AlphaFold2 structure prediction. Although AlphaFold2 solely predicted OxlT in the outward-open conformation, mutation of the identified critical residues made it partly predict the inward-open conformation, identifying possible state-shifting mutations.

Electrostatic Ratchet for Successive Peptide Synthesis in Nonribosomal Molecular Machine RimK.

A nonribosomal peptide-synthesizing molecular machine, RimK, adds l-glutamic acids to the C-terminus of ribosomal protein S6 (RpsF) in vivo and synthesizes poly-α-glutamates in vitro. However, the mechanism of the successive glutamate addition, which is fueled by ATP, remains unclear. Here, we investigate the successive peptide-synthesizing mechanism of RimK via the molecular dynamics (MD) simulation of glutamate binding. We first show that RimK adopts three stable structural states with respect to the ATP-binding loop and the triphosphate chain of the bound ATP. We then show that a glutamate in solution preferentially binds to a positively charged belt-like region of RimK and the bound glutamate exhibits Brownian motion along the belt. The binding-energy landscape shows that the open-to-closed transition of the ATP-binding loop and the bent-to-straight transition of the triphosphate chain of ATP can function as an electrostatic ratchet that guides the bound glutamate to the active site. We then show the binding site of the second glutamate, which allows us to infer the ligation mechanism. Consistent with MD results, the crystal structure of RimK we obtained in the presence of RpsF presents an electron density that is presumed to correspond to the C-terminus of RpsF. We finally propose a mechanism for the successive peptide synthesis by RimK and discuss its similarity to other molecular machines.

Coupling of Redox and Structural States in Cytochrome P450 Reductase Studied by Molecular Dynamics Simulation.

Cytochrome P450 reductase (CPR) is the key protein that regulates the electron transfer from NADPH to various heme-containing monooxygenases. CPR has two flavin-containing domains: one with flavin adenine dinucleotide (FAD), called FAD domain, and the other with flavin mononucleotide (FMN), called FMN domain. It is considered that the electron transfer occurs via FAD and FMN (NADPH → FAD → FMN → monooxygenase) and is regulated by an interdomain open-close motion. It is generally thought that the structural state is coupled with the redox state, which, however, has not yet been firmly established. In this report, we studied the coupling of the redox and the structural states by full-scale molecular dynamics (MD) simulation of CPR (total 86.4 µs). Our MD result showed that while CPR predominantly adopts the closed state both in the oxidized and reduced states, it exhibits a tendency to open in the reduced state. We also found a correlation between the FAD-FMN distance and the predicted FMN-monooxygenase distance, which is embedded in the equilibrium thermal fluctuation of CPR. Based on these results, a physical mechanism for the electron transfer by CPR is discussed.

Hydrophobic Surface Enhances Electrostatic Interaction in Water.

A high dielectric constant is one of the peculiar properties of liquid water, indicating that the electrostatic interaction between charged substances is largely reduced in water. We show by molecular dynamics simulation that the dielectric constant of water is decreased near the hydrophobic surface. We further show that the decrease in the dielectric constant is due to both the decreased water density and the reduced water dipole correlation in the direction perpendicular to the surface. We finally demonstrate that electrostatic interaction in water is actually strengthened near the hydrophobic surface.

Long-range coupling between ATP-binding and lever-arm regions in myosin via dielectric allostery.

A protein molecule is a dielectric substance, so the binding of a ligand is expected to induce dielectric response in the protein molecule, considering that ligands are charged or polar in general. We previously reported that binding of adenosine triphosphate (ATP) to molecular motor myosin actually induces such a dielectric response in myosin due to the net negative charge of ATP. By this dielectric response, referred to as "dielectric allostery," spatially separated two regions in myosin, the ATP-binding region and the actin-binding region, are allosterically coupled. In this study, from the statistically stringent analyses of the extensive molecular dynamics simulation data obtained in the ATP-free and the ATP-bound states, we show that there exists the dielectric allostery that transmits the signal of ATP binding toward the distant lever-arm region. The ATP-binding-induced electrostatic potential change observed on the surface of the main domain induced a movement of the converter subdomain from which the lever arm extends. The dielectric response was found to be caused by an underlying large-scale concerted rearrangement of the electrostatic bond network, in which highly conserved charged/polar residues are involved. Our study suggests the importance of the dielectric property for molecular machines in exerting their function.

Electrostatic balance between global repulsion and local attraction in reentrant polymerization of actin.

Actin polymerization depends on the salt concentration, exhibiting a reentrant behavior: the polymerization is promoted by increasing KCl concentration up to 100 mM, and then depressed by further increase above 100 mM. We here investigated the physical mechanism of this reentrant behavior by calculating the polymerization energy, defined by the electrostatic energy change upon binding of an actin subunit to a filament, using an implicit solvent model based on the Poisson-Boltzmann (PB) equation. We found that the polymerization energy as a function of the salt concentration shows a non-monotonic reentrant-like behavior, with the minimum at about 100 mM (1:1 salt). By separately examining the salt concentration effect on the global electrostatic repulsion between the like-charged subunits and that on the local electrostatic attraction between the inter-subunit ionic-bond-forming residues in the filament, we clarified that the reentrant behavior is caused by the change in the balance between the two opposing electrostatic interactions. Our study showed that the non-specific nature of counterions, as described in the mean-field theory, plays an important role in the actin polymerization. We also discussed the endothermic nature of the actin polymerization and mentioned the effect of ATP hydrolysis on the G-F transformation, indicating that the electrostatic interaction is widely and intricately involved in the actin dynamics.

Over-Destabilization of Protein-Protein Interaction in Generalized Born Model and Utility of Energy Density Integration Cutoff.

The generalize Born (GB) model is frequently used in MD simulations of biomolecular systems in aqueous solution. The GB model is usually based on the so-called Coulomb field approximation (CFA) for the energy density integration. In this study, we report that the GB model with CFA overdestabilizes the long-range electrostatic attraction between oppositely charged molecules (ionic bond forming two-helix system and kinesin-tubulin system) when the energy density integration cutoff, rmax, which is used to calculate the Born energy, is set to a large value. We show that employing large rmax, which is usually expected to make simulation results more accurate, worsens the accuracy so that the attraction is changed into repulsion. It is demonstrated that the overdestabilization is caused by the overestimation of the desolvation penalty upon binding that originates from CFA. We point out that the overdestabilization can be corrected by employing a relatively small cutoff (rmax = 10-15 Å), affirming that the GB models, even with CFA, can be used as a powerful tool to theoretically study the protein-protein interaction, particularly on its dynamical aspect, such as binding and unbinding.

Dielectric Allostery of Protein: Response of Myosin to ATP Binding.

Protein uses allostery to execute biological function. The physical mechanism underlying the allostery has long been studied, with the focus on the mechanical response by ligand binding. Here, we highlight the electrostatic response, presenting an idea of "dielectric allostery". We conducted molecular dynamics simulations of myosin, a motor protein with allostery, and analyzed the response to ATP binding which is a crucial step in force-generating function, forcing myosin to unbind from the actin filament. We found that the net negative charge of ATP causes a large-scale, anisotropic dielectric response in myosin, altering the electrostatic potential in the distant actin-binding region and accordingly retracting a positively charged actin-binding loop. A large-scale rearrangement of electrostatic bond network was found to occur upon ATP binding. Since proteins are dielectric and ligands are charged/polar in general, the dielectric allostery might underlie a wide spectrum of functions by proteins.

Piezoelectric allostery of protein.

Allostery is indispensable for a protein to work, where a locally applied stimulus is transmitted to a distant part of the molecule. While the allostery due to chemical stimuli such as ligand binding has long been studied, the growing interest in mechanobiology prompts the study of the mechanically stimulated allostery, the physical mechanism of which has not been established. By molecular dynamics simulation of a motor protein myosin, we found that a locally applied mechanical stimulus induces electrostatic potential change at distant regions, just like the piezoelectricity. This novel allosteric mechanism, "piezoelectric allostery", should be of particularly high value for mechanosensor/transducer proteins.

Intrinsic disorder accelerates dissociation rather than association.

The intrinsically disordered protein (IDP) has distinct properties both physically and biologically: it often becomes folded when binding to the target and is frequently involved in signal transduction. The physical property seems to be compatible with the biological property where fast association and dissociation between IDP and the target are required. While fast association has been well studied, fueled by the fly-casting mechanism, the dissociation kinetics has received less attention. We here study how the intrinsic disorder affects the dissociation kinetics, as well as the association kinetics, paying attention to the interaction strength at the binding site (i.e., the quality of the "fly lure"). Coarse-grained molecular dynamics simulation of the pKID-KIX system, a well-studied IDP system, shows that the association rate becomes larger as the disorder-inducing flexibility that was imparted to the model is increased, but the acceleration is marginal and turns into deceleration as the quality of the fly lure is worsened. In contrast, the dissociation rate is greatly enhanced as the disorder is increased, indicating that intrinsic disorder serves for rapid signal switching more effectively through dissociation than association. Proteins 2016; 84:1124-1133. © 2016 Wiley Periodicals, Inc.