Search

This computational study addresses the question of how large membrane-bound proteins of electron transport chains facilitate fast vector-based charge transport. We study electron transfer reactions following ultrafast initial charge separation induced by absorption of light by P800primary pair and leading to the electron localization at the A0cofactor. Two subsequent, much slower reactions, electron transfer to the iron–sulfur cluster Fxand reduction of the menaquinone (MQ) cofactor, are studied by combining molecular dynamics simulations, electronic structure calculations, and theoretical modeling. The low value of the electronic coupling between A0and Fxbrings this reaction to the microsecond time scale even at the zero activation barrier. In contrast, A0-MQ electron transfer occurs on a subnanosecond time scale and might become the preferred route for charge transport. We elucidate mechanistic properties of the protein medium allowing fast, vectorial charge transfer. The electric field is high and inhomogeneous inside the protein and is coupled to high polarizabilities of cofactors to significantly lower the reaction barrier. The A0-MQ separation puts this reaction at the edge between the plateau characterizing the reaction dynamical control and exponential falloff due to electronic tunneling. A strong separation in relaxation times of the medium dynamics for the forward and backward reactions promotes vectorial charge transfer.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein–nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9–2.5 Å) resolution. Three published maps were selected as targets: Escherichia colibeta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge-driven or artificial intelligence-driven and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary-structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models that contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly while refining them against the density maps. We introduce such an adaptive decision-making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two prototypical protein systems, adenylate kinase and carbon monoxide dehydrogenase. For these test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2–3 Å map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that are inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from the chimpanzee adenovirus. The overhead for decision making remains low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install R-MDFF on different computing environments, from local Linux-based workstations to high-performance computing environments.

Cytochrome c2(cyt. c2) is a major element in electron transfer between redox proteins in bioenergetic membranes. While the interaction between cyt. c2and anionic lipids abundant in bioenergetic membranes has been reported, their effect on the shuttling activity of cyt. c2remains elusive. Here, the effect of anionic lipids on the interaction and binding of cyt. c2to the cytochrome bc1complex (bc1) is investigated using a combination of molecular dynamics (MD) and Brownian dynamics (BD) simulations. MD is used to generate thermally accessible conformations of cyt. c2and membrane-embedded bc1, which were subsequently used in multireplica BD simulations of diffusion of cyt. c2from solution to bc1, in the presence of various lipids. We show that, counterintuitively, anionic lipids facilitate association of cyt. c2with bc1by localizing its diffusion to the membrane surface. The observed lipid-mediated bc1association is further enhanced by the oxidized state of cyt. c2, in line with its physiological function. This lipid-mediated enhancement is salinity-dependent, and anionic lipids can disrupt cyt. c2–bc1interaction at nonphysiological salt levels. Our data highlight the importance of the redox state of cyt. c2, the lipid composition of the chromatophore membrane, and the salinity of the chromatophore in regulating the efficiency of the electron shuttling process mediated by cyt. c2. The conclusions can be extrapolated to mitochondrial systems and processes, or any bioenergetic membrane, given the structural similarity between cyt. c2and bc1and their mitochondrial counterparts.

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As with any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the glycosylated spike protein, and identify the diversity of low-energy conformations in the vicinity of its open (so called 1RBD-up) state. The resulting atomic refinement reveal global and local molecular rearrangements that cannot be inferred from an average 1RBD-up cryo-EM model. Here we report varied degrees of “openness” in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these low-energy conformations. Five out of six binding sites we analyzed, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This apo-to-holo switching is reminiscent of a conformational preequilibrium. We found only one binding site, namely that of AB-C135 remains in apo state within all the sampled free energy-minimizing models, suggesting an induced fit mechanism for the docking of this antibody to the spike.

Predictive modeling of macromolecular recognition and protein-protein complementarity represents one of the cornerstones of biophysical sciences. However, such models are often hindered by the combinatorial complexity of interactions at the molecular interfaces. Exemplary of this problem is peptide presentation by the highly polymorphic major histocompatibility complex class I (MHC-I) molecule, a principal component of immune recognition. We developed human leukocyte antigen (HLA)-Inception, a deep biophysical convolutional neural network, which integrates molecular electrostatics to capture non-bonded interactions for predicting peptide binding motifs across 5,821 MHC-I alleles. These predictions of generated motifs correlate strongly with experimental peptide binding and presentation data. Beyond molecular interactions, the study demonstrates the application of predicted motifs in analyzing MHC-I allele associations with HIV disease progression and patient response to immune checkpoint inhibitors. A record of this paper’s transparent peer review process is included in the supplemental information.

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1?Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

