Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: A Symmetry-Based Perspective.

We report on a new method for the characterization of local structures in proteins based on extensive molecular dynamics (MD) simulations, here, 1 µs in length. The N-H bond of the Rho GTPase binding domain of plexin-B1 (RBD) serves as a probe and the potential, u(MD), which restricts its internal motion, as a qualifier of the local dynamic structure. u(MD) is derived from the MD trajectory as a function of the polar angles, (θ, φ), which specify the N-H orientation in the protein. u(MD) is statistical in character yielding empirical descriptions. To establish more insightful methodical descriptions, we develop a comprehensive method which approximates u(MD) by combinations of analytical Wigner functions that belong to the D2h point group. These combinations, called u(simulated), make it possible to gain a new perspective of local dynamic structures in proteins based on explicit potentials/free energy surfaces and associated probability densities, entropy, and ordering. A simpler method was developed previously using 100 ns MD simulations. In that case, the traditional "perpendicular N-H ordering" setting centered at Cα-Cα with (θ, φ) = (90, 90) and generally, featuring positive φ, prevailed. u(MD) derived from 1 µs MD simulations is considerably more complex requiring substantial model enhancement. The enhanced method applies to the well-structured sections of the RBD. It only applies partly to its loops where u(MD) extends into the negative-φ region where we detect nonperpendicular N-H ordering. This arrangement requires devising new reference structures and making substantial algorithmic changes, to be performed in future work. Here, we focus on developing the comprehensive method and using it to investigate perpendicular ordering settings. We find that secondary structures (loops) exhibit varying (virtually invariant) potentials with Ag, B2u, and B1u (Ag and B2u) D2h symmetry. Application to RBD dimerization and RBD binding to the GTPase Rac1 is described in the subsequent article. Applications to other probes, proteins, and biological functions, based on explicit local potentials, probability densities, entropy, and ordering, are possible.

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: 2. The Role of Symmetry in GTPase Binding and Dimer Formation.

The Rho GTPase binding domain of Plexin-B1 (RBD) prevails in solution as dimer. Under appropriate circumstances, it binds the small GTPase Rac1 to yield the complex RBD-Rac1. Here, we study RBD dimerization and complex formation from a symmetry-based perspective using data derived from 1 µs long MD simulations. The quantities investigated are the local potentials, u(MD), prevailing at the N-H sites of the protein. These potentials are statistical in character providing an empirical description of the local structure. To establish more methodical description, a method for approximating them by explicit functions, u(simulated), was developed in the preceding article in this journal issue. These functions are combinations of analytical Wigner functions, DL,K, belonging to the D2h point group. The D2h subgroups Ag and B2u are found to dominate u(simulated); the B1u subgroup contributes in some cases. The Ag (B2u) functions have axial or rhombic symmetry. For the first time, local potentials in proteins can be quantitatively characterized in terms of their strength (rhombicity) evaluated by axial Ag (rhombic Ag and B2u) contributions. Until now, the chain-segment [ß3-L3-ß4] and to some extent the α2-helix have been associated with GTPase binding. Here, we find that this process causes an increase (decrease) in the potential strength of ß3 and ß4 (the preceding L2 loop and the remote chain-segment [(α2-helix)-(α2/ß5-turn)-(ß5-strand)]), suggesting effects of counterbalancing and allostery. There is evidence for the L2 loop being associated with RBD-GTPase binding. Until now only the L4 loop has been associated with RBD dimerization. The latter process is found to cause an increase (decrease) in the potential strength and rhombicity of the L4 loop (the adjacent chain-segment [(α2-helix)-(α2/ß5-turn)-(ß5-strand)]), suggesting counterbalancing activity. On average, the RBD dimer features stronger local potentials than RBD-Rac1. The novel information inherent in these findings is mesoscopic in character. Prospects of interest include exploring relation to atomistic force-field parameters.

Microsecond MD Simulations of the Plexin-B1 RBD: 2. N-H Probability Densities and Conformational Entropy in Ligand-Free, Rac1-Bound, and Dimer RBD.

Orientational probability densities, Peq = exp(-u) (u, local potential), of bond-vectors in proteins provide information on structural flexibility. The related conformational entropy, Sk = -∫Peq(ln Peq)dΩ - ln ∫dΩ, provides the entropic contribution to the free energy of the physical/biological process studied. We have developed a new method for deriving Peq and Sk from MD simulations, using the N-H bond as probe. Recently we used it to study the dimerization of the Rho GTPase binding domain of Plexin-B1 (RBD). Here we use it to study RBD binding to the small GTPase Rac1. In both cases 1 µs MD simulations have been employed. The RBD has the ubiquitin fold with four mostly long loops. L3 is associated with GTPase binding, L4 with RBD dimerization, L2 participates in interdomain interactions, and L1 has not been associated with function. We find that RBD-Rac1 binding renders L1, L3, and L4 more rigid and the turns ß2/α1 and α2/ß5 more flexible. By comparison, RBD dimerization renders L4 more rigid, and the α-helices, the ß-strands, and L2 more flexible. The rigidity of L1 in RBDRAC is consistent with L1-L3 contacts seen in previous MD simulations. The analysis of the L3-loop reveals two states of distinct flexibility which we associate with involvement in slow conformational exchange processes differing in their rates. Overall, the N-H bonds make an unfavorable entropic contribution of (5.9 ± 0.9) kJ/mol to the free energy of RBD-Rac1 binding; they were found to make a favorably contribution of (-7.0 ± 0.7) kJ/mol to the free energy of RBD dimerization. In summary, the present study provides a new perspective on the impact of Rac1 binding and dimerization on the flexibility characteristics of the RBD. Further studies are stimulated by the results of this work.

