Search

The function of proteins is linked to their conformations that can be resolved with several high-resolution methods. However, only a few methods can provide the temporal order of intermediates and conformational changes, with each having its limitations. Here, we combine pulsed electron-electron double resonance spectroscopy with a microsecond freeze-hyperquenching setup to achieve spatiotemporal resolution in the angstrom range and lower microsecond time scale. We show that the conformational change of the C α -helix in the cyclic nucleotide-binding domain of the Mesorhizobium loti potassium channel occurs within about 150 μs and can be resolved with angstrom precision. Thus, this approach holds great promise for obtaining 4D landscapes of conformational changes in biomolecules.

F 1 -ATPase is a highly efficient molecular motor that can synthesize ATP driven by a mechanical torque. Its ability to function reversibly in either direction requires tight mechanochemical coupling between the catalytic domain and the rotating central shaft, as well as temporal control of substrate binding and product release. Despite great efforts and significant progress, the molecular details of this synchronized and fine-tuned energy conversion mechanism are not fully understood. Here, we use extensive molecular dynamics simulations to reconcile recent single-molecule experiments with structural data and provide a consistent thermodynamic, kinetic and mechanistic description of the main rotary substep in the synthetic cycle of mammalian ATP synthase. The calculated free energy profiles capture a discrete pattern in the rotation of the central γ-shaft, with a metastable intermediate located-consistently with recent experimental findings-at 70° relative to the X-ray position. We identify this rotary step as the ATP-dependent substep, and find that the associated free energy input supports the mechanism involving concurrent nucleotide binding and release. During the main substep, our simulations show no significant opening of the ATP-bound β subunit; instead, we observe that mechanical energy is transmitted to its nucleotide binding site, thus lowering the affinity for ATP. Simultaneously, the empty subunit assumes a conformation that enables the enzyme to harness the free energy of ADP binding to drive ATP release. Finally, we show that ligand exchange is regulated by a checkpoint mechanism, an apparent prerequisite for high efficiency in protein nanomotors.

There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A. , 2009 , 106 , 440 - 444 ) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large 1 H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pK a upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical 1 H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc. , 2014 , 136 , 3542 - 3552 ). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure.

G protein coupled receptors (GPCRs) transmit extracellular signals into the cell by binding and activating different intracellular signaling proteins, such as G proteins (Gαβγ, families Gi, Gs, Gq, G12/13) or arrestins. To address the issue of Gs vs Gi coupling specificity, we carried out molecular dynamics simulations of lipid-embedded active β2-adrenoceptor (β2AR*) in complex with C-terminal peptides derived from the key interaction site of Gα (GαCT) as surrogate of Gαβγ. We find that GiαCT and GsαCT exploit distinct cytoplasmic receptor conformations that coexist in the uncomplexed β2AR*. The slim GiαCT stabilizes a β2AR* conformation, not accessible to the bulkier GsαCT, which requires a larger TM6 outward tilt for binding. Our results suggest that the TM6 conformational heterogeneity regulates the catalytic activity of β2AR* toward Gi or Gs.

F1-ATPase (F1) is the catalytic portion of ATP synthase, a rotary motor protein that couples proton gradients to ATP synthesis. Driven by a proton flux, the F1 asymmetric γ subunit undergoes a stepwise rotation inside the α3β3 headpiece and causes the β subunits' binding sites to cycle between states of different affinity for nucleotides. These concerted transitions drive the synthesis of ATP from ADP and phosphate. Here, we study the coupling between the mechanical progression of γ and the conformations of α3β3. Using molecular dynamics simulations, we show that the nucleotide-free β subunit, initially in the open, low-affinity state, undergoes a spontaneous closing transition to the half-open state in response to the γ rotation in the synthesis direction. We estimate the kinetics of this spontaneous conformational change and analyze its mechanism and driving forces. By computing free energy profiles, we find that the isolated empty β subunit preferentially adopts the half-open conformation and that the transition to this conformation from the fully open state is accompanied by well-defined changes in the structure and interactions of the active site region. These results suggest that ADP binding to F1 occurs via conformational selection and is preceded by the transition of the active site to the half-open conformation, driven by the intrinsic elasticity of β. Our results also indicate that opening of the nucleotide-free β during hydrolysis is not spontaneous, as previously assumed. Rather, the fully open conformation observed in the F1 X-ray structure is enforced sterically by the γ subunit whose orientation is stabilized by interactions with the two other β subunits in the completely closed state. This finding supports the notion that γ acts by coupling the extreme conformational states of β subunits within the α3β3 hexamer and therefore is responsible for high efficiency of the coordinated catalysis.

The molecular mechanism of urea-induced protein denaturation is not yet fully understood. Mainly two opposing mechanisms are controversially discussed, according to which either hydrophobic, or polar interactions are the dominant driving force. To resolve this question, we have investigated the interactions between urea and all 20 amino acids by comprehensive molecular dynamics simulations of 22 tripeptides. Calculation of atomic contact frequencies between the amino acids and solvent molecules revealed a clear profile of solvation preferences by either water or urea. Almost all amino acids showed preference for contacts with urea molecules, whereas charged and polar amino acids were found to have slight preferences for contact with water molecules. Particularly strong preference for contacts to urea were seen for aromatic and apolar side-chains, as well as for the protein backbone of all amino acids. Further, protein-urea hydrogen bonds were found to be significantly weaker than protein-water or water-water hydrogen bonds. Our results suggest that hydrophobic interactions are the dominant driving force, while hydrogen bonds between urea and the protein backbone contribute markedly to the overall energetics by avoiding unfavorable unsatisfied hydrogen bond sites on the backbone. In summary, we suggest a combined mechanism that unifies the two current and seemingly opposing views.

Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy-lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited-state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier. The origin of this effect lies in the temporal electrostatic stabilization of dipole in the charge-transfer state in DNA.