Despite significant advances in resolution, the potential for cryo-electron microscopy (EM) to be used in determining the structures of protein–drug complexes remains unrealized. Determination of accurate structures and coordination of bound ligands necessitates simultaneous fitting of the models into the density envelopes, exhaustive sampling of the ligand geometries, and, most importantly, concomitant rearrangements in the side chains to optimize the binding energy changes. In this article, we present a flexible-fitting pipeline where molecular dynamics flexible fitting (MDFF) is used to refine structures of protein–ligand complexes from 3 to 5 Å electron density data. Enhanced sampling is employed to explore the binding pocket rearrangements. To provide a model that can accurately describe the conformational dynamics of the chemically diverse set of small-molecule drugs inside MDFF, we use QM/MM and neural-network potential (NNP)/MM models of protein–ligand complexes, where the ligand is represented using the QM or NNP model, and the protein is represented using established molecular mechanical force fields (e.g., CHARMM). This pipeline offers structures commensurate to or better than recently submitted high-resolution cryo-EM or X-ray models, even when given medium to low-resolution data as input. The use of the NNPs makes the algorithm more robust to the choice of search models, offering a radius of convergence of 6.5 Å for ligand structure determination. The quality of the predicted structures was also judged by density functional theory calculations of ligand strain energy. This strain potential energy is found to systematically decrease with better fitting to density and improved ligand coordination, indicating correct binding interactions. A computationally inexpensive protocol for computing strain energy is reported as part of the model analysis protocol that monitors both the ligand fit as well as model quality.

The mitochondrial respiratory chain, formed by five protein complexes, utilizes energy from catabolic processes to synthesize ATP. Complex I, the first and the largest protein complex of the chain, harvests electrons from NADH to reduce quinone, while pumping protons across the mitochondrial membrane. Detailed knowledge of the working principle of such coupled charge-transfer processes remains, however, fragmentary due to bottlenecks in understanding redox-driven conformational transitions and their interplay with the hydrated proton pathways. Complex I from Thermus thermophilusencases 16 subunits with nine iron–sulfur clusters, reduced by electrons from NADH. Here, employing the latest crystal structure of T. thermophiluscomplex I, we have used microsecond-scale molecular dynamics simulations to study the chemo-mechanical coupling between redox changes of the iron–sulfur clusters and conformational transitions across complex I. First, we identify the redox switches within complex I, which allosterically couple the dynamics of the quinone binding pocket to the site of NADH reduction. Second, our free-energy calculations reveal that the affinity of the quinone, specifically menaquinone, for the binding-site is higher than that of its reduced, menaquinol form—a design essential for menaquinol release. Remarkably, the barriers to diffusive menaquinone dynamics are lesser than that of the more ubiquitous ubiquinone, and the naphthoquinone headgroup of the former furnishes stronger binding interactions with the pocket, favoring menaquinone for charge transport in T. thermophilus. Our computations are consistent with experimentally validated mutations and hierarchize the key residues into three functional classes, identifying new mutation targets. Third, long-range hydrogen-bond networks connecting the quinone-binding site to the transmembrane subunits are found to be responsible for proton pumping. Put together, the simulations reveal the molecular design principles linking redox reactions to quinone turnover to proton translocation in complex I.

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Small, diffusible redox proteins play an essential role in electron transfer (ET) in respiration and photosynthesis, sustaining life on Earth by shuttling electrons between membrane-bound complexes via finely tuned and reversible interactions. Ensemble kinetic studies show transient ET complexes form in two distinct stages: an “encounter” complex largely mediated by electrostatic interactions, which subsequently, through subtle reorganization of the binding interface, forms a “productive” ET complex stabilized by additional hydrophobic interactions around the redox-active cofactors. Here, using single-molecule force spectroscopy (SMFS) we dissected the transient ET complexes formed between the photosynthetic reaction center-light harvesting complex 1 (RC-LH1) of Rhodobacter sphaeroidesand its native electron donor cytochrome c2(cyt c2). Importantly, SMFS resolves the distribution of interaction forces into low (∼150 pN) and high (∼330 pN) components, with the former more susceptible to salt concentration and to alteration of key charged residues on the RC. Thus, the low force component is suggested to reflect the contribution of electrostatic interactions in forming the initial encounter complex, whereas the high force component reflects the additional stabilization provided by hydrophobic interactions to the productive ET complex. Employing molecular dynamics simulations, we resolve five intermediate states that comprise the encounter, productive ET and leaving complexes, predicting a weak interaction between cyt c2and the LH1 ring near the RC-L subunit that could lie along the exit path for oxidized cyt c2. The multimodal nature of the interactions of ET complexes captured here may have wider implications for ET in all domains of life.