Microsecond MD Simulations of the Plexin-B1 RBD: N-H Probability Density as Descriptor of Structural Dynamics, Dimerization-Related Conformational Entropy, and Transient Dimer Asymmetry.

Amide-bond equilibrium probability density, Peq = exp(-u) (u, local potential), and associated conformational entropy, Sk = -∫Peq (ln Peq) dΩ â”€ln ∫dΩ, are derived for the Rho GTPase binding domain of Plexin-B1 (RBD) as monomer and dimer from 1 µs MD simulations. The objective is to elucidate the effect of dimerization on the dynamic structure of the RBD. Dispersed (peaked) Peq functions indicate "flexibility" ("rigidity"; the respective concepts are used below in this context). The L1 and L3 loops are throughout highly flexible, the L2 loop and the secondary structure elements are generally rigid, and the L4 loop is flexible in the monomer and rigid in the dimer. Overall, many residues are more flexible in the dimer. These features, and their implications, are discussed. Unexpectedly, we find that monomer unit 1 of the dimer (in short, d1) is unusually flexible, whereas monomer unit 2 (in short, d2) is as rigid as the RBD monomer. This is revealed due to their engagement in slow-to-intermediate conformational exchange detected previously by 15N relaxation experiments. Such motions occur with rates on the order of 103-104 s-1; hence, they cannot be completely sampled over the course of 1 µs simulation. However, the extent to which rigid d2 is affected is small enough to enable physically relevant analysis. The entropy difference between d2 and the monomer yields an entropic contribution of -7 ± 0.7 kJ/mol to the free energy of RBD dimerization. In previous work aimed at similar objectives we used 50-100 ns MD simulations. Those results and the present result differ considerably. In summary, bond-vector Peq functions derived directly from long MD simulations are useful descriptors of protein structural dynamics and provide accurate conformational entropy. Within the scope of slow conformational exchange, they can be useful, even in the presence of incomplete sampling.

Evolutionary Effects on Bound Substrate p Ka in Dihydrofolate Reductase.

In the present study, we address the effect of active site structure and dynamics of different dihydrofolate reductase (DHFR) isoforms on the p Ka of the bound substrate 7,8-dihydrofolate, in an attempt to understand possible evolutionary trends. We apply a hybrid QM/MM free energy perturbation method to estimate the p Ka of the N5 position of the bound substrate. We observe a gradual increase in N5 basicity as we move from primitive to more evolved DHFR isoforms. Structural analysis of these isoforms reveals a gradual sequestering of water molecules from the active site in the more evolved enzymes, thereby modulating the local dielectric environment near the substrate. Furthermore, the present study reveals a clear correlation between active site hydration and the N5 p Ka of the substrate. We emphasize the role of the M20 loop in controlling the active site hydration level, via a preorganized active site with a more hydrophobic environment and reduced loop flexibility as evolution progresses from bacterial to the human enzyme.

Effect of Asp122 Mutation on the Hydride Transfer in E. coli DHFR Demonstrates the Goldilocks of Enzyme Flexibility.

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF) in the presence of NADPH. The key hydride transfer step in the reaction is facilitated by a combination of enzyme active site preorganization and correlated protein motions in the Michaelis-Menten (E:NADPH:DHF) complex. The present theoretical study employs mutagenesis to examine the relation between structural and functional properties of the enzyme. We mutate Asp122 in Escherichia coli DHFR, which is a conserved amino acid in the DHFR family. The consequent effect of the mutation on enzyme catalysis is examined from an energetic, structural and short-time dynamic perspective. Our investigations suggest that the structural and short-time dynamic perturbations caused by Asp122X mutations (X = Asn, Ser, Ala) are along the reaction coordinate and lower the rate of hydride transfer. Importantly, analysis of the correlated and principle component motions in the enzyme suggest that the mutation alters the coupled motions that are present in the wild-type enzyme. In the case of D122N and D122S, the mutations inhibit coupled motion, whereas in the case of D122A, the mutation enhances coupled motion, although all mutations result in similar rate reduction. These results emphasize a Goldilocks principle of enzyme flexibility, that is, enzymes should neither be too rigid nor too flexible.

Improved Sugar Puckering Profiles for Nicotinamide Ribonucleoside for Hybrid QM/MM Simulations.

The coenzyme nicotinamide adenine dinucleotide (NAD+) and its reduced form (NADH) play ubiquitous roles as oxidizing and reducing agents in nature. The binding, and possibly the chemical redox step, of NAD+/NADH may be influenced by the cofactor conformational distribution and, in particular, by the ribose puckering of its nicotinamide-ribonucleoside (NR) moiety. In many hybrid quantum mechanics-molecular mechanics (QM/MM) studies of NAD+/NADH dependent enzymes, the QM region is treated by semiempirical (SE) methods. Recent work suggests that SE methods do not adequately describe the ring puckering in sugar molecules. In the present work we adopt an efficient and practical strategy to correct for this deficiency for NAD+/NADH. We have implemented a cost-effective correction to a SE Hamiltonian by adding a correction potential, which is defined as the difference between an accurate benchmark density functional theory (DFT) potential energy surface (PES) and the SE PES. In practice, this is implemented via a B-spline interpolation scheme for the grid-based potential energy difference surface. We find that the puckering population distributions obtained from free energy QM(SE)/MM simulations are in good agreement with DFT and in fair accord with experimental results. The corrected PES should facilitate a more accurate description of the ribose puckering in the NAD+/NADH cofactor in simulations of biological systems.