Transients in GPS time series can occur due to post-seismic deformation, seasonal hydrological loads, sea-level changes, flood and drought conditions, excessive groundwater withdrawal and recharge, etc. We report two new cases where the application of external loading, namely, earthquake loading and surface loading due to impoundment of hydroelectric reservoir, probably altered the local hydrological conditions to cause anomalous transients in the surface displacement. In the first case, moderate shaking due to the 2015 Gorkha earthquake at Patna (Bihar, India) caused transients in ground deformation in the following 50–60 days of the earthquake which are recorded by a continuous GPS site at Patna. In the second case, impoundment of the Tehri reservoir and its seasonal variations in the Garhwal Himalaya probably altered the local hydrological conditions which is causing anomalous biannual cyclic deformation at a site KUNR, near the reservoir.

Macromolecular assemblies often display a hierarchical organization of macromolecules or their subassemblies. To model this, we have formulated a space warping method that enables capturing overall macromolecular structure and dynamics via a set of coarse-grained order parameters (OPs). This article is the first of two describing the construction and computational implementation of an additional class of OPs that has built into them the hierarchical architecture of macromolecular assemblies. To accomplish this, first, the system is divided into subsystems, each of which is described via a representative set of OPs. Then, a global set of variables is constructed from these subsystem-centered OPs to capture overall system organization. Dynamical properties of the resulting OPs are compared to those of our previous nonhierarchical ones, and implied conceptual and computational advantages are discussed for a 100 ns, 2 million atom solvated human Papillomavirus-like particle simulation. In the second article, the hierarchical OPs are shown to enable a multiscale analysis that starts with the N-atom Liouville equation and yields rigorous Langevin equations of stochastic OP dynamics. The latter is demonstrated via a force-field-based simulation algorithm that probes key structural transition pathways, simultaneously accounting for all-atom details and overall structure.

The mechanism of rotatory catalysis in ATP-hydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V1-type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using X-ray crystallography, a new intermediate or “dwell” is identified, which enables the release of an inorganic phosphate (or Pi) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V1-ATPase to be the rate-limiting step of motor action.

Research is going on to synthesize materials for acting as sensors of metal ions and are also published in different journals but the cause for sensing has not been clearly explained. In this paper, the reason for selective sensing of Cu2+ions has been explained. A novel Schiff base fluorescent probe 3-[(2-hydroxy-naphthalen-1-ylmethylene)-hydrazono]-butan-2-one oxime (NPTH) was designed and synthesized to use as Cu2+sensor. From our recent experiments, it has been found that only Cu2+among different metal ions has been sensed by the synthesized ligand. The reason for selection and sensing of Cu2+by the ligand was established by different spectroscopic techniques. The detection limit of NPTH was calculated as low as 4.11 ​× ​10-4M for Cu2+. Therefore, these results indicate that sensor NPTH has great prospective to detect Cu2+ion in environmental analysis systems. Density functional theory (DFT) calculations have been done to ascertain the ground state geometry of NPTH.

Polymorphisms in MHC-I protein sequences across human populations significantly affect viral peptide binding capacity, and thus alter T cell immunity to infection. In the present study, we assess the relationship between observed SARS-CoV-2 population mortality and the predicted viral binding capacities of 52 common MHC-I alleles. Potential SARS-CoV-2 MHC-I peptides are identified using a consensus MHC-I binding and presentation prediction algorithm called EnsembleMHC. Starting with nearly 3.5 million candidates, we resolve a few hundred highly probable MHC-I peptides. By weighing individual MHC allele-specific SARS-CoV-2 binding capacity with population frequency in 23 countries, we discover a strong inverse correlation between predicted population SARS-CoV-2 peptide binding capacity and mortality rate. Our computations reveal that peptides derived from the structural proteins of the virus produce a stronger association with observed mortality rate, highlighting the importance of S, N, M, and E proteins in driving productive immune responses.